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PREFACE TO RUSSIAN EDITION 


THE theory of electron—atom collisions is given in many monographs, text- 
books and reviews. The *‘Theory of Atomic Collisions”” by Mott and Massey 
is the most substantial among these publications (the second edition was 
issued in 1950). Massey’s articles “‘Theory of Atomic Collisions’? and 
“Excitation and Jonization of Atoms by Electron Impact’’ in Handbuch der 
Physik (1956) are also concerned with this subject. Yu. N. Demkov’s 
‘“Variatsionnye Printsipy v Teorn Stolknoveni’’ (1958) is devoted to 
variational methods. 

However, the theory of atomic collisions has recently developed very 
rapidly in directions which are either not considered in these books or are 
only briefly mentioned. For example, substantial progress has been made 
in the application of dispersion relations to atomic collisions. Important 
results have been obtained on the extremal properties of the scattering length 
derived by variational methods. Computers have facilitated extensive 
calculations of scattering cross-sections with allowance for such factors as 
the polarization of the target atom by incident electrons, the coupling 
between different channels and so on. Calculations of this kind are being 
carried out at a very intensive rate both in the USSR and abroad. It has 
been found that the factors just mentioned are in many cases very important. 
The semi-classical theory of electron—atom collisions, first developed by Niels 
Bohr in connection with the interaction of fast particles with matter, has 
been developed still further. 

The present book is largely devoted to recent developments in the theory 
of electron—atom collisions. It expounds both the general theory and its 
applications to specific problems. Chapters 1 and 2 are devoted to central 
force problems. The general theory of electron—atom collisions is discussed 
in Chapter 3 in the simplest case, i.e. the case of the hydrogen atom. 
Equations describing the collision of an electron with an arbitrary atom or 
ion are considered in Chapter 8. The remaining chapters are concerned with 
concrete calculations. 

The general theory given in Chapters 3 and 8 is valid at any energy of the 
incident electrons, and covers both excitation and ionization processes. 
However, in concrete applications of the general theory we shall confine our 
attention to the interaction of slow electrons with atoms, i.e. elastic scattering 
at energies of a few electron volts, and inelastic scattering near the excitation 
threshold. It is precisely this energy range which is of the greatest interest in 
connection with the study of plasma processes and many other phenomena. 
On the other hand, approximate methods based on perturbation theory are 
not valid in this energy range, and a more complete analysis becomes neces- 
sary. For obvious reasons, most attention will be given to the interaction of 
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electrons with the simplest systems, 1.c. the hydrogen atom, the helium atom 
and the helium ton. 

The account which we shall give of concrete calculations of collisions 
between electrons and the simple systems mentioned above is not intended 
to be exhaustive. We shall be mainly concerned with recent calculations in 
which the various complicating factors are taken into account. The results 
of some earlier papers will also be summarized. Ionization processes will 
not be considered. 

I wish to thank Academician V. A. Fock, Yu. N. Demkov and L. D. 
Faddeev, whose comments have contributed to the improvement of this 
book. 


Leningrad 1962 G. DRUKAREV 


PREFACE TO ENGLISH EDITION 


As the study of atomic physics progresses, the scope of textbooks and mono- 
graphs appears to become more limited, but correspondingly their depth 
increases, Editorial failure to understand and keep pace with this involution 
process would lead to a paucity of books which describe physical behaviour 
in sufficient detail or depth. 

With depth of treatment in mind Professor Drukarev seeks only to 
discuss the quantum theory of electron—atom collisions, omitting detailed 
reference to other atomic and molecular collision processes. Since the 1949 
edition of Mott and Massey’s standard work “‘Theory of Atomic Collisions”’ 
few attempts have been made tocover the more general field, whilst all the time 
the subject has advanced, in principle, in magnitude and in refinement. 

Important contributions to these advances have come from the University 
of Leningrad, where, under the aegis of V. A. Fock, Drukarev and others have 
concentrated their studies on electron collisions, developing the brilliant 
team who have established the neighbouring school at Riga. 

“The Theory of Electron—Atom Collisions” should find a place not only 
with the professional quantum theorists but with those physicists who require 
merely to understand the effects produced by elastic and inelastic collisions, 
without personally manipulating the mathematics. No longer is it possible 
for the more classically inclined ionized gas physicist to aver that the quantum 
theory methods only nibble at the collision problem. Consider the achieve- 
ments as of 1964: the ability to predict elastic scattering and excitation 
functions for simple atoms, with pointers in the direction of success for 
many-electron atoms; the prediction of the existence of virtual and real 
energy levels, and the scattering effects they produce; the establishment of 
upper limits with confidence which is far greater than that usually applied 
to theoretical cross-sections; the prediction of scattering cross-sections for 
electrons with ions and with excited atoms, collisions which have hardly yet 
been attempted in experiment. Only in ionizing collisions is the theoretical 
position rather less satisfactory. 

The successes have come from such treatments as are described by 
Professor Drukarev: perturbation theory, distorted wave calculations, 
variational methods, allowances for exchange and for coupling, together 
with improved computational techniques. 

All this has made it well worthwhile to undertake the translation into 
English of ‘““The Theory of Electron-Atom Collisions”, and a debt of 
gratitude is owing to Mr. Chomet for his capable efforts. 


November 1964 J. B. HASTED 
University College London 
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THE SCATTERING OF A PARTICLE BY A CENTRE OF FORCE 


1. TIME-DEPENDENT THEORY OF SCATTERING 


THERE are two ways of approaching the scattering of a particle by a centre of 
force. Onc of these is to be found in ncarly all the existing textbooks of 
quantum mechanics; it involves the consideration of the continuous-spectrum 
wave function with a known asymptotic behaviour in a given field. This 
corresponds to the scattering of a stationary beam of particles with well- 
defined momentum at a large distance from the scattering centre, and 
represents an idealization of the real state of affairs. 

The other approach is used in the time-dependent theory of scattering 
{1-5} and involves the study of the development of the process with time. 
Here one investigates the collision of a single particle with the centre of force 
rather than the scattering of a stationary beam of particles. The state of the 
particle is represented by a wave packet with a certain momentum distribu- 
tion. This is closer to the conditions obtaining in real scattering experiments 
which nearly always involve the use of various collimating devices. 

Any scattcring experiment can be divided (somewhat schematically) into 
three parts: the formation of the state of the particle, its interaction with the 
centre of forcc, and finally, the determination of its state after the intcraction. 
When the mcasuring and forming deviccs are at a large distance from the 
region of the interaction between the particle and the centre of force, the 
momentum distribution in the wave packet is independent of the distance 
traversed by the particle. We shall simplify the real conditions by assuming 
that the statc of the particle is formed at time ¢ = —oo while the final state 
determination is carried out att = +00. A definite limiting momentum dis- 
tribution will be supposed to exist in the wave packet as t > —oo and 
t— +09. The problem 1s then to establish the relation between the initial 
and final distributions. 

The expressions for the scattering amplitude in terms of the wave functions 
are, of course, the same on either of the two approaches mentioned above. 
However, the time-dependent theory is in some respects more satisfactory. 
It is particularly advantageous in the case of scattering of a particle by a 
complicatcd system, for example an clectron by an atom; and also in the 
analysis of the asymptotic behaviour of the wave function for the system 
along differcnt directions tn configuration space. The time-dependent treat- 
ment is also more convenicnt in the analysis of certain general problems. 

We now proceed to the mathematical formulation of the theory, using the 
system of units in Which /; = m = 1. In this system the dimensions of all the 
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parameters can be reduced to the dimensions of length by multiplying them 
by the appropriate powers of fi and m. For example, the total energy - and 
potential energy V of a particle have the dimensions of [L]~* when multiplied 
by m/f*. With a suitable choice of the scale of length it is possible to reduce all 
the parameters to a dimensionless form. The length characterizing the range 
of action of the field V(r) is a natural unit of length. In the ensuing analysis 
of the scattering of an electron by an atom it will be convenient to take for it 
the Bohr radius ay; the system of units is then known as the atomic system. 
Let M(r,t) be the wave function for the particle. In the atomic system of 
units it will satisfy the equation 
cD 
i 71 H® (1.1) 
where the energy operator H may be written as the sum of the energy of 
free motion Hy = —43A and the energy of interaction V’ between the particle 
and the centre of force: 


We shall assume that V ts spherically symmetric and falls off with distance 
more rapidly than r~' as r>0o. The case V~r~' requires special considera- 
tion and is discussed in § 2.7 

Let @> be the eigenfunctions of the operator H, in the form 


PoE, v,r) = k*(2m)7* el" (1.3) 


The plane wave e!*"” describes the free motion of a particle with energy 
E = 3k? and momentum k = ky. The normalizing factor k4(2r)~# is chosen 
so that 


: oo(E’,v',r)0(E,v, r) dr = 6(E— E’) d(v—v’) 
We expand ® in terms of the eigenfunctions gp: 
O(r, 1) = | eg EW (Ey, 1) O(E,¥, dE dQ (1.4) 


We shall use the symbol dQ to indicate integration with respect to angles over 
all the directions of the unit vcctor ». This expansion represents the wave 
packet mentioned above. The description of a system with the aid of the 
function YW’, is usually referred to as the interaction representation. On this 
particular representation the operator //, is diagonal. The function 
VY (EZ, v, t)is the state vector on the interaction representation. It is physically 
reasonable to assume that when ¢ ~ —oo there exists a definite limit of this 
function, which we shall denote by ‘¥_ (the function ®(7, 7) has no limit when 
t+ —o). When t— +o, the state vector Yo(E. », 1) tends to a definite 
limit which we shall denote by ‘’,. The existence of these limits retlects the 
fact that @(r, ¢) satisfies the equation i@M/cr = #/, for a free particle when 


t Some of the results given below are valid only when the requirements imposed on V 
are more stringent. This will be pointed out at the appropriate places in our analysis. 
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{— —od. The relation between the final state ‘¥, and the initial state P_ 
may be expressed in the form 
W.(E,v) = | S(E, v, E', v’¥ _(E', v') dE’ da’ (1.5) 


The matrix S which is defined by this relation is of fundamental importance 
in the theory of scattering. It follows from 
[\w.[? ae dQ = IW _|? dEdQ (1.6) 
that the S-matrix is unitary, that 1s, 
| S*(E, v, E’,v’) S(E, v, E", v") SE dQ = SE’—E") S(v'—v") (1.7) 


By considering the time-reversed motion, i.e. the situation when ¢ is replaced 
by —¢t and M(t) by ®*(—1F), it 1s possible to show that the S-matrix 1s sym- 
metric with respect to time reversal; that 1s, 

S(E,v, E’,v’) = S(E’, —v’, E, -—v) (1.8) 


where the right-hand side corresponds to the time-reversed motion. 

In order to obtain more detail about the structure of the S-matrix, let us 
consider the Heisenberg representation in addition to the interaction repre- 
sentation, and establish the relation between them. 

Let ~ be the eigenfunction of the operator H which satisfies the cquation 


Ho = Eo 
or in another form 
Ag+(k?—2V)9 =0 


We shall consider those functions @ which have the asymptotic form: 


I 
Pro ~90+0(-) 
and are normalized to 6(E— E’) o(v—y’), as is @y. We shall see later that this 
does not define the function g unambiguously. 
Let us expand ® in terms of the functions g: 


(r,t) = { et WE, y) o(E,v,r) dE dO (1.9) 


where the nght-hand side includes summation over the discrete spectrum. 
The function is the state vector on the Heisenberg representation. 
The relation bctween Y, and is: 


(E,¥,t) = { e EE WE y, E’ vy’) P(E", v') dE’ da’ (1.10) 
Comparison of the expansions (1.4) and (1.9) shows that 
wW(E,v, E’, v’) =| po(E, v,r) @(E’, wv’, r) dr (1.11) 
It will be convenient to transform (1.10) and (1.11) to another,form. Let 
P=PotyZ (1.12) 
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so that in view of (1.2) we have 

CGI—E)y =—V@p (1.13) 
In order to write the solution of (1.13) in an operator form let us introduce 
the resolvent 

G(E) = (H-—E)"! (1.14) 


The operator G is uniquely defined for complex E, since in this case the 
operator H~—F has no zero eigenvalues. The ambiguity arises in the case 
of real E. If we regard G(E) as lim G(E+ie) then the result of the appli- 
e—> +0 

cation of G will depend on the sign of je. It will be shown later that this 
sign, and therefore the resolvent itself (and thus the function ¢ also), are 
unambiguously defined by specifying the initial condition for the state 
vector Y,(E, v, t) as f+ —oo. When expressed in terms of the resolvent G, 
the solution of (1.13) is of the form 


1 = —G(E)V Qo (1.15) 
Substituting (1.15) into (1.12) we have 

~ = Po— GV Qo (1.16) 
We can now employ the equation 

G=G6,-G,.VG (1.17) 


which can be obtained from simple operator algebra. In this expression G, 
is defined, analogous to G, by the formula 


Go =(H,—E)"! (1.18) 
Equation (1.16) can now be written 
~ = ~o—GyVop (1.19) 


which is the integral equation for g. Substituting (1.19) into (1.11) we 
have 


w(E,v, E'. vw’) = 0(E— E’) d(v—v')— } pr(E.y,r,)Go(E V(r) o(E’, v',r) dr 


(1.20) 
Since 
Raxer Vodr = | V(Go mo) dr 
and also (for complex E") 
: | ] 
Gop (Ev, 1) = oe PAE, Vr) = EEF! Pol E, ¥, r) 
we can, using the substitution 
T(E,¥, Ev) =| gS(E,y,) VN) QE nde (1.21) 
rewrite equation (1.20) in the form 
T(E. y, Ey 
HERE HeaShE Gaye (1.22) 


pap 
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We note that this expression contains all the singularities of win an explicit 
form. The matrix 7(E£, v, E’, v’) has no singularities at E= E’. If V(r) 
falls off more rapidly than r~? when r > oo, then T has no singularities 
when v= v’. It will be seen later that the matrix 7 plays an important 
role, since the effective scattering cross-section can be expressed in terms 
of this matrix. We now substitute (1.22) into (1.10), giving 

—i(E’-E)t 


W(E, v, 1) = YE, v) — | T(E, v, E’,v') « W(E’,v') dE’ dQ’ (1.23) 


E-E’ 
In this expression £ is real and E’ complex. When £’ becomes real, there 
is a Singularity in the integrand at the point E = £’. In order to establish 
the behaviour near this singularity, let us substitute the initial condition for 


t = —oo. We shall require that when t = —o the state vectors on the 
Heisenberg and on the interaction representations should be the same: 
Y .(E, v, — 00) = P(E, v) (1.24) 


In order to establish on this basis the contour of integration, let us transform 
(1.23) so that 
Y o(E, v, t) = Y(E, v)— 
: T(E, v, E’, v’) P(E’, v')— T(E, v, E, v’) P(E, Y), ~WE'- aE dO! 
E-E' 
eo E— E)t 
- | TE. V, E,v’) P(E, v) dQ) a (1.25) 
Cc 
The expression 
1 
Spo [ TCE, v, E’, v’) YE’, v’)— T(E, v, E, v’) YE, v’)] 

which enters into the first integral on the right-hand side of (1.25) exhibits 
no singularity at E = E’ and 1s integrable in the interval 0 < E’ < oo. 
Therefore, from the Riemann-Lebesgue theorem we have 


co 


lim | Oe Et dE’ =0 


t-toH 


—i(E’—E)t 
j= | ee | 5, 
C 


The integral 
E'-E 


tends to zero as t + — 00 if the contour of integration in the plane of complex 
E’ is chosen so that the point E = E’ 1s excluded, as shown in Fig. 1. In 
fact, assuming that for t < 0 

—(E’—E)|t| =z 


E.A.C. 


6 THE THEORY OF ELECTRON-ATOM COLLISIONS 


the integral J may be transformed to read 
E|t| 


J= | Sa: 
z 


where the singularity at z = 0 is by-passed, as shown in Fig. 2. When 
It] + co we have 


J= | Saz=0 
2 


This method of excluding the singularity corresponds to the definition of 
the resolvents G and Gy as the limits 


G,(E) = lim G,(E + ie) 


and 

G(E) = lim G(E + ie) (1.26) 

e> +0 
MZ 
ImE\ 
E Re E' Rez 
Fic. 1. Contour of integration in the Fic. 2. Contour of integration in the 
E’-plane z-plane 


We can use the above limits, together with (1.14) and (1.18), to rewrite 
equation (1.16) in the form 


o =-— lim ieG(E + ie)@o (1.27) 
e> +0 


We now consider the expression (1.25) in the limit t-— oo. The first 
integral on the right-hand side vanishes as before in view of the Riemann- 
Lebesgue theorem. However, the integral 


oe E"—E)t 
on | pop “ 
Cc 


does not vanish as f > oo. Substituting (E’— £)t = z when t > 0, we have 


oO 


e7 fF 
Jz= dz 
vA 


— E% 
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When t > o0 and the point z = 0 is excluded by deforming the contour 
of integration into a small circle above this point, we have J = 2zi. It 
follows that 


Ww, = P22 i T(E, v, E', v’) YE, v’) dQ’ (1.28) 
Comparison of (1.28) with (1.5) yields 
S(E, v, E’, v’) = 6(E— E’)[6(v—v’)— 22iT(E, v, E, v’) | (1.29) 


We have thus expressed the singularities of the S-matrix in an explicit form. 

Let us now establish the connection between the effective scattering 
cross-section and the 7-matrix [5]. We suppose that in the initial state the 
particle has a well-defined momentum near Ky which is parallel to the z-axis 
and is localized in the volume U, so that wy = U~*C exp (ikgz), where 
C = | and 


a, ay. a a a, a, 
—~><x<z}, eS SS), -37<7<¢5 
(-Feesg) (-Fers3) (-7<289) 


C = 0 when at least one of the coordinates lies outside this interval. 
In momentum space the wave function will be of the form 


3 . z . . ee 
W_(k) = (=) yond sintk,a, sin}k,a, sin 4(k,—ko)a, (1.30) 
7 k.. k, k, 
If we normalize along the energy scale, then 
dk\? 
Y_(E,v) = k{(—) W_(k 1.31 
Ey) =«(F) ¥-0) (1.31) 
Let 
1 sin4aé 
6,(¢) = a 7 
(1.32) 
2 \? sindaé 
A,(¢) = @ Se 
na € 
so that 
lim 6,(C) = 6(¢) 
pone 1.33 
lim A3(E) = 56) aa 
Since E = 4k? and U = a,a,a, we have 
k \t 
¥ (EY) = (=) 26, (kK) (Ke—ko) 34) 
ay ay 4 


Let j be the particle current density in the initial wave packet. The cross- 
section of an arbitrary packet is a,a,, therefore the total flux of particles 
per unit time 1s 

n_ = a,a,J (1.35) 


2—2 
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Next, the probability of finding a particle moving within the solid angle 
dQ in the final state in the energy interval dE is equal to 


\Y 4(E, v)|? dE dQ 


Therefore, the number of particles which have been scattered per unit time 
into the solid angle dQ within the energy interval dE 1s 


dn, = a,a, j|¥4(E, v)|? dE dQ (1.36) 
In view of (1.15) and (1.34), 
dn, = jl2n { S(E, v, E', v’) 54,(k) 5g, (kj) Ag (ki — ko) 
x(k’ dE’ dQ'|? dEdQ (1.37) 
Substituting (1.29) into this expression, we have, for v # v’, 
dns = (2n)4j| | T(E, ¥, E,¥) 5.,(K1) 50,5) Aag(K— Ko) 
x (E—E')Vk' dE’ dQ’? dEdQ (1.38) 
Consider the integral 
J =| T(E,v, E,v') 5,(K.) 59,(k,) be (ki — Ko) (E- Ek’ dE’ dQ! (1.39) 


It will be convenient to replace integration with respect to E’ and Q’ by 
integration with respect to k,, k, and k;, and replace E’ by 4(ki2 + ke + ki), 
The result of these substitutions is 


ee { T(E, ¥, E’, v') 5,,(K,) 5g, (Kt) Ap, (Kel —Ko) 
1 
x SLE —A(k2 +2 + ki2)] a tkdks dks (1.40) 
Ti 


If it is supposed that a, and a, tend to infinity, we have (provided |Ai| = k’) 


re @) 


J [Te E‘ vA, (k ko) (2 =) oe (1.41) 
— vo ,Y a,\% —K ~~  y : 
Wok) (E-S) F 


0 
and hence 


J = k7?T(E, v, E,v') A, (k— ko) (1.42) 


On substituting (1.42) into (1.38) and taking the limit a, > oo in A,., we 
have, in view of (1.33), 


dn, = j(2n)*|T(E, v, E,v’)|? 0(k —kg)k7 3 dE dQ (1.43) 


When (1.43) is integrated over a small interval of energy in the final state, 
and the result denoted by dN. then 


dQ 
IN, = j(2n)*|T|? re (1.44) 
0 
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Hence, by definition of the effective cross-section do: 
dN, (2n)* 
jo wR 
The total scattering cross-section can be obtained from this expression 
by integrating over all angles; the result 1s 


~~ | IT|? dQ (1.46) 


do = |T(E, v, E, v’)|? dQ (1.45) 


On substituting (1.29) into (1.7), and assuming that v” = vy’ then: 


i{ T(E, v, E, v)— T*(E, v, E,v)] = 22 i IT(E,v,E,v)P dQ’ (1.47) 
or 
~Im T(0) = n | |T/ dQ (1.48) 
where 7(0) = 7(£, v, £,v). The result (1.48) 1s known as the optical 
theorem. Finally, in view of (1.46) and (1.48): 


2*73 


52 
ko 


Im T(0) =o (1.49) 


2. THE WAVE FUNCTION FOR A PARTICLE IN A CENTRAL 
FORCE FIELD 


Let us now consider in detail the eigenfunction of the energy operator H 
defined by equation (1.19) and the relation (1.26). It will be convenient to 
replace @ by another function w defined as 


W(E, v,r) = (22)?k ? O(E, v, r) (2.1) 
In view of (1.3), (1.19) and (1.26) the integral equation for yw 1s 
y = el "_ lim G(E+ ie) Vp (2.2) 
e> +0 


On the coordinate representation the resolvent is an integral operator 
centred upon a Green’s function. The Green’s function corresponding to 
the resolvent Go(E+/e) will be denoted by G,(r, r’, E+ie). In some cases 
we shall omit £-+-/e and write simply G(r, r’‘). This will be done whenever 
there is no danger of ambiguity. In this notation (2.2) may be rewritten in 
the form 


Wer) = eh | Golr, V(r) Wr’) de (2.3) 
Let us determine the Green’s function G(r, r’,) corresponding to the 
resolvent Go(E+ie). In the momentum representation: 
1 eik'rw’ e ik’riv’ pte dk’ dQ’ 
aay | $k'?— (4k? + is) 
If we complete the integration with respect to Q' and substitute 
lr—r']|=p;  k’p=z';— (k? +2ie)p? = 27(e) (2.5) 


G,(r,r’) = lim (2.4) 


e~ +0 
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we obtain 
1 e” dz e' dz' 1. 
Gy = lim —— —_——_— ——§— | = lim —-e” 2.6 
° 0 4T pi | z— 26)” en 2np- — 
and finally 
1 
Giese ik|r-r’| 27 
© 2alr—r'| ° ) 
Therefore the integral equation (2. 7 becomes 
elklr— r’| 
gue: oe Vwedr' (2.8) 
por" 
For large r the asymptotic behaviour of w is 
ikr 
ww ek" + aly, ae where n = : (2.9) 
r 
and the amplitude g is given by 
1 eas 
q(v,n) = — 5 | enter ly, r) de! (2.10) 
I 
Comparison of (2.10) with (1.21) yields 
9) 2 
ee ee T(E, n, E, v) (2.11) 


With the aid of (2.11) and (1.45) we can express the effective cross-section 
in terms of g: 


do = |q\? dO (2.12) 
The optical theorem (1.48) may now be written in the form 
Im g(n, n) = i {lac v)]? dQ (2.13) 


The effective scattering cross-section given by (2.12) is obtained in the 
Stationary theory by considering the asymptotic behaviour of the wave 
function (2.9). This is done by finding the current density associated with the 
scattered particles, i.e. with the term g e™’/r in (2.9). This current density 
turns out to be equal to klg|?.. The quantity do/dQ which is equivalent to 
(2.12) is obtained by dividing this current density by the current density 
of the incident particles, which are described by the plane wave exp (ikv.r). 
The latter current density is, of course, equal to k. The time-dependent 
theory provides the justification for this procedure. The expression (2.13) 
can also be obtained directly by considering the asymptotic form of the 
wave function (2.9). This can be done by setting the current through a 
sphere of arbitrary radius equal to zero, i.e. 


[ Wvy*—yry) ds =0 
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This condition represents the conservation of the number of particles on 
scattering. 
Consider the wave function ¥ at r = 0. On setting r = 0 in (2.8) we have 


w=1- =a ers Vw dr’ (2.14) 
Tt 


The second term in this expression is a measure of the distortion of the 
wave function at r = 0 due to the field V(r). For large k the value of this 
integral is small owing to the rapidly oscillating factors exp (ikr’) and y. 
The maximum distortion occurs at k = 0. In the latter case (2.14) 1s. of 
the form 


WO) = 1-2 | VOW)’ dr’ (2.15) 
OQ 
Applying the mean value theorem to |/(0)—1]| we obtain 
0)—1 vs 
woes = 2 | [Yr dr’ (2.16) 
Wn , 


where w,, is the mean of |W|. Therefore the relative distortion of the wave 
@ 
function is characterized by the integral i Vr dr.t 


0) 
Let us now separate the angular variables in (2.8). To do this, we shall 
expand exp (iky.r) and w into partial waves: 


u ne r) 


er = Y [4n(2 + ty]? Yio(4) Gm) 
l=0 
where @ is the angle between v and r and 
+ 
mtr 
Mu) = ie Ji4g(kr) (2.18) 
k 
y= [42s DY” y,0(0 (2.19) 


By analogy with spectroscopic notation the functions fo(k,r), /,(A,r), 
f.(k,r)... are referred to respectively as the s, p,d... waves. It will also 
be convenient to expand the Green’s function G, defined by (2.7) in the 
form 

a ue 


Go(r,r’) = 2 2 Y,(9, 9) Y¥,,(8'9’) (2.20) 

Y,, are the normalized seiciiaed heals I’, is the one-dimensional 
Green’s function 

rrr) = —ujkr)o(kr’) ‘ r’ 

—okr)uf(kr’) re<r 


{+ This estimate is adequate only if y does not change too rapidly over the range of r 
in which V is appreciable. 


(2.21) 
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where 1, 1s given by (2.18) and v, 1s defined by 
1 /nkr\* 
pe (5) H{D, (kr) (2.22) 
Substituting these expansions into (2.8) we obtain the integral equation 


for f(r): 


r a 


f(r) = u,(r) + 20,(r) { u(r Vr) fir’) dr! +2u, | vr) Vr fir) dr’ — (2.23) 
0 r 
The asymptotic form of f; when r > 0 1s 
j | ree 
fi(r)~ join (11 — 3) + qe ( z) (2.24) 


where q, is the partial scattering amplitude given by 
gq =—2f (Var (2.25) 
0 


The amplitude g can be expressed as a sum of the partial amplitudes 
if we substitute the above expansions into (2.10): 


q =), [4n(2l +} *ar¥,, (4) (2.26) 

The total cross-section o can also be expressed as a sum of the partial 
cross-sections ¢;: 

Co = \ 0; (2.27) 

o, = 4n(21+1)Iq,]? (2.28) 


The optical theorem can be used to establish the relation between the 
imaginary part and the modulus of q,. Substituting (2.26) into (2.13) we 
have 


sow 
Since the amplitude g, may be written in the form 
qi = |qi (2.30) 
there follows 
stl, 2 1 
q, = eft sin 5) = re $11) (2.31) 
or 
= (2.32 
4 cot 5,—ik 22) 


We can now use (2.31) to transform the asymptotic expression given by 
(2.24) so that it includes the phase 6,. After some rearrangement we find: 


eft In 
fix ras kr — SO, (2.33) 
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It is clear that the relationships given by (2.30) and (2.33) define the phase 
only to within a multiple of zx. This limitation can be removed by means of 
the natural condition 

lim 6, = 0 (2.34) 


E-@ 


The partial cross-section a, is given in terms of 6, by 


og, = 4n % sin* 6, (2.35) 
It is evident from this formula that 
G,<(2i+ Ws (2.36) 
The function /, 1s the solution of the differential equation 
wo + (2-52 av) j.=0 (2.37) 


where 
f(0) =0 (2.38) 


The function f(r) has the required asymptotic behaviour at large r. 
In addition, let us consider the functions which have the given behaviour 
for small r. We shall denote them by x,(r). In some respects they have ad- 
vantages over f(r). Since both functions are solutions of the same differential 
equation (2.37) and both vanish when r = 0, their ratio must be independ- 
ent of rand 1s a function of k only. We may therefore make the substitution: 


x(k, r) = L(k) fk, 1) (2.39) 


and require that the functions x, should behave as u,kr) for small r. 
Substituting (2.39) into (2.33) we have, after some rearrangement 


r 


Xr) = u(r) + [o(r)u(r’) —u(r)o,r’)|2V(7')x,(r’) dr’ (2.40) 


0 


L(k) = 1-2 | »,Vx,dr (2.41) 
j 


We note that the function v, has a pole of the form r~' at r = 0. Moreover 
V has a Coulomb singularity at r = 0. It follows that v,V has a pole of the 
form r~'~! or r~'. However, since f; falls off as r’** when r > 0, it follows 
that the integrands in (2.40) and (2.41) remain finite at the lower limit. 
The function L, defined by (2.41) is identical with the so-called Jost 
function [6]. In view of the connection between the Bessel and Hankel 


functions, we have 
vp = k((-1)'* ' uy — in] (2.42) 
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We can use this relation to rewrite (2.40) and (2.41) in the form 


x, =u,+(- 1k | Luar’) are RV OIL) dr’ (2.43) 


L=14(-1)'k [ u_)-12Vx; dr-+ik | 4,2V%x, dr (2.44) 
0 0 


It follows from the properties of Bessel functions that at given r the expres- 
sions uk‘ and w_,_,k'*! are analytic functions of k*. Therefore k7 ‘x is 
an analytic function of k* at fixed r. ‘Hence, L, may be written in the form 


L, = Ak) + ik?!* Bik) (2.45) 


where A, and 8, are even functions of a real k, defined by 
A, = | +(— 1k | HW ecpe 1 2Vx,dr 
O 


. (2.46) 
B, = pet [ 12Vx,dr 
0 


The partial scattering amplitude g, can be expressed in terms of L,;. Using 
(2.39) and (2.46) we can rewrite (2.25) in the form 


k?"B,(k) 
= = ee ay 2.47 
q1 L, (2.47) 
Let 
L, = |L,| e* (2.48) 


The quantity g, is found to be closely related to the phase 6, introduced 
above. This can be shown by substituting (2.48) into (2.47) and utilizing 
(2.45): 


ae 
qu=- ,e sin Q; (2.49) 


Comparison of (2.49) with (2.31) then shows that go, =—6,+22N. The 
term 2x/N is unimportant and can be removed by a Suitable choice of origin 
of @,. We shall therefore set 


The properties of L, when k 1s complex are of interest. They can be 
conveniently investigated with the aid of (2.41). By using the asymptotic 
expressions for u, and v, and the expression for x, given by (2.40) we can 
verify that L, is analytic at least over the upper half-plane. Its properties 
in the lower half-plane depend on the behaviour of V(r) as ro. If V 
falls off as a power of r, that is V ~ r~", then in the lower & half-plane L 
is not analytic. If, on the other hand, V ~ e°*/r, as is often the case, 
then L, is analytic in the lower half-plane within the band |Im 4] < 3a. 
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If the potential energy V is negative, then the existence of bound states 


with energies —|E,|, —|E2|... is, in general, possible. Negative energies 
correspond to imaginary wave numbers k; = ik; = iV 2)E)). It may be shown 
that the function L,(k) vanishes at the points ik,, ik2,.... Let us substitute 


k = ik into (2.40) and require that x,(ix, r) vanishes when r > oo, as in a 
bound state. In view of the asymptotic behaviour of the Bessel functions 
of complex arguments, we find that as r > o the coefficient of u,(r) must 
vanish, which in turn implies that L, = 0. It may be shown that all these 
roots are simple ones. 

When V(r) behaves at infinity as a Coulomb potential, as for example 
in the case of the scattering of an electron by an ion, the above equations 
must be modified in view of the particular asymptotic properties of the wave 
functions for the continuous spectrum of the Coulomb field. To make 
this modification, let 


V(r) = W(r) - ~ (2.51) 


considering the attractive interaction in which we shall be interested later. 
We introduce the Green’s function Gc for the Coulomb field: 


Z 
(44 +4k* + =) Go = —d(r-r’) (2.52) 
Instead of (2.3), we now have 
W = Ve —| GoW War’ (2.53) 


Where the function Wc is the wave function for the continuous spectrum 
of the Coulomb field whose asymptotic behaviour for r > © is 


We oe gilkverty In Cky key) Ie 1 eikrty In 2kr (2.54) 
r 
where 
Z 


and qc is the scattering amplitude in the Coulomb field and is given by 
. 98 
i y (i+ y) 1 givin 2 sin? 5 (2.56) 

2k T(1—iy) . 48 

sin 7 


The problem is to find the solution of (2.53) which has the asymptotic form 
given by (2.54) but with an amplitude corresponding to the potential energy 
(2.51). In order to separate the angular variables, let 


ur(r) 


r 


Wo=) [4n(2l+ D]Fi Yio(9) (2.57) 
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where wr is the regular solution of 
d7 ul li+1)  2Z 
—~ + ais a= uf = 0 (2.58) 
dr r r 


which, when r > oo, has the asymptotic form 


ur ~ “sin (ir ae + yln 2kr + nee in 2kr ) (2.59) 
where gc , 1s the partial scattering amplitude in the Coulomb field: 
Goi = rating (2.60) 
We shall take the Green’s function Ge in the form 
Ger.) =T 2-2 ¥,0,0) ¥,(0'0) (2.61) 
where 
ranrefTBOED 18% am 


where ur is defined by (2.58) and (2.59), and vf is a non-regular solution of 
(2.58) whose asymptotic form is 


I 
(neo 44 In 2kr 


vp~ —e (2.63) 
We write the function w in the form of a series analogous to (2.57): 
ra ti) 
I 


For the functions f/f; we obtain 
[o 8) 


f(r) = uf(r) +2 D, r(ryul(r Wor) fie’) dro + 2 fu r(ryekr Wr f(r’) dr’(2.65) 


r 


For large r we have 


In 
r- ae In 2hkr 


| Ir 
fir) ~ > ;, sin (kr = 3 ——+ yn 2kr + (qceitaye ic (2.66) 


where 


oO. 


pe | u(r WO f(r) dr (2.67) 
0 


It follows that the partial scattering amplitude gq, in the field (2.51) ts 
given by 
N= Ici (2.68) 
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If we now let 


q, = = (er I), = dear = = 5 (et 1) (2.69) 
and 
Op = Nt ey (2.70) 
we obtain 
§, = = e2!M (92! _ 1) (2.71) 
2ik 


When r > © we have, using (2.66) and (2.69), 
|. 34 ! 
fi~ 7 e sin (i — % + yIn2kr + i) (2.72) 
The integral equation (2.65) will now be transformed to a form analogous 


to equations (2.40) and (2.41). To accomplish this, equation (2.65) will be 
rewritten in the form 


fir) = ur) [ v2f rr ir Wr’) in| ae 


+2 f [of (r)uf(r')— v6 ru Wr) dr’ (2.73) 


If we substitute 
yi(kor) = Li(Wf(k, 


we obtain 
vir) = uF) +2 [ [oF uF(r)— oF WENO ylr'dr’ (2.74) 
0 
LS(k) = 1-2 f ofiry, dr’ (2.75) 
Consequently, ; 
—2 fucwy, dr 
G,=— (2.76) 


2 { v Wy, dr 
0 


3. SCATTERING AMPLITUDE AND PHASE 
It is evident from equations (2.47) and (2.50) that the properties of the 
amplitude g, and phase 6, are determined by the quantity L,. Using equations 
(2.48) and (2.45) we can express cot 6, in terms of 4, and B;: 


fits 


cotd, =— mies (3.1) 
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Hence it follows that cot 6, is an odd function of k. Since 6(+00) = 0 we 
find that the phase 4, 1s itself an odd function of k, 1.e. 


6(k) = —6(—k) (3.2) 


On this definition 6, will in general experience a discontinuity at k = 0. 
Since L, vanishes on the imaginary axis at points corresponding to bound 
states, we can establish with the aid of (3.2) an important relationship 
between the number JN, of levels with given / and the quantity 6,(+0), 
i.e. the limit of 6,(k) for k —> 0 on the side of positive k. This relationship 1s 


6 +0) = N,2 (3.3) 
and is known in the literature as Levinson’s theorem. In order to prove 
this result we recall that, according to Cauchy's theorem, the change in 
phase of L,(k) for integration over the contour shown in Fig. 3, 1.e. the 
quantity —[6(—0)—6(+0)]/2z, is equal to the number of zeros of L,(k) 
within the domain enclosed by this contour in the plane of complex k. 
All the zeros of L, lie along the imaginary axis and their number if AN). 
Consequently, 

—[6(—0)-—6(+0)] = 2zN, (3.4) 
If we now combine equations (3.2) and (3.4) we obtain the result given 
by (3.3). 
Imk 


Rek 


Fic. 3. Levinson’s theorem. The points on the imaginary axis are the zeros of the 
function L 


A special case arises when L,(0) = 0. This situation occurs, for example, 
when one of the discrete energy levels merges with the boundary of the 
continuous spectrum, i.e. tends to zero as the depth of the potential well is 
reduced. Under these conditions L,(0) = 0 and the phase is given by 


§(+0) = Nx +> (3.5) 


where AN is the number of levels (not including the special level with zero 


energy). 
Let us consider the limiting value of g, as k ~ 0. We shall start with the 


expression given by equation (2.47). Substituting equation (2.45) into 


(2.47) we obtain 
Ak) \~? 
um (ik oo | (3.6) 
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Suppose that A,(0) 4 0 (when A,(O) = 0, it follows from equation (2.49) 
that L,(0) = 0,1.e. we have the special case mentioned above). It follows that 


lim q, =0 when / #0 
B,(0) (3.7) 


i.e. all the partial amplitudes except the zero-order amplitude vanish and 
the amplitude of the s-wave is, in general, the only one to have a non-zero 
value. The limiting value of this amplitude is usually denoted by —a, where 
ais the scattering length. Thus, 


_ Bo(0) ais 
A (0) 
According to equation (3.1) we have 
] 
lim k cot dg) =—- (3.9) 
k>0 a 


We can use the latter relation to show that the scattering length is equal, 
apart from the sign, to the derivative of the phase angle 6 with respect to 


the wave number X, 1.e. 
(28 310 
a=—{— 


If there is an s-level with zero energy, then the amplitude g(k) increases 
without limit as 1/k when k 0 (resonance scattering). Under these 
conditions, equation (3.10) is no longer valid. If there is a p-level with 
zero energy, then 

lim A) 

k-~0 k* By(k) 


is finite and therefore 
lim g, #0 
k-+0 
The partial amplitude for /= 1 has a non-zero value when & = 0. 
Partial amplitudes for / > 1 vanish independently of whether a zero level 
with given / exists or not. 
Let us consider in greater detail the sign and magnitude of the scattering 
length. From the definition given by equation (3.8) and the expressions in 
(2.46) we have 


2 | Vx,rdr 
q = —__—___ (3.11) 
142 | Vx dr 
0) 
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In a repulsive field V > 0. For zero energy, the function xg satisfies the 
equation 


2 r+2| (r—s)Vx¢ ds 
0 


and is therefore everywhere positive. It follows that the scattering length a 
is positive. In an attractive field V < 0, and if there are bound states the 
function XxX) has no nodes (apart from the origin) and is everywhere positive. 
Therefore the integral in the numerator of (3.11) is negative. The denomina- 
tor remains positive as V passes from positive values through zero to negative 
values. The denominator in equation (3.11), which is equal to A,(0), 
vanishes (and the scattering length becomes infinite) only when a bound 
State with zero energy appears. In an attractive field the scattering length 
is negative in the absence of a bound state. In the presence of bound states, 
the scattering length can be either positive or negative, since the sign and 
magnitude of the integrals in equation (3.11) depend on the number and 
position of the nodes of the function xg. All the possible cases can be 
conveniently illustrated by introducing a parameter a, which characterizes 
the intensity of the field of force, and considering the dependence of a on a. 
Let 
V(r) = aV,(r) (3.12) 
and assume that VY) > 0. On this definition a positive a corresponds to 
repulsion and a negative @ to attraction. 
On the basis of the general properties of functions satisfying Volterra’s 
integral equations, we can state that for given A and r the function Xg will 
be an analytic function of « It follows that @ will be a meromorphic 


function of «. 
Figure 4 shows the general character of the dependence of a on a. The 


0 
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Fic. 4. Dependence of the scattering length and the energy levels Eo, E,,... (broken 
curves) on the parameter a 
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points &, a,,... are the values at which the first, second, ..., levels with 
zero energy appear as the depth of the potential well increases. The broken 
curves represent the dependence of the energy of the first, second, and so on, 
levels on the depth of the potential] well. As can be seen from Fig. 4, the 
scattering length, and with it the scattering cross-section at zero energy, 
become equal to zero (Ramsauer effect) for certain values of the parameter «. 
As an example, consider the dependence of the scattering length of an 
electron in a Thomas-Fermi field due to a nuclear charge Z [7] (see Fig. 5). 
Figure 5 is a mirror image of Fig. 4, since —Z plays the role of «. 


~100 1 ° 
—200-F > I 
—300 


Fic. 5. Dependence of the scattering length for electrons (in units of 0-885/Z) on the 
nuclear charge Z in the Thomas-Fermi field 


Let us now consider the behaviour of the amplitude g, as a function of 

k for small k. If the potential energy falls off as r > co more rapidly than 

any negative power of r, then ZL, will be an analytic function not only in the 

upper half-plane of complex k but also in a band in the lower half-plane. 

Consequently, A, and B, can be expanded into a series in powers of k 

(within the limits of convergence), and in view of the fact that A, and B, 

are even, these will in fact be power series in k? when k is real. We can 
then write 

Atk) AlO+k7AfO+... 

Bk) BJO + kK? BEY +... 


If k is so small that B, vanishes in the interval between 0 and k, we can 
expand 1/(B(?+k7Bi + ...) in powers of k’, giving 


(3.13) 


Atk) Am 1 Fo _ APB” -| (3.14) 
Bk) BY LB BYP” | 
Let 
Aj? 1 Aj? ALO BP 
BO oe BO ~ BO2 =—4Hyr, (3.15) 
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so that equation (3.14) may be rewritten in the form 


A, 1 
Sr ea (3.16) 
B, a, 
On substituting the latter series into equation (3.1) we have 
] 
k?**! cot 6, =— ~ + 4k2r,4+... (3.17) 
ay 
We can now use equation (2.32) to obtain the following expression for gq,: 
Oke 
q, = (3.18) 


l 
—-~-+ dk?r, +Fik?!*! 
ay 


In the special case / = 0, the expansion given by (3.17) becomes identical 
with the well-known formula [8] 


1 
kcotdy9 =— Pa +4k?r, (3.19) 
0 


originally obtained by a different method. The quantity r, in this expression 
is called the effective range of the forces. 

If for given / there exists an energy level such that x, < 1, it is possible 
to obtain a relation between a,, r; and x, by using equation (3.16) and 
assuming that A = ik. Since L,(ik) = 0, it follows from equation (2.45) 
that 


Asin) _ 


Bix) (—1) 732! (3.20) 
Therefore we have 
(—1)ire77*} = -+ LK?r, (3.21) 
and in particular when / = 0 
- = K(1+4xkro) (3.22) 


0 


If the potential energy falls off as r~" when r > oo, then L, is not analytic 
in the lower half-plane of complex k, and therefore A, and B, cannot be 
expanded in powers of &. It follows that the expressions given by equations 
(3.19) and (3.22) will no longer hold. 

As an example, consider the energy of interaction between a point 
charge and an induced dipole which at large distances ts of the form 

V~— p (3.23) 

r 
This expression 1s encountered, for example, in the case of scattering of an 
electron by an atom when polarization of the atom by the incident electron 
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is taken into account. In this case, (3.19) is replaced by [9] 
1 7p 4p Bk 
kcotdg =—-+—5|kl + —k* In(—+) +... 3.24 
oa a aa 5, n(FP) + — 


When / # 0 we have instead of equation (3.17) 


_ 8043+ )I-4) | 
= . 


Therefore as k > 0 the amplitude g, falls off as A for all /. 

Let us also consider the behaviour of the amplitude as k — 0 when the 
field is of the Coulomb type. According to equation (2.68) the analysis 
can be reduced to the investigation of g,, since the scattering amplitude 
dic 1n the Coulomb field is known. We shall confine our attention to the 
case / = 0. 

The dependence of G on k can be analysed with the aid of equation 
(2.76). However, it will be useful to express the functions vy and uy (with 
given asymptotic behaviour) in terms of a function having a given behaviour 
at the origin. Let x, and x; be two linearly independent solutions of the 


k* cot 6, (3.25) 


equation 
d*x. 2Z 
ila («2 rn =) 20 (3.26) 
dr r 
where 
Xx,7r 
oe (3.27) 
when r > 0. We then have 
u& = eC x,(r) (3.28) 
fi 
py a e '7 (= s+ ikC9x,) (3.29) 
Co 
where 
2nZ/k \? 
ene ae) (3.30) 


When k < 2nZ, the quantity Cy is approximately equal to (2xZ/k)*. It 
follows from equations (3.28) and (3.27) that as A > oo the function u& 
increases as k~*. According to equation (2.74), the function u,(r) should 
also increase as kK~* when r—> 0. It is evident from equation (3.25) that 
v© falls off as k? as k > 0, so that the product v“y, remains finite as r > 0. 
The quantity q, increases as 1/k when k — 0 (this is evident from equation 
(2.76)). 

So far, we have discussed the general properties of the phase 6 and 
amplitude g without considering their numerical values. Let us now proceed 
to the calculation of the phase 6. In principle, this requires the solution of 
equation (2.37). However, this equation has an analytical solution only 
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for a certain special form of V. In all other cases it can be solved only 
numerically. Equation (2.37) is inconvenient for the evaluation of the 
phase, since the function /; is oscillatory and does not tend to any definite 
limit as r > ©O. 

However, it is possible to reduce the problem to the integration of a 
first-order non-linear equation for the function y,(r), whose limit as r > 00 
is equal to the phase 0, [10], and which exhibits monotonic variation for 
monotonic V. 

Consider the logarithmic derivative 


1 df, 
=-— 3.31 
~ ai (3.31) 
It satisfies 
d I+ 
gay (3.32) 
dr r 


and may be written in the form 


d du_)_ 

ert ae ni{r)+(- 1)! el! sin yr) 

_ dr dr (3.33) 
an u,cosn(r)+(— 1)'u_,-4 sinn,(r) | 


If we substitute equation (3.33) into (3.32) and bear in mind the equations 
satisfied by u,and u_,_,, and the fact that 


du _,- d 
ul, —— ae psy 7 = (-1)'*'k (3.34) 
we obtain 
dn 2V 
7” = —, Lin cos n(r)+(—1)' u_)_, sinn(r)]} (3.35) 


For large r the function n,(r) in equation (3.33) may be replaced by 1,,(oo), 
so that by using the asymptotic form of the Bessel functions we have 


d ! 
o~ 0 sin (ir — 3) cos n,(c0)+ 


l d ! 
+ cos (ir — >) sin n(co)| = — Insin kr 22 +n {0o)| (3.36) 
2 dr 2 


On substituting equations (3.36) and (3.31) we find that the asymptotic 
expression for f; which corresponds to equation (3.36) is 


Ir 
fi(r) ~ constsin kr aac n(o)] (3.37) 
from which it follows that 


(co) = 0, (3.38) 
Consider the behaviour of y,(r) for r > 0. 
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Since f, ~ r'*! for small r, we have 
lim gr =1+1 (3.39) 


r-Q 


On the other hand, it follows from equation (3.33) that for small r 
1—1(1+1)y,(r) (kr) 777! 


= (+1 4 
Or OT en en oe 
Hence in view of equation (3.39) we obtain 
lim 4,(r)(kr)7 777+ =0 (3.41) 
r-0 


i.e. numerical integration of equation (3.35) should be carried out subject 
to the condition n,(0) = 0 for all /. 


4. DISPERSION RELATIONS 


Consider the matrix element 7(E, v, £, v’). According to equations (1.21) 
and (1.16) 


T(E, v, E,v') = | o*Vo,dr — lim otVG(E+ieVo dr (4.1) 
e~ +0 


If we recall the explicit form of the function @ 9 given by equation (1.3) and 
replace the resolvent G by the Green’s function G(r, r’, E+i0), we have 
in the case of forward scattering (v = v’) 


k 1 ; 
T(E, v, E,v) = ral V dr — | kyr) V0r)G(r, r’)V(r') dr dr’ 


7 (22)? 4x 

(4.2) 
If we multiply the latter expression by —(2z)*/k and let — f Vdr = 27qp, 
we have in view of equation (2.11) 


| 
q-GQo= e72 eerient) Vir)G(r,r')V(r') dr dr’ (4.3) 
4 
Let 
G-g =F (4.4) 


and consider the analytical properties of F in the plane of complex energy 
z = E+ fe cut along the positive real axis; this corresponds to the upper 
half-plane of complex k. The properties of F are determined by the prop- 
erties of the Green’s function G(r, r’, z). For a sufficiently rapid decrease 
in V(r) at large distances (faster than r~ °*), the function G has the following 
properties in the z-plane [11]}: 

1. G(r, r’ z) is analytic in the z-plane with a cut along the positive real 
axis and has, in general, simple poles on the negative real axis E = —|E;j. 
These poles correspond to discrete energy eigenvalues in the field V. Using 
the eigenfunction expansion of G 

eo 
G(r,r’,z)=)> PAN (4.5) 


J Jj 
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or in the presence of degeneracy, 


* , 
G(r, r’,2) = 5 CO) (4.6) 
ja SE; 
we see that the residues at the poles z = E; are equal to 
2. Y ACW) 
2. The limit lim G(r, r’, E+ie) exists for all E > 0. 
e—>+0 

3. G(r,r’, E+i0) = G*(r, r’, E—i0) (4.7) 


4. For large |z| 
IG—G,| < Cue ene 


where Gz is the Green’s function for free motion (2.7); Ce ~ O when E > oo. 


Fic. 6. Derivation of the dispersion formula. £), Eo, ... represent the bound-state 
energy levels 


It follows from these properties that the integral on the right-hand side 
of equation (4.2) converges absolutely for all z and is an analytic function 
of z everywhere on the cut plane except for the poles of the Green’s function 
G; moreover, 


F(E+i0) = F*(E— i0) (4.8) 
lim F(z) =0 (4.9) 
|z[ +a 


We now derive the dispersion relation [11, 12]. Consider the contour C 
in the z-plane (Fig. 6). From Cauchy's theorem 


Lf F(2’)dz’ d, 
F(z) = — | 2 - 4.10 
@) mal vo} EL + oe) 
Cc 
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where 


1 | 
dj == ) { ery Vy, dr | er Vy dr’ (4.11) 


are the residues of F at the poles z’ = —|E;| (we recall that d; are real). As 
the radius of the circle tends to infinity, the contribution due to the integrand 
on the circle tends to zero in view of equation (4.9). There remains only the 
integral along the boundary of the cut: 


F(E’ + 10 I F(E'—i0 d, 
Fe) = 5 ae - [ae 4 (4.12) 
2ni E=Z 2m E'-—z j |Ej|+z 
0 O 
From equation (4.8) we have 
Im F(E’ ) AE 
F(z) =—- 4.13 
) ae E’-—z BD oy 
O 
Now let Im z > 0 so that z > E. Since 
lim ————- = FP ——. + ind(E'-E 4.14 
eto ee oe! 
where ¥ represents the principal part, we have 
Im a d; 
E’ — 4.15 
Re F(E) = to [ er LEI+E (4.15) 
Since 21g) =— f V dr is real, the expression given by equation (4.15) may 
be written in the form 
1 Im q(£’) dE’ 
Reg = qo dy [ maeae | YY (4.16) 
T E j zi] 


0 


We can now use the optical theorem (2.13) to express Img in terms of the 
scattering cross-section: 


I a 
mq 7 


so that 


reas = ao gia | SESE —y win 


ae 


It is interesting to consider the limiting case of equation (4.17) when £ = 0. 
Since the amplitude g at zero energy is related to the scattering length by 
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q(0) = —a, it follows from equation (4.17) that when E = 0: 


I d 
am ao~ za] odk+ TH (4.18) 


The importance of the dispersion relations in collision theory lies in that 
they can be used to derive rigorous formulae for quantities directly related 
to measured values, and can therefore serve as a criterion for the analysis 
of experimental data (see § 16). 


Chapter 2 


APPROXIMATE METHODS FOR THE DETERMINATION OF 


WAVE FUNCTIONS AND SCATTERING AMPLITUDES 


5. PERTURBATION THEORY 


IN § 2 we estimated the effect of the field of force on the wave function near 
the centre of force for k = 0, and saw that it could be characterized by the 


quantity { Vr dr. 


co 0 


If { |V|rdr < 1, this effect is small for k = 0 and even more so when 


0 
k # O, so that V can be regarded as a small perturbation. Let V = «Vp as 


in equation (3.12) and expand the solution of (2.8) into the series 
y= Dey 

Let us take a plane wave for yO 
yy — gikvr 

The expression for y™ is then 


gikir— r ies 


The scattering amplitude qg will also be expanded into a power series: 


q = s q@(or*} 
It follows from equation (2.10) that 


i 
gq” oe = |ew Vow” dr’ 
20 
The first term in the expansion (5.4), namely 


gq Sai | Vo po ike yr dr 
2i 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


represents the well-known Born approximation. If Vo is independent of 
the angles, we can carry out integration with respect to the angles in 


equation (5.6). The result is 


Pee 
g0 = — 2| Vo "2 dr 
r 


0 
29 


(5.7) 
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The quantity 
80 
K =2k sin5 (5.8) 


is the momentum (in atomic units) transferred to the centre of force as a 
result of scattering through an angle 0, where @ is the angle between v and v’. 
From the formula 
sin Kr ur(kr 
= 5 [4n(2l+ }*¥,, 0) 2 
Kr F 
and using equation (2.26), we obtain the following expression for the partial 
amplitude on the Born approximation: 


(5.9) 


r2 


qi =— 2 Vou? dr (5.10) 
0 


It is evident from the expression for qg, given by equation (2.47), and from 
equations (2.45) and (2.46), that the partial amplitude q, is a meromorphic 
function of the interaction constant «. Consequently, the expansion given 
by (5.4) (and the equivalent expansion for the partial amplitudes) is only 
valid for values of a smaller than the radius of convergence (a detailed 
investigation of the convergence of the series for g in powers of a@ is given 
in [13]). 

However, the expression for the partial amplitude can be written in the 
form of a ratio of two series of powers of a [14], which is valid for all a. 
With this purpose in mind, let 


v= 2 ocx (") (5.11) 


in (2.40) and assume that 
x) = Uy (5.12) 
Then, 


xi) = 2k(-1)'*! (ua (Nu lr’) —u_ yr u(r) ]xt"" D Vg dr’ (5.13) 
; 


It follows from the general properties of Volterra’s integral equations that 
the series given by (5.11) converges for all a. 

If we substitute equation (5.11) into (2.44), we obtain the following 
series for the Jost function 


b=) a? (5.14) 
where 
oe | (5.15) 
LD) = 2(- 1)'k | Vouju_,-,dr+2ik Vou dr (5.16) 
0 0 


1) = 2(—1)'k | wpxf uy, drt Qik | Voupxf? dr (5.17) 


Q 
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Finally, substituting equation (5.14) into (2.47) we obtain 
1) «Im Ly? 
a: yo” 
The approximate expression for g, can be obtained by cutting off the series 


in the numerator and denominator. In particular, if we retain first-order 
terms only, we have, in view of equation (5.15), 


(5.18) 


a / Im LP 
sie ee ate (5.19 
_ k ( Sa) 
or on substituting the explicit expression given by (5.16) 
2] Vuj dr 
a ee (5.20) 


1+2k | Vuuju yy dr-+2ik [ Vu? dr 
0 0 


Comparison of (5.20) and (5.10) shows that the Born approximation can 
be obtained from (5.20) if the integrals in the denominator are neglected 
in comparison with unity. 

As an example of the application of Born’s approximation consider 
the scattering of a charged particle by a distributed charge of density p 
where 


i pdr=Z (5.21) 
In this case, the amplitude g® can be directly expressed in terms of p. 


Suppose that Z, 1s the charge on the particle. Its potential energy V must 
then satisfy the equation 


AV =—4nZyp (5.22) 
Moreover, we can write 
A elk-r noe — kK2 eik-r (5.23) 


On multiplying equation (5.22) by exp (iK.r) and (5.23) by V, we have on 
subtracting and integrating over all space 


i [elr AV —VA(e!K")] dr = K? | eK V dr—4nZ, |e pdr (5.24) 


According to Green’s theorem, the integral on the left-hand side vanishes 
if V has no singularities and falls off sufficiently rapidly as r > oo. If this 
iS SO, 


4 
[ems Vdr=— =| er 5 dr (5.25) 
so that 


22 ; 
g = =| pdr (5.26) 
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The quantity 


1 
F(K) = =| er odr (5.27) 
is known as the form factor. It can be used to rewrite equation (5.26) as: 
2202 
gi = ae F(K) (5.28) 


If Vir) has a Coulomb singularity, this singularity must be considered 
separately. For example, when the scattering system is an atom of charge Z, 
we have 


Z,Zo 


V= +W (5.29) 
where W is finite at r = 0. We then have 
2Z 
gi) = ra [(Z,+ZF(K)] (5.30) 
For a neutral atom, 
Z = —Z,; ($.31) 
so that the formula given by equation (5.30) becomes 
=e a 
g = =8=TF(K)-1] (5.32) 


6. THE METHOD OF TRIAL POTENTIALS 


When the potential energy V cannot be regardcd as a small perturbation, 
one can try to choose an expression for V, such that there is a rcasonably 
simple solution for scattering in this ficld, and at the same time V—V, can 
be regarded as a small perturbation [15]. If we confine our attention to 
simple trial expressions for Vy) which contain only one or two parameters, 
we cannot, of course, expect to obtain a good approximation to the exact 
wave function and the exact amplitude over a wide range of encrgies and 
angles 0. 

However, the situation is more favourable when AR < 1, where R is 
the range of the ficld V and w is a slowly varying function of & and 0. Having 
chosen V, so as to achieve the best approximation to w at zcro energy, we 
can use the same expression for V, to obtain an adequate approximation 
in the entire intervalO < AR Sl. 

It is natural to choose V so that there are no effects of the first order 
in V—V5. 

Consider the limiting case of zero energy. The wave functions wy and 
Wo corresponding to V and Vg respectively satisfy the equations 


Ay —2Vy =0 (6.1) 
Ao —2Vo Wo =0 (6.2) 


2. METHODS FOR THE DETERMINATION OF WAVE FUNCTIONS 33 
and have an asymptotic behaviour of the form 


q 


w=1+ . (6.3) 
Uo =l+ = (6.4) 


Let us multiply equation (6.1) by Wo, (6.2) by w and then subtract and integrate 
over all space. In view of (6.3) and (6.4) we have 


2n(q—4o) =| (V—VoWbow dr (6.5) 

We shall require that g = qo, so that from equation (6.5): 
| (V—V,)Wowdr =0 (6.6) 

On a first approximation we can replace & by Wo, so that 
{ (V —V>)W2 dr =0 (6.7) 


which is the condition for first-order effects to vanish in perturbation theory 
(see, for example, [17] Section 10). 

Consider now the condition which must be satisfied by the difference 
OV = V—V, for it to be looked upon as a small perturbation. Let 


by = Y—Yo (6.8) 
The equation for dw follows from (6.1) and (6.2) and Is 
Ad —2V, dy = 26Vy (6.9) 
Furthermore, consider the irregular solution of the equation 
Ag—2V,9 =0 (6.10) 


which behaves as 1/r when r— 0. If we multiply equation (6.9) by g, 
(6.10) by dy and then subtract and integrate over all space, we obtain, taking 
into account the singularity in g at r > 0, 


i 
5Y(0) = = { 5Vow dr (6.11) 


The order of magnitude will be retained if @ is replaced by 1/r in the integrand. 
Hence, 


|5Y(0)| S(max |p])2 | [SV |r dr (6.12) 

0 
The quantity |dy(0)|/max |W| may be looked upon as a measure of the 
relative departure of w from Wo» within the sphere of action of the field of 


force. It is evident from (6.12) that the condition for this departure to be 
small 1s 


2 | \sV\rdr<1 (6.13) 
0 
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When V, = 0 this condition becomes identical with the condition for the 
applicability of perturbation theory for V (§ 5), as should be the case. 

If the trial expression for Vg contains one parameter, equation (6.7) can 
be used to define it. If equation (6.13) is then not satisfied, it is necessary 
to introduce additional parameters in order to reduce the magnitude of 


{ \OV|r dr. An example illustrating the possibilities of this method will be 
0 


considered in § 13. 


7. THE VARIATIONAL METHOD 


The problem involving the scattering of a particle by a centre of force 
may be formulated as a variational problem involving the determination 
of the stationary value or an extremum of a functional. This formulation 
is used to establish certain general relationships such as, for example, the 
virial theorem [16], and also in the development of the variational method 
for approximate calculations. 

In the variational method one chooses a suitable analytical expression 
for the continuous-spectrum wave function, 1.e. the so-called trial function 
which contains a number of parameters. The parameters are determined 
by requiring that the functional which includes the wave function should 
exhibit either an extremum or a Stationary value. 

Apart from the latter condition, there are a number of integral relations 
which the wave function must satisfy. The exact function satisfies all the 
conditions whereas the trial function may satisfy only a few of them. The 
various modifications of the variational method differ from each other by 
the number of integral relations satisfied by the trial function. 

In this discussion we shall confine our attention to Kohn’s method and 
the interesting consequences recently established on this method in the case 
of zero energy. 

Consider the functional 


Pa { o*(v>,r)(H —E)W(v,,r)dr (7.1) 


where w and ¢ have the asymptotic forms 
ikr 


ws e 
wy a eikver 4: a (7.2) 


(7.3) 


(p ae eikver a. q' 
r 


and g and q’ are not necessarily the exact amplitudes. The functional 
vanishes for functions w which are exact solutions of the equation 
(H—E)) =0 (7.4) 


We shall vary the functions W, ~ so that in the asymptotic expression q, q’ 
change by 6g and dq’ respectively. The details of the calculations are given 
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in {17}. We shall not reproduce them here but merely state the final result: 
6Lf = 27 dq(V,, V2) (7.5) 
It follows that the expression 


1 
O=-—-—L£+q (7.6) 
2 


is stationary with respect to variations in the wave function, and that its 
value is equal to the amplitude g when exact wave functions are employed. 
Suppose now that instead of the exact wave function W we use a trial 
function yw, containing a few adjustable parameters and having the asymp- 
totic form 

ikr 


: e 
ikver 
Wy, ~e + 4, : 


(7.7) 


Let the value of the functional (7.1) calculated with the aid of this function be 
denoted by &,. 

The approximate expression for the scattering amplitude in Kohn’s 
method is taken to be the stationary value of (7.6) obtained using the function 


W,, 1.¢. 


I 
an=3'(- 5 2,444) (7.8) 
1 
In particular, if yw, includes the parameters c,, C2, ..., c,, the conditions that 
SY, should be stationary with respect to c,, C2,.-., Cy, J are [in view of (7.5)] 
OL 
—=0 (7.9) 
0c; 
OL 
‘= 27 (7.10) 
04, 


In the above discussion, we have only used the tact that the functional Q 
is stationary. The problem as to whether (7.8) yields the upper or lower 
bounds remains unresolved. In order to determine the nature of the 
approximation which can be achieved by the variational method, consider 
[18] the phase of the s-wave. Let y = fV 4a r, The function f(r) will be 
normalized so that it has the asymptotic form 


fr osin kr+Acoskr (7.11) 
in order to avoid complex quantities. If wy is the exact wave function, 
i = (tan 6)/k. If, on the other hand, w is the trial function, then according 
to Kohn 

1 i? @) 
~tand = 4-2] f(H-E)far (7.12) 
0 


where we must take the stationary value of the functional. 
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Let us take the trial function in the form 


f=9 + 2 Cid (7.13) 
where 
1 
oa sinkr+Aq,(r) cos kr (7.14) 


and @,(r) is a function which vanishes at r = 0 and is equal to unity when 
r — oo, but is arbitrary in all other respects; y; are functions which fall off 
rapidly as r — oo and ¢; are the constants to be varied. We then have 


J 

ptand =A+ ) cic;Mij+2 ) CR, +B (7.15) 
i,j i 
where 
Mi = My = 2] x(E-H)y,dr (7.16) 
0) 
R, = 2 { x(E-H)g dr (7.17) 
0 

B=2{ 9E-Hodr (7.18) 


0) 


Variation of the constants c; yields 
N 
2, Mycp=—Ri (7.19) 
= 


or in matrix form 
MC=-R (7.20) 


The matrix M7‘ may turn out to be singular when a determination is 
made of A and the other parameters entering into y;. The solution given 
by (7.20) will then not exist. Near these values the phase expressed as a 
function of A undergoes rapid changes and may be either much greater 
or much smaller than the true value. 

Specific calculations in a particular example [18] (see also § 16) show 
that as the number of terms in equation (7.13) is increased, there is an 
increase in the number of singularities in M~’, but that the phase fluctuations 
near these points are gradually smoothed out. 

However the limit of zero energy is a special case. It turns out that 
Kohn’s expression for the amplitude [equation (7.8)] gives the upper bound 
for the scattering length. The case when there is no bound state of the particle 
in the central force field and the case when the bound state does exist must 
be considered separately. 


2, METHODS FOR THE DETERMINATION OF WAVE FUNCTIONS 37 


Let w, be the trial function at zero energy which, in view of equation (7.7), 
has the asymptotic form 
Ay 


Pitas (7.21) 


r 


and consider the difference between w, and the exact function w: 


g=V,-wW (7.22) 
We shall use the following identity which can be verified directly 
| I 
a=a,t | tide ~ = | gta ar (7.23) 
2n ai 


where a is the exact scattering length. The quantity 
1 
ay =a,+ | W AW, dr (7.24) 


1s, according to equation (7.8), the Kohn scattering length if y, 1s taken to 
be the function giving the stationary value of ax. Consequently, 


1 
a =ag—=- | gig ar (7.25) 
Consider the integral 
ven | gHg dr ° (7.26) 


In the absence of bound states, H is a positive definite operator in the 
space of square-integrable functions. However, f is not square-integrable. 
We therefore introduce the integral 


J(A) =| 9 e-""Hge “dr (7.27) 


which includes the square-integrable function g exp(—Ar); for 2 = 0 this 
becomes identical with (7.26). The integral J(,) is such that JU) > 0. We 
shall show that J remains positive when 24 — 0. To prove this, it will be 
sufficient to show that J(/) is continuous at 4 = 0. 

Consider the difference 


JQ) —J(0) = | g(e~?*"—1)Hg dr+72 i g? eo 2*" dr—2) | gVge72** dr (7.28) 


Since [g(r)| has an upper limit for all r and Vg > 0 when r > o9, it can be 
shown that each of the three terms on the right-hand side of (7.28) becomes 
as small as desired when 1-0; this proves the continuity of J(A). 
Therefore the integral given by (7.26) is positive, or in the limiting case zero. 
It then follows from (7.25) that 


agdg (7.29) 
i.e. the Kohn scattering length is the upper-limit approximation for the 
scattering length. 


E.A.C. 4 
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Consider now the case when there are bound states. We shall assume 
for simplicity that there is only one such state. Let 


Hy, =| viHy, dr i,k — 1,2 (7.30) 


be the matrix elements of the energy operator H calculated with the aid 
of the orthogonal functions ¥, and w2, where wW, is an adequate trial function 
for the single bound state, so that H,, < 0. 


Let us reduce the matrix 
i a 
H,, A, 


( 0 

0 ES 

One of the eigenvalues, £,, will be negative and the other, £3, positive 
or zero, so that 


to the diagonal form 


E,E, <0 (7.31) 


Since the product £, F, of the eigenvalues is equal to the determinant of 
the matrix a sag it follows that 
Hy, fy, 

* H,,H1.—-Hi, <0 (7.32) 
We note that since ¥, enters twice into each of the terms in the inequality, 
the normalizing factor of Y, is unimportant. The inequality (7.32) remains 
valid when yw, is replaced by another function WY; which is not orthogonal 
to w,. This may be verified by substituting 


Wo =Ui-W [WW de (7.33) 
The inequality (7.32) then becomes 
H,,H5.-H3 <0 (7.34) 


where the dashed matrix elements represent the replacement of YW. by >. 
If we now take a suitable bound-state function for w, and let 


= —E, (7.35) 


where E, is the value of #/,, calculated using ¥,: and if we replace wy, by 
g, we have instead of (7.34) 


5 


] ~ 
| otigar <~=(| vs tor) (7.36) 


f 


Since Hg = H(y,—w) = Hw, (because Hy = 0) we may write 


~ | atigar <=( | vi tysdr) (7.37) 
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Then (7.25) becomes 


2 
a<agt =. (| WH, ir) (7.38) 


Therefore the right-hand side of (7.38) is the upper bound for the scattering 
length. 

It is important to note that (7.38) includes the trial function and the 
variational energy but it does not include the exact bound-state wave 
function and the exact value of the binding energy. 

The inequality given by (7.38) may be rewritten in the form 


as ax (7.39) 
where ax is Kohn’s scattering length calculated using the function 
| 
Viz Wear | Wi Hbede (7.40) 
{ 


The orthogonality condition [ ¥,y; dr = 0 is satisfied. It is evident from 
this that in Kohn’s method the trial function can reasonably be taken in 
a form containing the variational function of the bound state in addition 
to other expressions. 


8. THE METHOD OF TRIAL FUNCTIONS 


In the variational method the optimum parameters of the trial function 
are determined from the condition that the functional should be an extremum. 
However there are other methods for determining the optimum values of 
the trial-function parameters. One method is based on the use of certain 
relationships which are satisfied by the exact wave function (20, 21]. Depend- 
ing on which particular relationships are employed in the determination 
of the parameters, it 1s possible to modify the method in various ways. 
In the simplest case the trial function contains only one parameter. According 
to [20], the parameter may be determined from the integral equation for the 
scattering amplitude given by equation (2.25), It 1s convenient to choose 
the parameter so that it plays the role of the amplitude in the asymptotic 
expression for the trial function 


1 l i putt 
f~tsin (kr ~ 2) + qe z) (8.1) 


The equation for q, is then of the form 


q=~-2 i uw, Vf, dr (8.2) 
0 
If the trial function F, depends on the n parameters c,, c,...., ¢, In addition 
to g,, and has the asymptotic form 
l l i fc 
fi~ ), Sit (11 — 5 + q,e (+ 2) (8.3) 
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then one of the possible methods of determining these parameters is as 
follows. Substitute f(c,,¢3,...,C,,9; 7) into the right-hand side of the 
integral equation (2.23). This yields the function 

GC 15 C25 0063 Cae Ip *) = u(r) +2 | fo(nu(r)—u(nor) Vr yfi(r) dr’ (8.4) 

0 

Since f; is not an exact solution of (2.23), it follows that g, 4 f;. We expand 
g,and f; into series in powers of r and require that the first m terms be identical. 
This yields a system of equations for the coefficients c,. In order to determine 
the amplitude g,, we substitute f; into the integral relation for the ampli- 
tude (8.2). 

The application of one of the variants of the method of trial functions 
to the scattering of electrons by atoms has been described in a number of 
papers [21]. 

The disadvantage of the method of trial functions is the absence of a 
criterion of the type given by (7.29) which would determine the nature of the 
approximation in the limiting case of zero energy. Moreover, determination 
of the coefficients c, by means of the expansion of the trial function f,(r) in 
powers of r becomes ineffective when the potential well is so deep that the 
wave function has a number of zeros inside the well. 


Chapter 3 


COLLISIONS OF ELECTRONS WITH HYDROGEN ATOMS: 
GENERAL THEORY 


9. BAsIC RELATIONSHIPS 


THE collision of an electron with a hydrogen atom is one of the simplest 
processes, particularly when the energy of the incident electron is insufficient 
to excite the atom and the scattering 1s elastic. Nevertheless, even the 
description of elastic scattering of electrons by hydrogen atoms 1s theoretically 
quite a difficult problem. 

We now consider the main factors which characterize the scattering of 
a Slow electron by an atom. Suppose for the sake of simplicity that the 
electron shell of the atom is “‘frozen”’ so that it does not exhibit any changes 
during the interaction. The effect of the atom on an electron can then be 
determined by considering an equivalent centre of force. The interaction 
between the electron and the atom is in fact found to be attractive in 
character. 

It sometimes happens that when the atom 1s regarded as a centre of 
force there exists a shallow zero energy level for the electron, i.e. a negative 
ion can be formed, or that there is a virtual level. As a result, low-energy 
scattering approaches resonance scattering. Even this single fact complicates 
quantitative discussion of the problem quite considerably (for example, 
it prevents us from using the perturbation theory). 

When the electron shell of the atom is “‘unfrozen’’ a new factor arises, 
namely, polarization of the atom by the scattered electron. At large distances 
between the electron and the atom the influence of polarization may be 
described by an interaction of the form U = —f/r*, where f is the polariza- 
bility of the atom. This expression represents the interaction between the 
point charge and the induced dipole. However, at small distances, a single 
coefficient f is insufficient to describe the polarization effect. In contrast 
to the static polarization of an atom in a homogeneous field, we must now 
deal with the more complicated dynamic polarization. Moreover, the 
indistinguishability of the electrons has a major effect on the various detailed 
features of the scattering process. As a result, the so-called exchange- 
scattering, in which the incident electron changes places with the atomic 
electrons, becomes indistinguishable from simple scattering and both must 
be regarded as a single process. Mathematically, the indistinguishability 
of the electrons is reflected in the symmetry properties of the wave function. 
When these symmetry properties are taken into account, it is found that the 
effective interaction between the electron and the atom becomes non-local, 

41 


42 THE THEORY OF ELECTRON-ATOM COLLISIONS 


and the equations for the wave function F(r) of the scattered electron include 
an exchange operator of the form 


WF = { W(r',r) F(r’) dr’ 


in addition to the central-force term VF. 

When excitation is possible the theory becomes even more complicated. 
In elastic scattering calculations it is necessary to solve only one equation 
when polarization is neglected. When excitation is possible the corresponding 
calculation (on the same approximation) requires the solution of at least 
two equations even when there is only one inelastic scattering channel. 
When there are several channels it fs necessary to consider a complicated 
set of equations, since the cross-sections for processes in the different channels 
are related, in view of the conservation of the number of particles, which 
is expressed mathematically by the unitary property. It follows that the 
elastic scattering 1s also modified under the influence of inelastic scattering. 

The effect of all these factors on the collision between an electron and 
atom, their relative importance, and the connection between them, have 
not as yet been adequately explored, although it is possible to write down 
the equations which formally represent the exact theory including all the 
above factors. 

We now consider the theory of time-dependent states, which is analogous 
to the theory in § 1.+ The energy operator is of the form 


H = Ao(ry)+Holr2) + Volri) + Volr2) + V(r1, v2) (9.1) 
where H,(r,) and H,(r.) are kinetic energy operators for the first and 
second electrons respectively, Vo(r,) and V(r.) are the energies of inter- 
action between the electrons and the nucleus, and V(r,, r,) is the energy of 
interaction between the electrons: 


Vir)=-—, Vr=-—,  Venr)=———) 0.2) 

ry r2 Ir oe r'>| 
The wave function for the electron + atom system depends on the coordi- 
nates and the spin variables of the electrons. In the non-relativistic theory 
the total spin of the system is a constant of motion, so that we can separate 
out the spin variables and confine our attention to the coordinate part of 
the wave function. Since the total wave function is antisymmetric with 
respect to the interchange of coordinates and spins, there is a connection 
between the symmetry properties of the coordinate part of the wave function 
and the magnitude of the total spin of the system. 

In the case of the two-electron system which we are considering here, 
there are two possible spin states, namely, a triplet and a singlet state. In 
the triplet state the coordinate function is antisymmetric while in the singlet 
state it 1s Symmetric. 


+ The ensuing analysis is a variant of the so-called theory of rearrangement collisions 
for two particles in the field of a fixed centre of force. It will be assumed that the two 
particles do not form a bound state. The various problems in the theory of rearrangement 
collisions are discussed in [22~25]. 
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Let M(r,,#,,¢) be the coordinate part of the wave function satisfying 

the equation 

o® 

at 
As has been pointed out above, the function ® should be symmetric or 
antisymmetric with respect to the interchange of r, and r,. However, it is 
in many ways more convenient to begin without assigning symmetry 
properties to ® and to symmetrize the result in the final stages of the 
analysis. 

Proceeding as in § 1, we shall express ® on the interaction representation, 
that is, we expand ® in terms of the eigenfunctions for a non-interacting 
electron-atom system. Two such systems can be formed. One is described 
by the operator 


H® 


H, = Ho(ry)+ Volr1) + Holr2) (9.3) 
and the eigenfunction 
QL (ry, 2) = kF(2n)7F eh" WC) (9.4) 
where yw, is the atomic function. The subscript « represents the set of 
quantum numbers which characterize the state of the atomic electron, 1.e. 
the energy E£,, the angular momentum, and the component of the angular 
momentum in a special direction (in the present case, the direction of 
motion of the incident electron). The plane wave exp (/kv.r,) describes 
the free motion of the electron with energy E = 3k? in the direction of v. 
The functions g{’ are normalized to d(E~—E’)d(v—v’)dga and satisfy 
the equation 
Hy, py? = &oy” (9.5) 
where 
S=Htk*+E, (9.6) 
The eigenvalues & of the operator H, form the continuous spectrum 
beginning with the negative quantity —|F |, which is equal to the ground- 
state energy of the hydrogen atom. 
The other system is described by the operator 


H, = Ho(ry) + Holr2) + Volr2) (9.7) 

and the eigenfunction 
QO(ry,12) = KF(2n) 77 es, (r2) (9.8) 
whose significance is analogous to equation (9.4), except that r, has been 
replaced by r,. Either the functions given by (9.4) or those given by (9.8) 


can be used as the basic functions. Accordingly, we shall consider the 
following two equivalent expansions for ®(r,, r,, f): 


O(r,,r2,0 => [YOC, vy, te Qn, 7.) db dQ (9.9) 


Pr r,0=) | POE y te! oP Ur, 12) dé dO (9.10) 


a 
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As in § 1, we shall assume the existence of the limits 
yO = lim POE, v, t) r= 12 (9.11) 


t~+o 


and determine the S-matrix from the relation 


Oe =>, | S*(6,v, 0, 6, v’, w'PO.(6', wv) de'dQ’ i,t’ =1,2 (9.12) 


Let us express ® in the Heisenberg representation. Let (r,,r2, 6, v) 
be the eigenfunction of the operator H corresponding to the eigenvalue &, 
so that 


Hp =6€@ (9.13) 


where v is a unit vector which enters into the argument of (as in § 1). 
The eigenvalues & of the operator H form a continuous spectrum lying 
in the same region as the continuous spectrum of H, and H,, that is, 
—|E,| < & < o. Moreover, the region & < —|E,| includes a discrete level 
of the negative hydrogen ion, but this level is of no significance in the present 
context. 

In order to determine the function g, we must also specify its asymptotic 
behaviour. We shall suppose that an electron of energy 3X9 is incident on 
an atom in the state a, the position vectors of the incident and atomic 
electrons being respectively r, and r,. Accordingly, the function o will 
asymptotically tend to o{(r,,r,) aS rz > 0. Therefore the function @ 
depends on the parameters defining the initial state of the atom, and this 
will be implied in the ensuing analysis. 

Let us write ® in the form 


O(r,,r,1) = | eS WE v) o(r,,72, 8, de dO (9.14) 


The relation between ‘P“ and ¥ can be established by comparing equation 
(9.14) with (9.9) or (9.10), and is given by 


POE», 1) = e~F - Dt WO’ y &', yw") P(E’, vd" dQ’ (9.15) 


where 
WS y,8',v') =| p*(r,,173,86,v) O(',,1°2,6,v') dr, dry i=1,2 (9.16) 
Let us transform equation (9.16) by substituting 


P(r 1,12, 65) = OL +x (9.17) 
where, in accordance with (9.6), 
6’ = tk5t+E,, 
We first consider the case i = 1. Since 
(H—-&')o =0 (9.18) 
and 
(1, —& )gl) = 0 (9.19) 
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we obtain 
(H —6')y = —(H—-H,)o (9.20) 
or using the explicit form of H and H, 
(H—8')y = —[Vo(r2)+ V(r, rr) |p (9.21) 


By analogy with the method used in § 1, let us define the resolvents 
G=(H—6)', G,=(H,-6)7',  i=1,2,... 

as the limits 

G(é +i0) = lim (H—€— ie)! 

Z leis (9.22) 

G(&6 +i0) = lim (H,—6—ie)7' 

e> +0 
With the aid of equation (9.22) we have from (9.21) 
Y= Jim oS" + ie) [Vo(r)+Vr1, 172) Joo (9.23) 


On substituting this fate é 17), 
p= - jim oC + ie) [Vo(r2) + V(r, 2) oh (9.24) 


Again, following the alia al Chapter 1, the latter relation may be written 
in the form 


po =— lim ieG(6' + ie)? (9.25) 
r- +0 
The equation for the resolvent 
G = Gy—G,[Vo(r2)+ V(r. 72) IG (9.26) 
can now be used to transform (9.24) so that it reads 
1 
Q = Day’ — 26210 [Yo(r2)+ V(r1, 12) Io (9.27) 


This result can be used to reduce the expression for w{ to the form 
TLE, v, 6", V9) 


WB, ¥, 6, ¥9) = (6 — &) H(v— v9) 6 a 


.28 
aao Pam & a eee (9 ) 


where 
TUE, ¥, 6", 9) 
=| or. 82.6 V)[Vo(r2) + V1, 12) ] 0011, 2,6, Vo) dry dr, (9.29) 


Substituting equation (9.28) into (9.15) and repeating the steps analogous 
to those given in § 1 we obtain 
PDE, v) = VE, v)d,,, 


POE v) = PONE, v)—2ni { TU (8. 8, Vo) PLAE, v9) dQ’ 


aoa 


(9.30) 


The S-matrix element is therefore given by 
S'(8,v, a, 8", ¥9, 09) = 56’ —&)[5(vp — v) 5,,,— 221 T, OE, ¥,8,V9)] (9.31) 


aga 
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The presence of the factor 6(é’—6&) in (9.31) ensures the conservation of 
energy (6 = 6’). In accordance with equation (9.6), the energy of the final 
state is © = 4k7+4 E,, and therefore 


Ak7 +E, =4k$+E,, (9.32) 


Consider now the case / = 2 and suppose that the function ¢ is, as before, 
defined by equation (9.24). The vector V%” describes the so-called exchange 
scattering. It is now convenient to express the resolvent. G in (9.25) in terms 
of the resolvent G, as in (9.26) rather than in terms of G,: 


G= G, — G,[Vo(r))+ Vr, r,)|G (9.33) 
Substituting (9.33) into (9.25) we obtain 
yp =— lim ieG, eg? — lim G,[V(r,)+ V4, 72) ] 9 (9.34) 
e>+0 E> +0 
We shall show that 
lim 6G,(& + ial =0 (9.35) 


To do this let us expand g{)(r,, r.) in terms of yY: 


gay { a(8')o? dé’ dQ’ (9.36) 


If the atomic function w,(r,) entering into ~\” refers to the discrete spectrum, 
then a,(&’) is a continuous function of &’. But if the atomic function P(r) 
refers to the continuous spectrum, then the coefficient a,(é’) has a 6-function 
type of singularity. However for this part of the spectrum the symbol 
), represents integration over the continuous spectrum with respect to «, 
which removes the singularity and leads to a finite expression. Let us operate 
with G,(6+ie) = (H,—€—ie)~' on both sides of equation (9.36). The 
result is 


Goh) = = Efe Coe. = Oa de QQ’ 
The integral 
16 : 
, = don | dQ’ a AC ‘\opl?? 
has a finite value, so that lim eJ = 0 (we note that if we evaluate G,yY 
e-Q 
instead of G, ef), the result becomes G,~%) =—g@/ie, so that 
lim ieG,g& = pg). In view of equation (9.35) we have trom (9.34) 
e> +O 
\ 
ae Hr ol + Viry, ra) ly (9,37) 


It follows that gy is the solution of the nonhomogeneous equation with the 
resolvent G, and the solution of the homogeneous equation with the 
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resolvent G,.f Substitution of equation (9.37) into (9.16) yields 


T!2(8, v, 8’, Vo) 


CUE VE Vi) Sa 9.38 
wi(8, Vv, 6, Vo) ae. (9.38) 


where 
T2(8, v, &', Vo) 


= [Mr re 8 Molt Vr DIO Evo) drs dry 0.39) 


Let uS now substitute (9.38) into (9.15) and take r to +o. Since according 
to equation (9.30) 
US vie PS) (Ey) 
we have for discrete « and a 
P2(8,v) =0 


P28 y) = —2ni T28,v, 8", Vo) PY (8, V9) dO 


aoz — a 


(9.40) 


For continuous «, corresponding to integration over the continuous spectrum 
of the atomic states in (9.10), the element 7)2? of the 7-matrix regarded 
as a function of « has a singularity of the form 6(E,—4k$). Consequently, 
for continuous a, pg, v) will also have a singularity. We shall not 
consider the nature of the singularity here. 

Substituting equation (9.40) into (9.12) we obtain the following expression 
for the S-matrix element: 


S216, 0,0, 6", Vos Uo) = — 2215(8’ — €)T.26, v, &, Vo) (9.41) 


aoa 


In order to introduce the symmetry properties of the wave functions let us 
construct the symmetric and antisymmetric combinations 


I 


Y= 
: J2 


CPO) 4. yO) (9.42) 


| 
t= PAG ae be (9.43) 
V2 


In view of equations (9.9) and (9.10) it is evident that the state vectors given 
by the latter two expressions correspond to the symmetric and antisymmetric 
functions @. 

Transitions from the initial to the final state for ‘YW, and ‘¥, are described 
by the symmetrized S-matrix which can be constructed from S!! and S$?! 
by analogy with (9.42) and (9.43): 


St =(S'! 452) (9.44) 
S7 =(S!— $21) (9.45) 


+ This point is discussed on the basis of the time-independent theory in [26]. 
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From (9.41) and (9.31) we have 
S*(E,v, 0,8", Vo, %) = (P—6') [5(V9 — V) 5 ga — 271 T(E, V,E',Vo)] (9.46) 


where 
10 = Te an (Bes (9.47) 


The same result can be obtained in another way. Instead of the functions 
gs) and ~ we can use as the basis in the interaction representation the 


Rs and antisymmetric combinations (op) + a?) //2 normalized to 
6(€ —&')6(v—v’), and introduce the symmetrized function 


04 =— lim icG(6 + ie)(O) +92 (9.48) 
7 e> +0 


which may be rewritten in the form 


(1+P,,)9 (9.49) 


Pi = - 
where P,, is the exchange operator. Consider the integral 
1 
/9 | (pi + OS)" 4 dry dr, 
Using equation (9.49) we have 


“5 | OW HOY es drs dr, aa 5 | ae Prdetl Prado dr dr, (9.50) 
2 


However, 


| (LEP,,)e0 *Q+P,,)gdr,dr, = [oQ+P, 2) p* dr, dr, 


=2| PULP, 2)9Q™ dr ,dr, (9.51) 


and therefore 


=| we piy* p,dr, dr, = | (of? + oY)” pdr, dr, 


= | tod, dr, + | 92 ode, dr, (9.52) 


This corresponds to the representation of the matrix 73, in the form given 
by (9.47). It is possible to establish a relation between 75, and the effective 
cross-section by analogy with $1. But it must be remembered that the 
symmetrized wave functions correspond to the superposition of two physi- 
cally identical situations: scattering of the first electron by the atom con- 
taining the second electron, and scattcring of the second electron by the 
atom containing the first electron. The expression for the diffcrential cross- 
section for elastic scattering 1s of the form 
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= (nyt 21Q (9.53) 


do 


OC Fas ké Le 2920 


while the expression for the diffcrential cross-section for inelastic scattering 
in which the atom undergoes a — from the state a, to the state « 1s 
dQ (9.54) 


a (2n)* 3 | Taoal” 


doz, = 
When comparing theory with caeleians it must be remembered that in 
practice one 1s usually dealing with the scattering of an unpolarized beam 
of electrons forming a statistical mixture. Accordingly, equation (9.54) 
must be averaged over the triplet and singlet states. The result is 


(d6yoa> = 2do,,, +4do;,, (9.55) 
There is an important relation between the elastic and inelastic scattering 
cross-sections which is a consequence of the unitarity of the S-matrix and 


forms a generalization of the optical theorem given by (1.48). The matrices 
S* and S~ are unitary, i.e. they obey a relation analogous to equation (1.7): 


3 { S**(6,a,v, 6’, a’, v)S*(6,0,¥,8", 0", v") de dQ 
= 58-8") S(v"—v') 5yrqr (9.56) 


With the aid of (9.46) we can use the unitary property to deduce a relation 
analogous to (1.48), namely 


—Im T,,,(0) = x y | \T.2,|? dQ t (9.57) 


In this formula the summation is carried out only over those levels which 
can be excited at a given energy in accordance with equation (9.32). When 
the initial state of the atoms is the ground state and the energy of the incident 
electron 1s insufficient for excitation of the atom, the sum in (9.57) consists 
of one tcrm only. 

The S-matrix is also invariant with respect to time reversal. This leads 
to a generalization of equation (1.8) and may be written in the form 


SioalO,¥, 8", Vo) = Saal’, — Vo. 6, —V) (9.58) 
Accordingly, we have the following relation for the 7-matrix elements 
TE, VY, 6, Vo) = T,3,(é, —Vo; 6, aes v) (9.59) 


It is understood that in the time-reversed transition described by the right- 
hand sides of (9.58) and (9.59), the signs of the orbital and spin angular 
momenta have been reversed. The transition of an atom from one state 
to another which we have been considering is, in fact, invariant with respect 
to reversal of all the components of the angular and linear momenta. It 
follows that we can use equations (9.59) and (9.54) to write down the 
detailed-balance relation: 

kode. =k d0,;, (9.60) 
(S, Vo, &, Po). 


t T+ (0) represents T+ 


2020 2020 
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10. EQUATIONS FOR THE WAVE FUNCTION OF THE 
SCATTERED ELECTRON 


If the eigenfunction g(r,,1r,) of the energy operator for the atom + 
electron system is expanded in terms of the atomic functions @,(r,), the set 
of expansion coefficients will represent the wave function of the scattered 
electron. 

In this section we shall establish the equations and asymptotic behaviour 
for the wave function of the scattered electron. 

To start with, let us change the normalization by multiplying g and g, 
by (2z)?k9?. The corresponding functions will be denoted by w& and w,. 
For the & functions we have the integral equation 

W(r,,ry) = — lim ieG(& + ie) evo"? Wi ("1) (10.1) 
e+ +0 
As has been shown above, from equation (10.1) can be obtained two integral 
equations (one through the resolvent G, and the other through G,): 


W(r1,12) = eer Wao 1) = G\(E+ i0) LVo(r2) + V(r, |p (10.2) 
Wr r= —G (+ i0)[Vo(r,) + V(ry, ry) | (10.3) 
or, using Green's functions 


W(r 1,12) = ekorer2 Walt) — 
— [GG yrar i ror dt Vir DW.) dry dr, (10.4) 


* 


Wrist) =— [ Garasra rier) Vor +V DW rpdridr, (10.5) 


» 


The latter equations may be used to investigate the asymptotic behaviour 
of Wir,, r2), as follows. Consider the asymptotic behaviour of the Green's 
functions G, and G,. These two functions are expressed in terms of the 
free-motion function Gy. By expanding G, and G, in terms of the eigen- 
functions of H, and Hy, it is easy to show that 


Gary roti. ty) = > War Wr(r )Golkas 2.12) (10.6) 
and 
Go(ry.%2,11,%2) = > Warr Wa (1 )Go( Kas M514) (10.7) 
where 
Gk, rer) == | etkalr—r'| 
2n|r—r'| (10.8) 


k, = J2(E,,-E,) +k3 
Whenr > r’ 


Go eee gikar eT Kane (10.9) 


For those values of « for which F,,+ 345 < £,. the function Go(r) falls off 
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as exp (—|k,|r). Substituting equation (10.9) into (10.6) and (10.7) and 
retaining only those terms in the sum in which k, 1S real, we have 


G~ = Tver te) — 


r2~> @ os kz 2>0 


Substituting equation (10.10) into (10.4) and ak 
gikore eikar2 
y (ry, r>) es (« NOSE een = Wa('1) 2 > nae W(r;) (10.11) 
ro 


r27> 0 (k2> 0) 


where 


] —ik nurs’ , r , / ; , / / , 
ian == a een wiry) Vo(ra) + V(ri ro) wry, ro) dr} dry (10.12) 


gree (10.10) 


When r, > ©, 


tke 


y (11,12) ~ y Dag — Wal2) (10.13) 


where 
I —ik,n.r,’ / ! rf , ! / / ' 
d= — 5 | € are WF) LVoiri) + V(r), 1) Wnt, 73) ] dr, dr (10.14) 


Comparison of equations (10.12), (10.14) with (9.29) and (9.39) yields the 
relation between f,,2, Zao2 and (ae j Bee 


TOE, i, é, Vo) = Qn (2x)? ° ka k3 es aeals Yo) (10.15) 
! 
Tyoa (8,2, E, Vo) = — a kG KF gagalM, Vo) (10.16) 
) (27) Q 
Let us now expand yw(r,, r,) in terms of the atomic functions: 
Wr1,%2) = >; W(r1)F (12) (10.17) 


where the symbol ) includes integration over the continuous spectrum. 


We substitute equations (10.17) and (10.6) into (10.4), multiply by wz and 
integrate. For the coefficients F,(r,) which represent the wave function of 
the scattered electron we then have the set of equations 


F(a) = ef" 5. | WIC) Vora) + VOL 1] 
Wal )Go(Kas Pas 0)Falr) dry dr’, (10.18) 


where the integration is with respect to r;. Let 


Via (12) =| Wi (r1)[Vo(r2)+ V(ry, ro) \Wa(r;) dr, (10.19) 
Using the explicit form of V i and V(r,,r,) given by (9.2), we obtain 
V,A(r>) oe oe jt [ Perea (ri )Wa (ry) a (10.20) 
Ir —r,| 


$2 THE THEORY OF ELECTRON-ATOM COLLISIONS 


Substituting (10.20) into (10.18), we obtain 
F(r) = 5... [ Gok rer Veer Fr)dr’ (10.21) 


where the subscript 2 has been omitted. This system of integral equations is 
equivalent to the set of differential equations 


AF, +(k?-—2V,,)F, = 2 00 Vay Fy (10.22) 
a’ #2 
with the asymptotic condition 
(kyr 
; ; e'"* 
— pi sovo:t Op. + Jaca r37 ow (10.23) 


The system given by (10.22) can also be obtained by substituting the expan- 
sion given by (10.17) into 


[¥iC)H-SOr r2)dr, = 0 (10.24) 


It is shown in [17] that the latter relation is a consequence of the fact that 
the functional 


a) [ VIC DH- Sdn Fara) dr dra (10.25) 


is stationary with respect to arbitrary variations of F,, which vanish when 
r, > oo. On this basis we shall suppose that equation (10.24) is valid even 
when wW(r,, 72) 1S not an exact solution of (H—¢)W = 0 but represents an 
approximate expression, for example, the sum of a finite number of terms 
in (10.17). This leads to the equations for F, which give the best approxi- 
mation in the variational sense. 

Consider the part of the sum in (10.17) which corresponds to the con- 
tinuous spectrum. It 1s given by the integral 


Hrs r2) = | Fs, ky ra), 14) de (10.26) 

There is also a relation between A and x which follows from the law of 
conservation of energy 

k* +x? = ko +2E,, (10.27) 


It turns out that F(«, &, r) in (10.24) has a singularity at x? = &$ [27]. Let 
us consider the nature of this singularity in greater detail and expand F, 
in terms of w: 


Ko) Can (10.28) 
Substitution of this expansion into (10.22) yields 
Yd Caa(A tk? —2V,.)0e = 2 ¥ cg, Vag hy (10.29) 
a’ A, BF 


The function w,- satisfies the equation 


Aw, = —2 (E, i -) Wi (10.30) 
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From the latter equation and from equation (10.21) we have 


1 
by Cax(ke —2E Wy =2 » (a) Cpa WY; (10.31) 
a’ aa\lt—F'|/ pe 
Multiplying wz, and integrating, 
l 
caulk =2E,) = 20 ers (a) Ws (10.32) 
p, a PO] pa 


Substituting « = @ and using (10.27), 
ca (KB—HA) = 2D ey VOWEL Wed — mont dry (10.33) 
em) 


In order to asa iaatieetad connected with the specific asymptotic 
behaviour of continuous-spectrum Coulomb functions, we shall cut off 
the Coulomb interaction by introducing the factor exp (—wr) where @ 1s 
a finite positive number which can be as small as desired. The right-hand 
side of (10.33) will then be finite for Ay = x, so that 


const 


Cal) = Bae (10.34) 
Therefore, F(x, r) can be written in the form 
F(x, r) = rae Wr) + Br) (10.35) 


Where A(r) is a function having no Teale tke = kK: 

When we substitute equation (10.35) into (10.26) we must determine the 
path of integration with respect to k. This can be found from the asymptotic 
form of W(r,,r,) given by equation (10.15) when r, > co. We substitute 
(10.35) into (10.26) and separate the angular variables. Assuming that r, 
is large, the radial part W,(k, r.) has the asymptotic form 


| l 
Wir) ~ : sin («0 -5 + 5) (10.36) 
so that 
l 
t 2 sin(xry — 5 +41) 
cons 
Jirist2) ~ Maal 2) |e K?dk (10.37) 
can be reduced to the form 
a } re (10.38) 
~ con r,)— Cdk 
ae ry k2—K? 


The path of integration can now be chosen so that the pole k = ko 1s excluded 

by taking a small circle around it in the upper half-plane, while the pole 

Ko = —ky is excluded using a circle in the lower half-plane of complex x. 
E.A.C. 5 
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The result is 


eikors 


(10.39) 


J ~ const W,,(r2) 
My 


In the absence of a purely Coulomb field the divergence will be of the form 
(k§—K?)~1*4« instead of (kg—Kx?)71, and instead of the asymptotic form 
(10.36) it will have the Coulomb form: 


1 l l j 
Wy =-sin er =F 4 In dur + arg (I+ 1 -*)] 
r 2 «K K 


The analysis is then very much more complicated and has not in fact been 
completed. However, the general result will be similar to that in the simpler 
case of a truncated Coulomb interaction. 

Let us now consider symmetrized function w,. As with equation (9.49), 


. l 
this function can be obtained by applying the operator ri + P,2) to 
the function wy 


1 i 
Y a Ba Ra (10.40) 


We now deduce the integral equation for ¥,, transforming the first compo- 
nent on the right-hand side of equation (10.40) using the inhomogeneous 
equation (10.2), whilst the second component will be transformed with the 
aid of the homogeneous equation (10.3). Since P,,G, = G, and 


Py V o(r2) = Vo(r1)s 
1 ‘g 
= MOA Cr) — | Gris 25 FP: 
Ws V2 1) Wh 15 "2.01 >) 
x [Vo(r3) + Viri, ro) JW a(r.ro) dr, dr, (10.41) 


Using the expansion (10.6) for the Green’s function, 


1. 
ba = eer alr) — Daley | y(r') 


N a 
x [Vo(r2)+ Ver. r3)|Golka ta taWelri.re)dry dr, (10.42) 
from which we obtain the asymptotic behaviour of w, for large r,: 


1 »tkor2 oikar. 
Ws = Wa (11) (=: ai 1! 7 a ny GW) = (10.43) 
/2 Si Dra’ eas \2r, 


[2 
\ aoa | n.r’> , , i , ? ’ , ./ .f 
aca = — 2 fe camer (r LMo(ra) + VO rar rg) dry dry (10.44) 
l 
Writing Ww, in the form Ww, = Fa (Wry, ro) & Wrz, ",)], and bearing in 


mind (10.13) and (10.15), we have 
Gaoa = Styx £ Geoa (10.45) 
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We might expand wy in terms of the atomic functions, writing it in the form of 
(10.19). However, in this form the symmetry properties of w, are difficult 
to see (they are ensured by that part of the expansion which refers to the 
continuous spectrum). Instead we writc yw, in a form in which the symmetry 
properties are immediately evident: 


i 
Ws = 5 Dd LPs (raWalry) £ Fe (ri alr2)] (10.46) 
\ a 


This representation has an important property: it is invariant with rcspect 
to the transformation 


Fears) D Dias 


where Dj is symmetric or antisymmetric, depending on which of the two 
signs 1s taken: 

Dip = + Dja 
This result may be used to remove the singularity in F;(r) in the integral 
over the continuous spectrum which is contained in (10.46) (this does not, 
howevcr, definc F* unambiguously). We now proceed to derive the equations 
for FZ. Substituting equation (10.46) into (10.42), multiplying by W2(r,) 
and integrating with respect to r,, we obtain 


Fy (r2) =F 3 W (12) | ry )Fe(ry) dr, tere”? Ong 
< Df WR Vor) + Vr ra) Wal FE (0) 
x Golkas tata) dri dra FD | wry) (Vora) + V(r, 79) 


x Wa) FEC )Go(ka t2,62) dr dr, (10.47) 


This equation can be transformed to a more convenient form. Consider 
the following part of the last term on the right-hand side of equation (10.47): 


J=FY | GColkes ras raVolrWelra) drs | WIFE) dr, (10.48) 


Using equation 
~~ tA, Wa: ate Vo(r2)Wa = Ey Wa 
we can rewrite (10.48) in the form 


J=FL | GA +E Wel s)Golkastasradrs | UM )FEr)dr, (10.49) 


and hence, bearing in mind the explicit form of Go given by (2.7), we have 
after some rearrangement 


J=£Y [Welrs)+ [GK -EGolka tata) dra} [VIO FE@i)dr, (10.50) 


The term 
wir )Fe(ry) dry 


+ ry Wrz) | 
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is of exactly the same form as the first term on the right-hand side of equa- 
tion (10.47) but with the opposite sign; therefore they cancel out. It follows 
that we can reduce (10.47) to the form 


Fi (r2) Set Saao— >, | Golka Vo; P53) Vaal) Fz (13) dr 
~Y [ Golkys ro. ra)dry | WAC r)FE) dr, (10.51) 


where V,,- is given by (10.20) and ¥,; represents the kernel of the exchange 
operator 
Wray, 02) = +We(r (Vr, ra) + Ey —tky Wer) (10.52) 
On reversing the order of integration in (10.51) and substituting 
KE 954) = — Vaal?) Gol kas ¥252)— | WEG, 1)Golke ra ri)dr, (10.53) 


e 


we have finally 
FE(r2) = e072 5,0, + [ Kiera, ta) FE(14) dry (10.54) 


The integral equations (10.51) and (10.54) are equivalent to the following 
system of integro-differential equations 


AF# +(k?—2V,,)F+—2 i Wir’, r)F,(r') dr’ 
=2y (Vie F2 + | Ww, Ft dr’) (10.55) 


a’ xa 
with the asymptotic conditions 
ikar 


Fi(r) pinata Saxo 7 qi (Nn, Vo) 


ro Fr 


(10.56) 


In the ensuing analysis we shall, in view of equation (10.55), look upon 
the term V FE + f W3i.F dr’ for a # a’ as expressing the relation 
between the states o and a’. 

Utilizing equations (10.54) and (10.16) we have 


Ea, n, G, Vo) —— 


I 
Substitution of (10.57) into (9.54) yields the following expression for the 
inelastic scattering cross-section 


k 
dois = |Mine|” dO (10.58) 
0 


The differential cross-section for elastic scattering can be obtained from 
this by substituting « = a: 


ioxo = |For] IQ (10.59) 
These expressions can also be obtained directly from the asymptotic 
form (10.56) by dividing the flux in the solid angle dQ (ie. k,!qia|? dQ), by 


the flux density in the incident wave ko. 


do 
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On substituting (10.57) into the unitarity condition (9.57), 


4 k 
=F 1m gt,(0) = Y=! | g#,|? dO (10.60) 
ko 7 Ks 
The quantity 
k 
= -Yoee-DE| g2,|? dO (10.61) 


is the total cross-section for all the elastic and inelastic processes possible 
at the given energy. Therefore, 


4nImqi,,(0) = ko* (10.62) 


11. EQUATIONS FOR THE PARTIAL WAVES 


Let us separate the angular variables in equations (10.55) and (10.56) 
and expand FZ in terms of the spherical harmonics as in § 2. In order to 
distinguish the quantum numbers associated with the angular momenta 
of the atomic and incident electrons, we write /,, 4, for the atomic and 
l,, ft for the incident electrons, respectively. The symbol n represents the 
principal quantum number and the symbols + will be omitted. With all 
the indices written out in full, the expansion is 


Tota 


Frit) =, V4n(2l, + Di? ¥,,,(0, 9) ai) 
1 r 


The summation over fl, iS unnecessary in view of the conservation of the 

z-component of the angular momentum of the atom + electron system: 
ty ti, = const 

Since the z-component of the angular momentum of the incident clectron 

in the initial state is 3 = 0 it follows that h,+p5 = nf. If we now sub- 

stitute equation (11.1) into (10.56), multiply by Y“2o(0, y) and integrate, 


we obtain the following asymptotic expression for f;742.: 
. lon 
I On 0 el (kar ~“$*) 
ile. Ay = Sea fe eee 0lnw2 FN 2) 
nly pty ak ke sin Kno? os 9) Ong Ot, 19 Oy sp? a 4 igi pOnt sp (11.2) 
no r 
wheré 
SOlan2 : ok 
Unol wont yw ra | Y 130 Ta aria Mos 2) Lives dQ (1 1.3) 


Next, let us substitute equation (11.1) into (10.55), multiply by Y7;%,.-(0, ~) 
and integrate with respect to the angles. It will also be convenient to 
substitute 


Vow ie =e [ Yi5n348; () Vig") re 'y19(O, ~) AQ (1 1 .4) 
and 
yf twas ne rp) = rye { Yi ,(0', @') Wer’ Ye, (0, 0)d2dQ’ (11.5) 


In view of the presence of [r—r’|~' in W,,-, the form of (11.5) will be 
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Ul 


different forr <r’ andr2vr’. In order to take this into account, we shall 
use the symbol W“ for r < r’ andW™ for r > r’. The system of equations 
(10.55) will then be of the form 

Cohan + (ia ita 0) in 


I 
dr2 rH 


le 8) 

=2 3 [vin Sa + | ORE Seteedr | (11.6) 
n’l,’la’ py’ 6 

where f,742(0) = 0 when r = 0. 

The function f,74? describes the state of both electrons on the /,, j,. 
!,, 4, Tepresentation. It is convenient to replace this by the /,.1,.L, M 
representation, where L is the total angular momentum and 4M its z-compo- 
nent (this is known as the total angular momentum representation [28)). 
Since in the nonrelativistic theory the angular momentum is a constant of 
motion, on the total angular momentum representation the system of 
equations (11.6) splits into a series of unconnected sub-systems. The transi- 
tion from one representation to the other is carried out with the aid of 
the transformation 
en = » CLM fipt (11.7) 


S ntits lant nly uy 
Mi. HH? 


where Cyiv\1,,, are the Clebsch-Gordon coefficients. Substituting equation 
(11.7) into the asymptotic expression (11.6), 


1 Lo x 
LM : vee 
nil, ~ i sin Knot — “> Cie p01090 9 nno oe 
no 
jie") 
(urlays CEM a: 
+ 3 nolonOnty ny © Ey uylop , (11.8) 


Hi, H2 


We now divide both parts of (11.8) by C75 10190 and substitute 


1 
Lf 
nhl, i niyl> (11.9) 
C19 120 


LM 


——<- 


On the total angular-momentum representation the function pane has the 
standard asymptotic form 


I9 (tar) 

Tl e ”n > 
LM : brs aah nae LM 2 ee 
nlyly ™~ rR sin (Fn! _ 4 Dang Otyto + I nol 19nt yt ; (11.10) 

No ~ } 
where 
LM — LM asbaiis 
Fngl01ont sts C19 n0190 = \ Ci aie: Fig wOnt yp (1 1.1 1) 


Kel, $82 


Using the orthogonality of the Clebsch-Gordon coefficients 
2. Chin Cia 0) GMO AAO ig PO pant (11.12) 


Pigeylau2 “ly uy lo’ py Pola’ “yams “y2H2 
we can solve panto (11.11) to yield 


oO 0 
Lop2ly py oa LM LM LM 
A nol u8nlyuy = C1640190 Fol Int ta Crisis (11.13) 
i 
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Consider now the transformation of the matrix elements ™,}421242.,,. 
and W7y2!2 42, In accordance with the general formulae for the transfor- 


mation of matrices 
Ariat =) Ctinstana Citna'tatay Atta ay (11.14) 


Salona hahaa aah! 
Misa B2 
Commencing with ¥, we have, using the definition (10.20), 
1 1 
Via(1) a = Onn’ Ont! a (<7) (11.15) 
P Ir—r'| miyuan'ly ey’ 
We expand |r—r’|~* in terms of the Legendre polynomials 
lr—r'|o} = > y,(r, r’)P, (cos w) (11.16) 
A=0 
where w 1s the angle between r and r’ and 
nr? 
watt P 
re as (11.17) 
r 
r>r 
pAtt 


We then have 


I t t 
Via") are Onn’ Ont a \ Odanwn: | ye: ~')P (cos W) Yn (0 : ty) dQ 
A 
(11.18) 


where 
(Ya ntin'ty! = {Rural r'YRag (rr? dr’ (11.19) 
0 


(R,;, and R,-,,, are the radial wave functions for the atomic electron.) 
Substitution of (11.18) into (11.4) yields 


I 


Ta 5 15 5 | > re 
¥ nly pegn’'ly ey’ — = Onn’ Onn’ + 


5 y (Yantin'ty’ Yiua(O, (~) ) eg (p)P (cos a) tye (Os p') 
x Yi (0',@’')dQdQ’ (11.20) 


We now use the transformation (11.14), for simplicity denoting by v the 
group of indices 

v=nl,l, 
As a result we obtain 


1 
rLM _ + LM LM 
y yy 7 Oy + » (Yantyn'ty’ > Cr stous Cty’wa'la’na! 
A 


Mita’ 22’ 


x { ¥*, (0, 0)Yis-,,(0, @)P,(cosw) ¥*, (0', 0’) 
x Yio, (0',')dQdQ’ (11.21) 
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The expression 

h,= >. Cane Cin eae | Yea’, P)Y5*u,'(9, P)P (Cos w) GA Ce Q) 
—— x ¥, 4,0, g')dQdQ’ (11.22) 


which enters into (11.21) is independent of Af and can be expressed in 
terms of the Racah coefficient w: 


hay. 4, i, L) = (- 2 QA t+ 2 wh i, be, 1 A) 
x CP ro Chono( 2h, + 1)4(21, + 1)4(21, + 18215 +1)? (11.23) 


Therefore, 
it , , 
VA) = or 7 Ow’ + » (Yantai Aah, l,, lh, I), L) (11.24) 
ri 
For large r we have 
te oe 


so that (11.24) is a multipole expansion of ¥\. The quantity 4, has a 


non-zero value only under the following conditions: 
(1) The triads /,, /,, L and /;, /5, L satisfy the usual rules for the com- 
bination of angular momenta. 
(2) /,, 1, and /;, /; satisfy the condition expressing the conservation of 
parity, 1.e. 
(— jt = (ye th (11.25) 
(3) The parameter A satisfies the inequalities 
L-l)<4<d,4+] 
ie oaie aide als (11.26) 
Jy—-| <4 < (+4) 
which follow from the properties of Ci oie and Cr aio: when 2 = 0 


No = 611,'9 (11.27) 


Tal’ 


The quantity / is tabulated in [28] for certain values of /,, /,. /;, /3, and 4. 
We can transform 7 in a similar way. Starting with (10.52) and using 
(10.60) we obtain the following result 


We’, rh=t 2 Gay, 12,14, 12, L) 


X [Sr Eo— Fk) +yA Ra (Rey) (11.28) 
where 
Gli th 1, L) = (- I (20+ 12) wy, be, 1 4 AD 
x CPG 170 Chronyo(2h +) F(21, + IA(2H; + F215 + 1)# (11.29) 
and g, # 0 when 
Wy -i| <4 <(U, +1): 


The quantity g, 1s also tabulated in [28]. 


1-li| <A <(1,4+h) 
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It follows from (11.17) that when r < r’ 


t ks ie ’ 
WDE pr . r) = as > g, 50, (E, — >) + omi| Rat (vr Rn) (11.30) 
A 


/ 


whilst when r > r 
tat] 


WIMP) =D 93 | Bosley 4K) + | Ral Reig(t) (31) 
A 


On the total angular momentum representation, the integro-differential 
equations (11.6) for the partial waves assume the form 


d7fh (1, +1 r 
r+ | -AGFE| eazy ecto +f wr fe) in 
r r 7 ; 


(11.32) 
Henceforward we shall omit the superscript L whenever this leads to 
ambiguities. 


The possible values of /, for different 1,, /, and Z are given in Table 1 
for (—1)*" = (—1) and in Table 2 for (—1)""*" = (—1)™"*. 


TABLE 1 
m3 
L 1;0 2;0 22] 3; 0 a aa | 322 
0 0 0 l 0 | 2? 
l l I 0; 2 l 0:2 1; 3 
2 2 2 1; 3 2 L3 Q0;2;4 
3 3 3 2;4 3 2;4 Le 335 
TABLE 2 
md 
L pe | 3; 1 
0 
] l 1 
2 2 2 
3 1; 3 Ls 


The system of integro-differential equations (11.32) can be transformed 
into a system of integral equations: 


fo = Cryets( +O | Myer fer dr't+ D [Nyy fo(rdr’ (11.33) 
vy’ 0 v’ 0 
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where ¢,,.,, M@,, and N,,- are defined by 
M,y(r, r’) = 2fo (nur) —- v(r’ u(r) JY -(r') 


+2] [o(du(r")-2,0' 0] 
x Pwr”, AW Or", rr’) dr” (11.34) 


N,,-(r, rr’) = 2 | [our —o (ru (NW Or’, rv’) dr” (11.35) 
0 


Coe Sy t 2D | 0H re | Hole r’)v,(r") ir f(r’) dr’ (11.36) 
“0 0 


u,(r) as Jz Jina(ky, r) 


nk,r Z 
v(r) =-i y) Hi (k,, r) 


The integral equation for the partial amplitude is of the form 


(11.37) 


G2 22 > | seo ntrd + | U(r) Wr’, r’) i fr’) dr’ (11.38) 
vO 0 


We can carry out a transformation on (11.33), analogous to (2.39). In 
particular, we can introduce new functions x,,- which are related to f, by 


A) = y Coo Xyy() (1 1.39) 
We then have 


Xyyn = 0,4 2 [ Myx Mewar +E [Ns Xyye dr (11.40) 
v’ 0 
From equation (11.36) are obtained the following equations for c,,,: 


Cyov SO a C voy’ | Donen ] vr) Wr", r’) in| Nyye(r') dr’ 
: (11.41) 


On substituting (11.39) into (11.38), we obtain an expression for q,,, in terms 
of the functions x,,-: 


Yd vov ait 2 > Cyoy" oane +| Un Cle bane Ce r’) i Xyyelr’) dr’ 
= 0 (11.42) 
Let 
bee 8in=25 {| Wed fol rnde)de| sed 
2 (11.43) 


The matrix L,-, 1s a generalization of the Jost function given by (2.41). 
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We can use the properties of Bessel functions to write the matrix element 
Ly-, in the form 


LoS kp Ak, *Pik2 Bak? (11.44) 
where 
AiO ae eek key | (ua, Vt | tnd x ner 
v0 0 
(11.45) 


Ses a er ] , ae +| We yy ly, ir Recall (11.46) 
vO 0 

The quantities A,-, and B,-, are functions of k? and k%. Use can be made of 

(9.32), which represents the conservation of energy, to express all the k? 

in terms of kg. The diagonal elements are of a form analogous to (2.45): 


Ly, = Ay + ik??27'B,, (11.47) 
The system given by equation (11.41) may be written in the form 
yy C yoy” | = Oo (1 1.48) 
and the integral expression (11.42) in the form 
Yd vov Sa » Cyov" KE K2B yxy (11.49) 


In (11.48) the summation is carried out over all the states v”, whereas 
in (11.49) it is carried out only over those states which are consistent with 
the conservation of energy. If we neglect virtual transitions and retain in 
(11.48) only those states over which the summation is carried out in (11.49), 
the number of constants in (11.48) and (11.49) will be the same. 

If we consider all the possible values of vo, the quantities c,,, will form a 
matrix, which we denote by c. Equation (11.48) can then be written in 
the matrix form 

CL 4 (11.50) 


As with c,,,, the amplitudes q,,, depend on the subscript vp) and form the 
matrix q,,,. Let us introduce the matrix k whose elements are k,6,,. With 
its aid the expression (11.49) may be written in the form 

q = —ck'Bk' (11.51) 
Substituting c = L~', we have 
q = —L7'k’Bk’ (11.52) 
This is a generalization of (2.47) which expresses the scattering amplitude 
in terms of the Jost function. 

From the unitary relationships (10.60) we can obtain the corresponding 
relationship for the partial amplitudes. When the atom is in the ground 
state prior to collision we have 


Im abo—kolabol? = ¥ kulabsl (11.53) 
and/) = 0,/2 = L. 
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If the energy of the incident electron is too low for excitation to occur, 
the sum on the right-hand side will vanish and we return to (2.29). Let 


L _ fia! 


11.54 
doo ik, ( ) 
The partial cross-section for elastic scattering is then 
7 
Foo = (2L+ 1) 75 |Bx— 1? (11.55) 
0 


and the partial cross-section for all the inelastic processes 1s 
7 
» Cope Le Ia —|B,]7) (11.56) 
v#0 0 


which follows from (11.12). Hence the partial cross-section for all the 


inelastic processes is at all times smaller than 2(2L + 1)/K<. 


12. ELASTIC AND INELASTIC SCATTERING CROSS-SECTIONS. 
BEHAVIOUR NEAR THE THRESHOLD 


Let us now return to the general properties of the amplitudes. To begin 
with, we shall transform (11.52) to a more convenient form, using (11.44). 
In matrix notation, the latter expression is of the form 


L=k'Ak7'+ik'Bk'*? (12.1) 
Substituting 
A =—BM (12.2) 
into (12.1), we have 
L = —k'B(Mk~!—ik!*?) (12.3) 
Using (11.52) together with the latter expression, we obtain 
q = [k’B(Mk !—ik'*")]~*k'Bk! (12.4) 


or 
q = (Mk! —ik'*")7 'B™Pk7k'BE'! 
and hence finally 
q = ki(M—ik*!*4)7'k! (12.5) 
Bearing in mind (10.57), we obtain the following expression for the 
T-matrix 
(22)°T = k't4(M —ik*!* !)7 kt (12.6) 
This representation of the T-matrix is used in the general theory of multi- 
channel processes [29]. It is evident from the preceding discussion that this 
representation 1s a natural consequence of the properties of the Jost matrix. 
The formula (12.2) establishes the relation between the M-matrix and the 
Jost matrix. The matrix M is real and symmetric, which ensures that 


ToS oer (12.7) 
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We recall one further representation [30] also used in the theory of multi- 
channel processes involving the so-called reaction matrix K, which is given 
in terms of M by the expression 


M = k'*?K7!k'*2 (12.8) 
The amplitude q is 
q=k~?#K(1—iK)7'k7? (12.9) 
and the T-matrix 
(22)°T = K(1—iK)7! (12.10) 
It follows from (12.2) and (12.5) that the reaction matrix 
K = —k'*t?BA7 IK! +3 (12.11) 


The expression (12.5) 1s convenient for the investigation of the dependence 
of the cross-section on energy, particularly near the excitation threshold. 

As an example, consider a simple two-channel system [29] and suppose 
that the incident particle is in the s-state both before and after the collision. 
In spite of the fact that this is a much simplified example, it does exhibit 
certain general features characteristic of real multi-channel systems. 

By (12.5) the elements of the matrix q are of the form 


ey cn 12.12 
Qi dD ( ) 
M, 
= aos 12.13 
Gi2 = 421 D ( 3) 
M, 1 a ik, 
ae eee 12.14 
G22 dD ( ) 
where 
D =(M,,—ik,)(M,.—ik,)— M7, (12.15) 
k, = Vk?—-2AE (12.16) 
and AF 1s the excitation threshold. Below the threshold let 
ky = ix = i/2AE—K? (12.17) 
The effective cross-sections are then given by 
ky : 
Cin = 4 lial’ i,k = 1,2 (12.18) 


When the energy of the incident particle lies above the threshold, 1.e. 
when E, = 4kj > AE, the scattering and excitation cross-sections are 
respectively 

; M3, +k3 
— SSS ee ee ee 
(My1Mo.—Mi.—ky ka)? + (ky Maa thy My)? 
Ms Mi, 
ky (MyyMy2.—Myj2.—ky ka) +h, My2+k,M,,)° 


O11 (12.19) 


fp (12.20) 
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When the energy of the incident particle lies below the threshold, 1.e. 
when E, < AE, we have 


(M..+x)’ 
a,, =42n—-———-—-—-— (12.21) 
[M1 1(Ma.+K)—Mi2)7 +k{(M 2.44)’ 
which can also be written in the form 
4 
O41 — sin? 5 (12.22) 
ky 
where 
k?(M ¢)? 
sin? 6 = ma) (12.23) 


[M.1(M..+K)—M7{, |? +hk{(M2). +k)’ 
Consider the behaviour of the cross-sections near threshold. 


To within linear terms in &, in the amplitudes |q; ,| we have above the 
threshold 


Ci 5 = 4n|qh? |? (1 — 2K kg’) (12.24) 
k|q4o/ (0) (0){2 
1 


where the superscript 0 indicates that the particular quantity is taken where 
k, = 0, 1.e. at threshold. 
Below threshold 


O11 = 4a qh — 2k laa? x cot 6) (12.26) 
O11 
O74 
Ets = E:n 
(a) (b) 


Fic, 7. Elastic scattering cross-sections near the excitation threshold: (a) cot 6 < 0, 
(b) cot d > 0 (E,,, is the threshold energy) 


Comparison of equations (12.24) and (12.26) shows that do,,/dE, has 
a singularity at threshold.t At energies above threshold, ¢,, decreases with 
increasing £,, and the character of the variation depends on the sign of 
cot 6 below threshold. If cot 5 > 0, o,, decreases with decreasing £,. 
But, if cot 5 < 0, o,, increases with decreasing E). 

The variation of o,, with £, 1s illustrated in Fig. 7. Apart from the 
singularity at threshold, the cross-section ¢,, can also exhibit a resonance 


+ Elastic scattering at the threshold for inelastic scattering has been investigated by 
many authors (sec, for example, (31, 32)). 
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maximum in the sub-threshold region. As can be seen from equations 
(12.22) and (12.23), the maximum will occur when sin? 6 = 1, i.e. when 


M,,(M,,+%)—Mj,=0 (12.27) 


The width of the resonance peak depends on the ratio M,,/M,,. It can be 
shown [29] that the peak is sharp when |M, ,/M1,,| > 1. 

Figure 8 shows a typical variation of o with £ for the two kinds of 
threshold singularity and in the presence of resonance. 


(a) (b) 


Fic. 8. Elastic scattering resonance near the excitation threshold: (a) cotdéd < 0, 
(6) cotd > 0 


Numerical calculations carried out for a potential well [29] show that 
for some values of the parameters one does in fact obtain resonance peaks 
near the threshold. 

Equations (12.24)}-(12.26) are valid in the immediate neighbourhood 
of the threshold. It is possible to obtain approximate formulae which hold 
in a broader region by retaining in the amplitudes g;, terms which are 
quadratic in k,. With this intention, the quantities A, and B, in the Jost 
matrix (12.1) are expanded in powers of k?. An approximate formula 
for o,, has been obtained [33] in this way. A basically similar approach 
is used in [29], where the matrix MM; is expanded in powers of k3. Calcu- 
lating o,, by means of equation (12.21), in which the matrix elements M;, 
are represented by the two-term expansion 


M,,(«) = M,,(0)+ Ry, K? (12.28) 


it is possible to carry out an approximate study of the behaviour of o near 
resonance. This method is used in [34] to investigate the resonance in the 
elastic cross-section of the hydrogen atom for electron scattering, assuming 
a single excitation channel (excitation of the 2s-state). 

If the angular momentum of the particle after collision has a non-zero 
value and is equal to, say, /, then according to (12.5) the excitation amplitude 
9,, near the threshold varies as k5. The excitation cross-section near the 
threshold is 


O12 ~ ky"! (12.29) 
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The elastic scattering cross-section at the threshold will then no longer 
have a singularity of the form shown in Fig. 7a or 76 since do,,/dE 1s 
continuous at the threshold. 

The resonance in o,, below the threshold can also exist when / # 0. 
In the presence of three or more connected channels, the resonances can 
appear not only in the elastic scattering cross-section but also in the inelastic 
cross-section between two neighbouring thresholds. 

The case when there is degeneracy, 1.e. two or more channels at the 
same energy, occupies a special position. It is encountered in collisions 
with hydrogen atoms. Since the matrix elements V’,,, include elements 
which fall off as r~* as r > o, the behaviour of o near the threshold will 
be different from that indicated by equation (12.29) [35]. 


Chapter 4 


THE ELASTIC SCATTERING OF AN ELECTRON BY A 
HYDROGEN ATOM 


13. THE SCATTERING OF AN ELECTRON BY A HYDROGEN ATOM 
REGARDED AS A CENTRE OF FORCE 


As has already been pointed out, the elastic scattering of an electron by an 
atom can be described in various approximate ways as the superposition 
of three effects: the central-force interaction between the electron and the 
nucleus, the exchange interaction, and polarization of the atom by the 
electric field of the incident electron. All these effects are allowed for formally 
in the exact equations for the function Fp. If we take for F, the system of 
equations (10.55), the terms Vo9 Fy and [ Wo) Fo dr represent the central- 
force and exchange interactions between the electron and the atom in the 
ground state, while the sum 


: (VosFs + | WoaF dr) 


represents the effect of selnueaiba and the connection between the elastic 
and inelastic processes. By neglecting one of the factors it is possible to 
obtain simple equations which admit of a practical solution. Furthermore, 
the relative importance of the various factors may be estimated by comparing 
the results of different approximate calculations. 

The crudest approximation involves the rejection in the equation for Fo 
of all the terms containing F, and the term describing the exchange inter- 
action. In that case the system of equations reduces to a single equation. 
For simplicity, we shall omit the subscripts so that the equation for Fo 
becomes 


AF +(k*?—2V)F =0 (13.1) 
According to equation (10.20), 
Z 
Vir) ae 4 


Substituting the expression for ia wave enn of the hydrogen atom in 
the ground state 


Yo =e’ 
0 - 
and completing the integration, we have 


V(r) =—- m*(1 + 7 (13.2) 
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The asymptotic expression for F, 1s 
ikr 


Fw elt 4 a— (13.3) 


This approximation is equivalent to assuming that the atom behaves as a 
centre of force. 


On separating the angular variables 


fit G z ee : 2V | fi=0 (13.4) 


for which the systematic solution is _ 
id; | 
tim ~~ sin (1 - . ze i) 


Equation (13.4) has been solved numerically for /=0 and /=1 for 
0O< k? < 1-75 [36,37]. According to [37], the scattering length is 
a =—9-44, The large value of the scattering length and the negative sign 
indicate the presence of a virtual level near zero energy. 

A reasonable approximation to the numerical calculations [36] can be 
obtained by the method of trial potential (§ 6). We shall take the trial 
expression for Vo in the formT 


1 1 
SS r<R 
Vo = R r (13.5) 
0 r2>R 
where R is determined from the condition given by equation (6.7). 
For the s-wave at zero energy, we have 
1 1 
3-2(5- 7) fo=0 r<R 
: (13.6) 


fo =0 r>R 
When r < R, the solution given by (13.6), which satisfies the condition 


fo(0) = 0, is 
fo= cre F(t oe 2, 2/27) (13.7) 


where F is the degenerate hypergeometric function and C is a normalizing 
constant. When r 2 R the function f9 is of the form 


fo=rtq (13.8) 


t Expressions of this form were used in the 1930's to investigate certain general 
properties of the phase 6 [38]. In this paper the parameter R was determined with the 
aid of the semi-empirieal rules of Slater for the atomic radii (developed in the theory 
of atomic speetra). However, the phase 6 is very sensitive to changes in R, particularly 
at low energies. Therefore intuitive considerations leading merely to an order of magnitude 
for R are inadequate in the present context. 
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or on introducing the scattering length ay = —q, 
So =1—-ap (13.9) 
The continuity condition for the function and its derivative at r = R yields 
relations giving C and ap in terms of fo(R) and /6(R). In particular, 


The parameter R can be determined from equation (6.7). In the present 
case, this relation assumes the form 
R ¢) 


{(v- Vi) f2 dr +{ Vf2dr =0 (13.10) 


This equation has been solved numerically, and the result obtained for R is 

R= 1-43 (13.11) 
For this value of R, the magnitude of 2 [ |5V|r dr turns out to be equal to 
0-33. Therefore the condition given by (6.13) is satisfied and one can expect 


good agreement between wave functions calculated with V and V5. 
The results of the calculations are as follows. The scattering length is 


a =—10-3, quite close to the exact value of a = —9-44 [37]. 
Iq f 
12 fy 
te) 
8 
: O05 1 2 3 4 is) 


Fic. 9. The form of the functions fo and f for k = 0 


Figure 9 shows a plot of fo and a plot of the values of f found in [37] by 
numerical integration of the equation 


f"-2Vf=0 
In order to estimate the degree of approximation reached with the expres- 
sion for Vg with k # 0, consider the s-wave for kK # 0. The corresponding 
equation is of the form 
9 +(k? —2Vo) fo = 0 (13.12) 


6—2 
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The solution forr < R is 


1 
a cre™ F(1 a 240 (13.13) 
Ll 


2 
padeut (13.14) 


For r > R we can write the solution in the form 

9 = sin(kr+ do) (13.15) 
Both the constant c and 6, can be determined by matching equations (13.13) 
and (13.15) atr = R. We note that if gis an integer, the hypergeometric 


function becomes identical with the corresponding Laguerre polynomial 
and the calculation becomes particularly simple. For example, when 


where 


2 
hoe 1 = 0:40 (13.16) 


and jt = 1, we have 


fo=cre (13.17) 
A function f suitable for the comparison (k* = 0-45) is given in [36]. 


Figure 10 shows the values of fo and f. The phases 6, and o are, respectively, 
equal to 1.11 and 1.01. 


Fic. 10. The form of the functions fo and f for k? = 0-40 and k? = 0-45 respectively 


A considerable number of calculations have been carried out by the 
variational method. 


As in the analysis considered in § 7, one constructs the functional 
L, = [Wer rH Err) dry dry (13.18) 


where // 1s the encrgy operator for the entire electron + atom system. The 
trial function wy, should contain the atomic wave function which is not varied. 
The expression 


20 
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is Stationary with respect to variations in the trial function, and when the 
exact wave function yw 1s used the parameter Q is equal to the scattering 
amplitude. If the trial function contains mn parameters c, then by Kohn’s 
method these parameters are defined from the stationarity condition 


OL OL 
-= 0, ——=?2 13.20 
de aq = 
while the scattering amplitude is given by an expression analogous to (7.8): 
1 
Gk : a ie (13.21) 
TL 


If the atom is taken as equivalent to a centre of force, the trial function 
is of the form 


Wily, 72) = Vol DF A2) (13.22) 


where wo is the atomic function and F, is a function describing the scattered 
electron and containing the varied parameters. 
In order to obtain the partial scattering amplitude, g,, the function F, 
in (13.22) must be replaced by 
Ii) 


Fa = Yio(0,9) (13.23) 


For example, consider the variational calculation of the scattering 
length in [37]. Various forms of trial function were studied including the 
two-parameter function 


f,=r—a(l—e *)+bre" (13.24) 
The Kohn scattering length was found to be 
a = —9-38 (13.25) 


14, EXCHANGE EFFECTS 


A better approximation is obtained by retaining the exchange operator 
on the left-hand side of equation (10.55). This leads to the following 
integro-differential equation for Fo: 

AF# +(k?—2V)FE =2 { Wr’, n)Fe(r’) de’ (14.1) 
where V is given by (13.2) and 
‘ 1 1k? 
Wen) = (Eo 5-5) vol Mol? (14.2) 
jr’—r| 2 2 


If we separate the angular variables using the notation of (11.5) and omit 
both the superscript + and subscript 0, then 


I” +1 o 
4 a G = ” = 2¥, | f= 2 | Wr nA’) dr’ (14.3) 
0 
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where 
5 see. k? 1 
Wilr rn) =+4rr' e&—” lyr, r’) — ae (14.4) 
subject to the conditions 
f,(0) = 0 
‘ (14.5) 
fy~ evysin (19 acy + 3) 
The seattering length in the triplet and singlet states 1s [35] 
a, =8058, © a_ =2:347 (14.6) 


The wave functions are substantially modified by the introduction of the 
exchange interaction. For example, Fig. 11 shows the functions f*(r) and 


Fic. 11. The form of the functions f + and ffor k = 0 


f(r) for k = 0 in the absence of exchange. The phases 6g, 5* and 5} were 
calculated in [39] for 0 < A < 1. The results of the calculation are illus- 
trated in Fig. 12. Plots of 69, 6, and 6, which do not include exchange 
effects are given for comparison. It is evident that in many cascs the 
phase 6 lies between 5* and 67 in the absence of exchange. 

It is shown in [17] (§ 15) that in sufficiently accurate variational calcu- 
lations the phase 6 in the absence of exchange should tend to the half-sum 
4(6* + 57). 

The numerical solution of (14.3) must be carried out by the method of 
successive approximations owing to the presenee of the infinite integral 
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3 
(a) 
2 
50 
1 
Te) 
56 
00501 025 O5 LO 
ke 
0-4 5; 
(b) 
0-2 6, 
0 
oF 
-O:-2 
O25 os 1-O 


Fic. 12. The phase shifts 6.+ and dp (a), 6; + and 6, (6), and 52+ and 6 (c) 


on the right-hand side. It is shown in [39] that nine iterations are necessary. 
In this connection it is convenient to transform equation (14.3) to a form 
permitting direct solution without iteration. This transformation was 
carried out in [40] for the integral equation (11.35) which is equivalent to 
(14.3), and later by basically the same method for equation (14.3) [41]. 

We now consider a version of this involving the integral equation for 
the partial wave / = 0. Equation (11.33) assumes the form 


sin kr 


(aC | M(r,s)f(s) ds + | N(r, s)f(s) ds (14.7) 
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where 
r 


2 2 
M(r,s) = Fi sink(r—s)¥(s) + al sink(r—)N[WO(s, N—-W@X(s, NJ dt (14.8) 


AY 
r 


2 
N(r,s) = al sink(r—1t)h ts, Dd (14.9) 
0 
ene) i oe V(r) + [ W (r,s) ef | f(rydr (14.10) 
0 ~d 
To begin with, we consider the singlet case. Let f= cx, so that 
° k z 
= = a + | Mr s)x(s) ds +| N(r, s)x(s) ds (14.11) 
0 O 
] 
C= a (14.12) 
1-2 | G eth 4 | W (r,s) el is x(r) dr 
0 0 


Inspection of the expression for M(r, s) given by equation (14.9) shows 
that it is the product of the function of r: 


r 


8 
R(r) = = | sin k(r—tye'tdt (14.13) 
0 
and the function of s: 
ee 
O(s)=se *(-—k*—-1 (14.14) 
S 
Mey 
N = R(r)Q(s) (14.15) 
Bearing this in mind we shall rewrite equation (14.11) in the form 
sin kr 
x= P + | Mads + RO, (14.16) 
O 
Co = | xO ds (14.17) 
0 
Let 
X =X, +CyNXo (14.18) 


If we substitute (14.18) into (14.16) and require that x, satisfy the equation 


r 


Ls 
x, = ae kr + | Mr s)x,(s) ds (14.19) 


O 
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we have the following expression for x): 
2 
Xo = R(r) + -| M(r, s)xX(s) ds (14.20) 
0 


Equations (14.19) and (14.20) can be solved numerically without iteration. 
Therefore the problem formulated above is solved with the aid of trans- 
formation (14.18), which is possible because of the property indicated by 
(14.15). The unknown constant cg can be found in terms of x, and x,. 
Substituting equation (14.18) into (14.17), we obtain 


[ Ox, ds 
¢) = > (14.21) 
1—| Oxods 


The triplet case has a specific property associated with the fact that the 
antisymmetric function 


Wr 1,72) = fri Wolr2) —f(r2) Wor 1) (14.22) 
is invariant with respect to the transformation 
f(r,) 7 f(r) + const f(r) (14.23) 


It may be easily verified that equation (14.16) is also invariant with respect 
to (14.23) in the triplet case. This allows one of the two constants ¢ or Cg 
to be chosen arbitrarily. If, for one we set Cg = O then 


__ sin — 


‘| M(r, s)x(s) ds (14.24) 


We shall now show how to ieaiinate equation (14.3) to a form not requiring 


iteration. 
We have 


" I¢+1) 
eg G eae 


The right-hand side of this equation may be rewritten 


Es ay {22 i Wr, s)fds (14.25) 


i Wf ds = ( (WO wR) fds + { HOF ds (14.26) 
@) 0 0 


Since 2//*) may be written in the form 
20) = A(r)B(r’) (14.27) 


we have 


/ | WAOFds = A(ryc (14.28) 
O 
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where 
C= | Bfds (14.29) 
OQ 
Therefore, 
ldi+1 
I" + G — ( . — aa f =cA(r) (14.30) 
Substituting 
S=x+cy (14.31) 
we have two equations for x and y which can be solved without iteration: 
i+ 
x" + H? — ( _ — y x = | (#7 — wx ds (14,32) 
0 
(+1 : 
y" + G is 5 De y y=A(r) +{(WP—W My ds (14.33) 
0 


The unknown constant c can be found as before by substituting (14.31) 
into (14.29). The result is 


c =| Bx is[( ~ [By ts (14.34) 


In variational calculations, exchange effects are taken into account 
by using a symmetrized trial function y, in the form 


UW, = Fir, Wolr2) EF i(r2)Wolrs) (14.35) 


For example, we mention a variational calculation of the scattering length 
carried out in [37], where the trial function f(r) = rF(r) was chosen in the 
form given by (13.24). The scattering lengths were found to be 


a, = 8-11, OQ. = 2°38 (14.36) 


15. THE EFFECT OF POLARIZATION OF THE ATOM ON THE 
SCATTERING OF AN ELECTRON 


The effect of the polarization of an atom on the scattering of an electron 
may be investigated by including the term 


Y (VosFa+ | Woe Fe dr) 


in the equation for Fo [see (10.55)]. In spite of the fact that this term is the 
sum over all the excited states a, the polarization effect is not additive since 
in addition to the terms 


VioFo + | WaoFodr 


in the equation for F, there are also the terms 


>» (Va Fy +{ Woo F dr) 


a’ #2 


4, ELASTIC SCATTERING OF AN ELECTRON BY A HYDROGEN ATOM 79 


which represent the coupling between the different states. The effect is 
additive only when this coupling is neglected and the coupling 0 — «& is 
taken into account. 

We shall show that in the limit of large r the effect of excited states is 
equivalent to an additional interaction between the electron and the atom 
which can be characterized by a potential energy of the form U =—/r* 
[27], where f is the polarizability of the atom. We shall confine our attention 
to the approximation in which only the 0 > «& coupling is taken into account. 
Under these conditions the system (10.55) becomes 


AF, +(k2—2Vo,)Fo—2 { WooF dr’ =2 Y (VoaFe 2 { WF, dr’) (15.1) 


AF, +(k2—2V,,)F,—2 | W..F,dr' = 2(VeoFo i: { W.oFo dr’) (15.2) 


If the energy of the incident electron is such that excitation can take place, 
the right-hand side of (15.1) contains both the effect of polarization and the 
effect of coupling between excitation and scattering due to the unitary 
relation. But if the energy of the incident electron is insufficient for excita- 
tion, then the right-hand side of (15.1) contains only the polarization term. 
We shall begin with the latter case, i.e. we shall suppose that k2 < 0 for 
all « # 0 and only kj > 0. Consider equation (15.2) for large r. If k, < 0, 
the solution of the homogeneous equation AF,+k2F = 0 which is finite 
for r > oo will have the asymptotic form exp (—|k,|r). Under these con- 
ditions, the asymptotic behaviour of (15.2) is determined by the right-hand 
term. Since 


k? = k2+2(E,—E,) (15.3) 
we have 
AF, +(k3—2Va.)Fa—2 | Walt", r)F{r’) dr’ + 2(Eo—E,)F, 
= 2V.,F +2 { Wor nF o(r')dr’ (15.4) 


For large r the expansion of V,, in terms of the reciprocal powers of r 
begins with r~‘~!, where / is the angular momentum characterizing the state « 
(when / = 0, V,, falls off exponentially with increasing r). The kernel of 
the exchange operator W,, falls off exponentially as r + oo. Therefore for 
large r the right-hand side of equation (15.4) has the asymptotic form 


eikor 
yoind oa f q 
r 


The leading term on the left-hand side of equation (15.4) for large r is then 
(E,—E,)F,, so that we may assume that the asymptotic expression for F, 1s 


Viol o 
E,—Eo 


Pcs (15.5) 
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It is easy to show that for this asymptotic behaviour of F, the expression 
(A+k2—2V,,)F,—2 { WF dr’ 


in equation (15.4) may be neglected when r > od. In point of fact, according 
to equation (15.4), the term in F, which exhibits the slowest variation is 


r~'~! exp (ikgv.r), whilst in 


(A+k2—2V,,)F,-2 | W,, Fr’ 
the corresponding term is r~'~? 
(15.1) and let 


exp (ikyv.r). We substitute (15.5) into 


V,o(r)|? 
U(r) = -y | E.- e (15.6) 
WoCr wor 
W,(r',r) =— y ee (15.7) 


We then have 

AFo +(k2~-2Vo9 -2U) Fy = } (Woot W,)Fo dr’ (15.8) 
In this equation U enters as an addition to the potential energy due to 
excited states, while W, is an addition to the exchange operator. V9. JV 
and W, fall off exponentially as r > oo. As regards U, we need only con- 


sider the slowly varying part, 1.e. the dipole terms (/ = 1). For large r, 
we then have 


tw aa 5 
V0 2S ae aQ 


(the factor 2/3 appears as a result of the integration with respect to angles 
within the matrix element), so that 


Un~-—--— ot 15.9 
37 UE EF es 


2 
The quantity Poa 2 y ee ol is the polarizability of the hydrogen 
at#O0*aq  *0 


atom. Therefore, for large r, the function U(r) tends to —f/r*, which is 
the interaction energy between the electron and the induced dipole. Cal- 
culations [27] show that about 66% of f is due to the 2p-state. The discrete 
states contribute 81-4% in all, the remainder being due to the continuous 
spectrum. The numerical value of f 1s 9/2. 

Let us now return to the case where excitation is possible, i.e. one or 
more of the A? are positive. We can substitute 


: re Vio Fo 
=f. (15.10) 
E. =P, 
in equation (15.2) for the corresponding a. The leading term in the 


eikar 


asymptotic form of F, will be the diverging wave gq, a On substituting 
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(15.10) into (15.1), we have the following equation for Fj(r) 
AF, +(kz—2Vo9—2U) Fo —2 i (Woo + W,)Fo dr’ 
=2> (VoF, - i W,F, dr’) (15.11) 


a#0O 
When r —> oo, the right-hand side of this equation expresses the coupling 
between elastic and inelastic scattering channels. Similarly to equation 
(15.8), the left-hand side includes the polarization interaction U. 

The presence of the slowly decreasing potential (~ r~*) in (15.8) or 
(15.11) has an important effect on the asymptotic behaviour of the wave 
function. This can most easily be established in the limiting case of zero 
energy. For large r, the exchange terms fall off exponentially and are of 
little significance as far as the asymptotic expressions are concerned. We 
shall therefore consider the equation 


AF-—V(r)F =0 


The function F has the asymptotic form 


const 
Fw~i+——— 
For large r we shall set 
AF = V(r) 
and hence 
per- 2 [ Year 
An jr—r | 


Completing the integration with respect to the angles we obtain 


r ro 6) 


I 
F=1- “| Vr’? dr’ -| Vr' dr’ 
r 
0 r 
or 
I 2 ’ I 2 
F=1-——]| Vr'’dr’ —| Vr' dr’ +—- | Vr’* dr’ 
r r 
0 r r 
which amounts to the same thing. The function 
12 , 
J=|Vr'dr' —-| Vr'*dr (15.12) 
r 


determines the nature of the approach to the asymptotic form. 
If V falls off as exp (—4r) when r > oo, J will also fall off exponentially: 


- 
J=e “(1 +=) 
Ar 
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If, on the other hand, the polarization potential V ~ r~* is present, then 
1 
2r? 
In the absence of the polarization potential we have for large r 


Fr~l+ 


Cc t 
one + const e77” (15.13) 
} 


i.e. the approach to the asymptotic form is exponential. In the presence of 
the polarization potential, 


const const 
ee : rips (15.14) 


r 


and the asymptotic function is approached slowly in accordance with 
the r~7 law. 

It is possible to introduce polarization formally at finite distances, 
not only for r— co. However, it will then no longer be a function of 7 
as given by equation (15.6) but an integral operator [42, 43]. We shall 
show this for the case of the simplified system of equations in which the 
exchange terms have been neglected. Thus we have 


AF, +(k?—2V,,)F,, = 2Vio Fo (15.15) 
AF +(ki—2Vo9) Fo = 2 \ Vogl a (15.16) 


We shall confine our attention to the case k? < 0. Suppose that G, is the 
resolvent defined by 


G, = (A+k?+i0—2V,,)7! (15.17) 
The solution of (15.15) is then 
F, = 2G6,VoF 5 (15.18) 
Substituting this into (15.16) and introducing the operator 
U = » Vor. Vio (15.19) 
we have 
AF) +(k2 —2Vo9—2U0)F 5 = (15.20) 


Let us now consider the numerical calculations which include polariza- 
tion effects. Calculations of the phase shifts for L = 0, 1, 2 with allowance 
for the 1s—2s coupling in the interval 0 < k* < 0-75 are given in [44] 
(A? = 0-75 corresponds to the threshold for the excitation of the 2s-level). 
Calculations in which the ls—2s—2p coupling was taken into account 
are reported in [45]. When the ls—2s coupling is allowed for, it is found 
that 

ay = 75, A/S 92733 
while the result for 1s—2s—2p coupling is 


a, = 6:74, a_ = 1-89 
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The phase shifts 6g and 69 are plotted in Fig. 13. Calculations without 
polarization and with polarization allowed for by the variational method 
(see below) are given for comparison. 


02 04 O06 08 O O02 04 O06 O8 


k2 k* 


(a) (b) 


Fic. 13. The phase shifts d9* (a) and d9~ (6): 1—without allowance for polarization, 
2—with ls — 2s coupling, 3—with ls — 2s ~ 2p coupling, 4—-variational method with 
polarization 


The relative importance of the various couplings may be deduced from 
Table 3 which gives 63 for two values of k*. As can be seen, Ils—2s and 
ls—2p coupling give rise to roughly the same contributions. The simul- 
taneous inclusion of 2s and 2p states leads to a much greater contribution 
than the sum of the 1s—2s and ls—2p contributions. Therefore, the 2s—2p 
coupling is significant and the effect of polarization is not additive. 


TABLE 3 


+ ra 
States included 35 for k? equal to 


0:55 0-60 
Is 0-7004 0-6704 
ls — 2s 0-7352 0:7115 
ls — 2p 0:7338 0-7040 
Is — 2s — 2p 0-7846 Q:7707 
ls — 2s — 2p — 3s 0:7894 0:7738 
ls — 2s — 2p — 3p 0-7920 0-7738 
ls — 2s — 2p — 3d 0:7841 0-:7694 
ls — 2s — 2p —3s — 3p 0-7975 0:7814 
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Another feature evident from Table 3 1s that the coupling with the 3s, 
3p and 3d states adds very little to the contribution due to the ls—2s—2p 
coupling. On the other hand, it cannot be considered that the inclusion of 
only the 2s and 2p states is sufficient since the contribution to static polariza- 
bility of the hydrogen atom due to the 2p state is 667%. It should therefore 
be assumed that the convergence is very slow and a large number of excited 
states must be taken into account to improve the approximation which is 
achieved with the Is—2s—2p coupling only. 


O O02 04 06 08 !10 12 14 
ke 


Fic. 14. Peak in the elastic scattering cross-section [44] 


As can be seen from Fig. 13, the phase shift increases very rapidly near 
the threshold for the excitation of the 2s-level. The variation in the phase 
shift in the region of rapid increase can be described by the Breit-Wigner 
formula with resonance energy E = 9-61 eV (k* = 0-711) and a full width 
at half height [ = 0:109 eV [45]. As can be seen, this is a very narrow 
resonancc whose presence is predicted by the general theory ($12). In 
addition, there is also a threshold peak. Figure 14 shows the scattering 
Cross-section as a function of energy with a clearly defined pcak [44]. 
Analogous results, ic. the narrow resonancc, were reportcd in (34). 
It should be noted that detailed analysis of the variation in the 
phase near the 2s excitation threshold was madc difficult in (44, 45] 
because the convergence of the iteration method used to solve the 
integro-differential equations deteriorated as threshold was approached. 
It is evident from Fig. 13 that the phase shift 69 is not modified to any great 
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extent when the Is—2s—2p coupling is taken into account. Figure 15 
shows k cot 6* as a function of k*. As can be seen, k cot 59 is appreciably 
non-linear for small A*. This is due to the terms & and k? Ink in equation 
(3.24). In the case of k cot 59 the departure from linearity is very slight. 


> 
0 


k cotd 


0:05 O10 0:15 
42 
(a) (b) 


Fic. 15. Dependence of k cot d9* (a) and k cot dg~ (6) on k?: 1—with ls — 2s — 2p 
coupling, 2—without polarization 


This is due to the fact that the coefficient of & 1s inversely proportional to 
the square of the scattering length and is therefore much smaller in the 
singlet state than in the triplet state. 

Figure 16 shows the phase shifts df and 63.+ Here again one can see 
a rapid increase in the phase near the threshold for the excitation of the 


(b) 
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Fic. 16. The phase shifts 6,*+ (a), 6:7 (6), d2* (c) and 627 (d): 1—without polarization, 
2—with ls — 2s coupling, 3—with ls — 2s — 2p coupling, 4—by the method of polarized 
orbitals 


+ For comparison, the results of calculations including the polarization effects and 
using a different method are also given (see below). 
E.A.C. 7 
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2s-level. As has been pointed out in [45], the phase shifts 53 and 6; approach 
each other at low energies and are adequately represented by equation (3.25) 
if B is replaced by 66% of the true polarizability of the hydrogen atom, 
i.e. by the contribution due to the 2p state. Thus, when k* = 0-1 the values 
of 5; calculated in [45] are 63 = 0-0084 and 6; = 0-:0096, whereas equa- 
tion (3.25) yields 6, = 0:0089. 

There is another way of including the polarization effects. In this method 
the basic functions for the expansion of the wave function of the atom + 
electron system are taken to be the wave functions of the atomic electron in 
the field of two centres of force, i.e. the nucleus and the incident electron, 
rather than the atomic wave function as in (10.17) or (10.46) [46]. This 
method is a development of the so-called method of perturbed stationary 
states described in the book of Mott and Massey. In [46] it is referred to as 
the method of polarized orbitals. 

We first consider the elastic scattering in the absence of exchange. We 
shall write the wave function for the atom + electron system in the form 


W(ri,%2) = F(ry)P(r,, 2) (15.21) 


where F(r,) describes the incident electron and M(r,, r,) describes the atomic 
electron in the field of the nucleus and the incident electron. The function 
@(r,, r,) can in its turn be written as the sum 


P(r, 12) = Wolr2) +901 ¥2) (15.22) 
where y,(r,) 18 the wave function of the atomic electron in the isolated atom 
and g(r,,1r,) is the change in the atomic function at r, duc to the incident 
electron at r,. The equation for F can be obtained from equation (10.24): 


[ Vo(ra(H—EW(r1.r2)dr2 = 0 (15.23) 
where & = 1k*+E, and the cnergy of the ground state of the atom is 
Eo = —}, 

We substitute the expression for F/ given by (9.1) and the expression for 


W(r,, 72) given by (15.22) into (15.23). We shall, moreover, require that 
g(r,,"2) and w,(r,) should be orthogonal for allr,: 


[oc r2 )Wolr2) dr, =0 
We then have 
AF +(k?-2V—-2V,)F =0 (15.24) 
where V is the potential cnergy of the clectron in the field of the undisturbed 
atom given by (13.2) and V, is the polarization correction to the potential 
energy and is given by 
Wolr saris) , 


rir] 


V,(r;) = r, (15.25) 


The problem is therefore reduced to the determination of the function 
Jr 1, V2). 
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The equation for g(r,, r.) can be obtained by substituting (15.21) and 
(15.22) into (H— &)W(r,, r.) = 0 and using (15.24) and the equation satisfied 
by wW,(r,). In the present method, the determination of g(r,, r.) is carried 
out on the adiabatic approximation in which the position of the incident 
electron r, 1s regarded as fixed and all the derivatives of g(r,, r.) with respect 
to r; are equated to zero. Moreover, if itis assumed that g is a small correction 
to Wo, the determination of g(r,;,r,) can be simplified still further and g 
can be neglected in all the terms into which the sum W)+g enters. We then 
have the equation for g 


I 1 
(142 +4) oro) = | — Mees) | vols) (15.26) 


An accurate solution of this equation can only be found by numerical 
integration. However, the asymptotic behaviour of g for r; > r, can be 
found in analytical form. Whenr, — oo the potential V(r,) falls off expo- 
nentially and we can neglect V(r,) on the right-hand side of (15.26). Next, 
if we expand |r, —r,|~* in terms of the spherical harmonics and retain the 
dipole term only, we have 


1 
—— —5 COS 0 (15.27) 


where @ is the angle between r, and r,. When these simplifications are 
introduced, equation (15.26) 1s replaced by the asymptotic equation 


1 r 
(442 += = 4) g0lrisra) = 273 voles) cosd (15.28) 
2 1 
where gy is the asymptotic function. The solution of this equation is 
cos@ _| 
Jo(r1,'2) = 5 Te (ry +4713) (15.29) 
rea) n 


This solution can be augmented by the solution of the homogeneous equation 
obtained from (15.28) by rejecting the right-hand side. However, since g 
and Wo are orthogonal, the additional function is equal to zero. 

As regards the approximate solution of (15.26), we must first mention 
[47] which gives a variational method of solving this equation. The method 
is based on the fact that the function 


t= {ater |[(d2 += = 1) ore) = vole] dra (15.30) 
r pet 
is a minimum with respect to an arbitrary variation of g (for any r,), and 
its stationary value is equal to the polarization potential V, given by equation 
(15.25). Another approximation is described in [46], where the polarization 
potential is calculated by replacing g by its asymptotic form gp for r, > rz 
es 
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and assuming g = Oforr, < r,: 
g = &r1,62)9o 
1 eee a (15.31) 
E= 
0 PT 
Having determined g(r,,r,) by any of these approximate methods we can 
calculate the polarization correction to the electron—atom interaction energy 
from equation (15.25). 


In order to take the exchange effect into account we must construct the 
symmetrized function 


Wry.12) = F(ry) [Wol(r2) + 9(11 ryt Fe alWolri)+9(r2,71)] (15.32) 
where g satisfies (15.26) as before. Substituting equation (15.32) into (15.23), 
we obtain an equation of the form 


AF +(k?-2V —2V,)F = 2 | (W+W,)F dr’ (15.33) 


where W, 1s the polarization correction to the exchange operator. As before, 
the angular variables of this equation can now be separated and equations 
with radial functions obtained. 
The equations derived from the approximate expression (15.31) for g 
are [46]: 
fo 6) 


i? (i+1 
me am G -- r? 27 | fi =| murs ar’ 


9 
=2 5 ht NOhEPD (15.34) 


where ¥ and //’,, are given by equations (13.2) and (14.4). 

The right-hand side contains terms describing the polarization effect. 
The first term contains 9/2r*, i.e. it is the asymptotic form of the polarization 
potential with which we are already familiar. The remaining terms are: 


N =— 20° 7 (r? +374 4-9r? 4429? 4227422 (15.35) 
d 
PU) = 3 +k)5y [owe ago 
l ee ren f(r) dr 
a ane ee i-1 To OIE 
- =— pre 


(i+1)r'*! e~"fi(r)dr 
(21 4-1)(214+3) pore 


r 


| e~'(r+4r?) 


+ 4r7 eo" diy rt 


r 


me r (15.36) 
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The equations (15.34) are transformed to a form not requiring iteration 
and solved numerically for / = 0, 1,2 and & between 0 and 1. The wave 
functions are calculated up to r = 20 after which they are assumed to be 
identical with the asymptotic form 


In 
sin (11 ers + 7 


Figure 16 shows plots of df and 63. The scattering lengths are found to be 
a” = 5:8,a~ = 1-9, whereas the corresponding values without the polariza- 
tion corrections are respectively 8-1 and 2:35. However, subsequent calcu- 
lations [48] show that even for r > 20 there 1s an appreciable contribution 
due to the polarization potential. The change in the phase shift 6, between 
r= 20 and r = © can most easily be determined from equation (3.35) 
which for / = 0 is of the form 


dy _ 
dr 


For r > 20 all the exchange terms in (15.34) are practically zero and 
this equation is valid. At such large distances we need only retain the 
asymptotic part f/r* in the potential V. If, in sin? (kr+65,), we replace 
Oo(r) by its limiting value 6,(00) which is valid for small changes in the 
phase shift, we have 


2V 
— = sin’ (kr +59) (15.37) 


d9()—d0(R) = 


fe 6) 9 
k 
2B | OA 15.38) 
k r 
R 


Hence, to within terms of the order of k, 


| a? a? 1 Tt 

z Lo0( 20) — oR) = Ge — Pos + p,) + a + o(k? Ink) (15.39) 
where a is the exact scattering length. This formula was derived in another 
way in [48]. 

We note that equation (15.39) contains the linear term Azf}/3 which is 
identical with the linear term in equation (3.24). When & = 0 we obtain the 
relation between the accurate scattering length a and the scattering length 
a(R) which is obtained from the asymptotic form for R = 20: 


fat a 1 
a=a(R)-P (555-7 +R] (15.40) 
Calculations yield a@~ = 1:74. The correction to the singlet scattering 


length a* is less significant because a” is small. 

In variational calculations, polarization is taken into account by including 
in the trial function terms containing r,;, = [r,—r,|, so that W(r,, r,) is not 
reduced to the symmetrized combination of the products of functions of r, 
and functions of r, of the form given by equation (14.35). 

Determination of the scattering length is particularly important in 
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variational calculations. In the case of scattering of an electron by an atom, 
it turns out, just as in the simpler case of scattering by a central force (§ 7), 
that the scattering length calculated by Kohn’s method is the upper bound 
on the true scattering length [49]. This provides us with a criterion for 
estimating the accuracy of the variational calculations. In the triplet state 
we have a@~ < ag where dg is Kohn’s scattering length: 


az = ar +2 i Wy (H —E,)W, dt (15.41) 


H is the energy operator for the electron + atom system, Ep 1s the ground 
state energy of the hydrogen atom, yw, is the trial function whose asymptotic 
form 1s 


Wo = (1—P,)Wolr) (2: 7 . (15.42) 


and Wo is the atomic wave function. 


In the singlet state, we have the following inequality which is analogous 
to (7.29): 


(15.43) 


<a, 42 {yt (H—E,)W; dt — a || wu E,W; ea 


where W’ and E” are the eigenfunction and energy eigenvalue of the bound 
state of the electron + atom system (it is assumed that there 1s only one 
such state). The function w’ must be a sufficiently good approximation to 
the eigenfunction for the negative hydrogen ion in order to ensure that 


E' =| W'Hy' dt < Eg (15.44) 
The inequality given by (15.43) can also be written in the form 
at <a; (15.45) 
provided that one uses the function & varied by Kohn’s method and including 
the bound-state function w’, so that 
Wi, = y+ Ay’ (15.46) 
where A 1s the varied parameter. 

It is shown in [49] that even if the function yw does not have the structure 
indicated by equation (15.46), the inequality given by (15.45) 1s still obeyed 
under certain conditions. 

We now briefly review the numerical calculations. The most systematic 
are given in [50]. The s-wave phase shift was calculated from the wave 
function w for the electron + atom system of the form 

y= pty 
where 


= (LEP 2207 mo flr) (15.47) 
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and 
1 k _X 
(ej Si ee) (15.48) 
k, k 
_ ttn 
X= bY Comne 2 riyo(rirgtrir3) (15.49) 
I,m,n>0 


| ; : 
The function — tan 6 was defined as the stationary value of the functional 


k 


1 

itgd =A - | W(H —E)w dr, dr, (15.50) 
The set of numbers /, m,n in equation (15.49) was chosen so that the 

sum /-+2+7 varied between 2 and 7 in the triplet state and between | and 6 


in the singlet state. Correspondingly, the coefficients C;,,, ,, formed matrices 
of rank 3, 7, 13, 22, 34 and SO. 


, I 
Figure 17 shows the stationary state values of — ; tan O, for the triplet 


s-wave with k = 0-4 (about 2 eV) as a function of the scale parameter kx. 
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Fic. 17. Dependence of — (1/k) tan d9~ on the scale parameter x for k = 0-4; the 
numbers indicated on the curves correspond to matrices of the following rank: 34 (1), 
15 (2), 7 (3) and 3 (4) 


When xk = 1-55, the curve corresponding to the matrix of rank 7 exhibits 
a discontinuity. Similar properties are exhibited by all matrices of higher 
rank but they are so narrow that they have not been resolved. 

When the number of parameters is large, the curves become smoother 
and approach each other. If it is assumed that as the number of parameters 
is increased all the curves tend to a single common curve, then it is possible 
to determine the corresponding phase shifts by extrapolation. Figure 12 
shows the extrapolated phase shifts for the triplet and singlet states. 
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For the calculation of the scattering length the function m was chosen 
in the form 


225 ao 


p= (1epiad2e" [1-2 4, - fs -E fyeoso] (15.51) 
1 


2 


where 


dp _.. 
nn | Ete ee (n = 1, 2,3) 


The dependence of ~ on r, was chosen so as to reflect correctly the effect 
of the polarization of the atom on the asymptotic form of the wave function. 
Figure 18 shows the scattering length as a function of the parameter x. 


0:8 16 08 16 
xX xX 
(a) (b) 


Fic. 18. Scattering lengths a, and a_ as functions of the parameter x. 


Table 4 gives a summary of the scattering lengths calculated by the various 
methods. 


TABLE 4 


ls — 2s —2p Method of polarized Variational 


coupling orbitals method 
a, 6:74 5°8 5:96] 
Qe 1-89 1-74 1-768 


16. DISPERSION RELATIONS FOR THE SCATTERING AMPLITUDE IN 
THE CASE OF THE HYDROGEN ATOM 


The dispersion relations for the scattering amplitude in the case of the 
scattering of an electron by the hydrogen atom can be deduced by means of 
an analysis analogous to § 4 [5], 52]. As in $4 the dispersion relations 
are based on the properties of the Green’s function G(E+£,)+i0) for the 
atom + electron system, where £ is the energy of the incident clectron, 
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and Eo is the energy of the hydrogen atom in the ground state. We shall 
assume, by analogy with the scattering of a particle by a centre of force, 
that G(E+£,+i0) is an analytic function in the plane of the complex 
variable E with a cut along the positive real axis, and has simple poles on 
the negative real axis at the points E = —|£,| which correspond to the bound 
states of the negative hydrogen ion H~. The residues at these points are 
equal to wi(r,, PU (ry, 72) where w,(r,,r2) is the eigenfunction of H™. 
Moreover, there are branch points on the real positive axis which correspond 
to thresholds for the excitation of the atomic levels. 

Using the expressions given by (10.14) and (10.16) for the direct and 
exchange amplitudes, and substituting the expression for yw through the 
resolvent G into these expressions, we have from (10.1) 


W(ri.r2) = ene Wolr1)— G(E+Eg+ i0)[Vo(r2)+ V(r, r)] one Wolrs) 
(16.1) 


Since we are interested in the forward scattering amplitude, we shall set 
v =n. This yields 


fa fors-| eM ravers) [Votes + Cr. 19) 

x G(E + E, + i0) [Vo(r5) + V(r), r5)] et 2Wo(r}) dri dr, (16.2) 
and 
9 = 90+ | eM WSS LVOCr)+ MH] 

x G(E+ Ey +10) [Vo(r) +V(r4, 75) ] 2 Wo(ri) dr dr, (16.3) 


where 


fo=— 5- | Walrad + Verio rll volri)P dr ds (16.4) 


I t t ikv.(r2’—r,’ t ’ / t 
Jo = — 5 | [Voli + Verio] Vole WSC) dri drs (16.5) 


The functions fj and gy are real and can readily be given in an explicit form. 
For f—fo and g—gy we have expressions analogous to (4.16): 


Imf(E’) d 


= J 16.6 
E-E' L E+E oo) 


l 
Ref =f += 9 | dE 
O 


and 


1 Im g(E’) C; 
= —Ff | dE’ —— - J 16.7 
pee aor |. E—E' LENE one 
e) 
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where 
d,;= | eI" (ro) [Vo(ri) +V (ry, ro) Jury, r2) dr; dr. 

x [mmole [MOF +V ri rDuH(ri, rDddr{ dry (16.8) 
Cc; = | dr, dr,e ""? Wo(ry)[Vo(r2) + V(ry, r2)]uj(ry, r2) 


x | ary dr es o(rs) [Vole + VC ru fri.) 16.) 


and F represents the prineipal value of the integral, while x; = J 2|E;|. 
The singlet ground state of the negative hydrogen atom (binding energy 

~ (0-7 eV) is well known at the present time. There are grounds for supposing 

that H~ does not have exeited states. We shall continue our analysis on the 

assumption that there is only one singlet state. In this state the funetion 

u,(r,, 2) is symmetric and therefore we have from (16.9) and (16.8) 

2 


et Wo(ry) [Volri)+ Vr, 62) Ju(ri,r2) dr, dr, (16.10) 


Substituting the symmetrie and antisymmetric eombinations of f and g 
and using (10.45) we have 


k'o ©) ap 2c 
Req* = qo +73 39 | ~ E+E, (16.11) 
1 


1 
SF 
: 2m 


poe ) 


Req” =40 +73 ne fa (16.12) 


As was to be expected, the singlet eee state of the negative hydrogen ion 
has an effeet on the singlet scattering amplitude gq” but has no effeet on the 
triplet amplitude g°. For eomparison with experiment one must use a 
formula for the eross-seetion <a> whieh is equal to 307 +407. It is readily 
shown that 


k 
7, <9 = Im(f-49) (16.13) 
TU 
The dispersion relation for f—4g 1s 


k! ct c 
—E 2(E+|E,|) 


Re(f—49) = (f-4g)o + 3 ai dE’ (16.14) 


At zero energy (16.14) yields the foitewins expression for the seattering 
lengths: 


ae eer mee ee sa esieca (16.15) 


a 


4. ELASTIC SCATTERING OF AN ELECTRON BY A HYDROGEN ATOM 95 


An interesting example of the application of (16.15) is the analysis of the 
two contradictory experiments on the scattering of electrons by hydrogen 
atoms [53, 54] which is given in [52]. The analysis is based on the comparison 


of the values of | (o> dk from (16.15) and the same quantity determined 
0 


from experimental data. 
First we note that because a* and a™ are positive we have immediately 
the inequality 


= c 
[ <o> dk < 2n? aE do. $40] (16.16) 


0 


If we now calculate fo and gy from (16.4) and (16.5) we have 
fo-490=-2 (16.17) 


The quantity c/2/E,| cannot be calculated exactly since the wave function 
for the ground state of the H™ ion is not known. It is pointed out in [51] 
that the approximate values available for this quantity vary by about 20%, 
depending upon which approximation is used for the wave function of H™. 
Using calculated values for the quantities on the right-hand side of (16.16) 
we have 


{ Ya) dk < 7-52 (16.18) 


On substituting the values of a* and a~ calculated in (46, 48] into (16.15) 
we obtain 


| {a> dk = 4-77 (16.19) 


oO 
Let us now return to the determination of i <a> dk with the aid of 


0 
experimental data. This was done in [51] by the following procedure. At 
low energies (S 10 eV) use was made of the experimental data on elastic 
scattering of clectrons by hydrogen atoms. At high energies (2 500 eV) 
the determination was based on the values calculated on the Born approxi- 
mation. The cross-section <a> in the intermediate region was dctermined 
oO 


by interpolation. The uncertainty in | <a») dk due to interpolation errors 


0 
can be estimated quite easily, and the final result is 


i 
Ee Yo) dk = 10:340-5 (16.20) 


0 
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when we use the results of [53], and 
| 
xa | <> dk = 5:1+0-°5 (16.21) 
TL 
0 


when we use the results of [54]. 
As can be seen, (16.20) lies well above (16.18), whereas (16.21) is in 
agreement with (16.18) and (16.19). 


17. ANGULAR DISTRIBUTION OF ELECTRONS AND THE SCATTERING 
CROSS-SECTION. COMPARISON WITH EXPERIMENT 
In this section we shall considcr the differcntial scattcring cross-sections 
averaged over the spins and calculated from the phase shifts 6g, 57 and 
6; with allowance for the 1s—2s—2p coupling for four values of the 
energy [45]. The corresponding curves are shown in Fig. 19. The points 


in this figure represent the experimental data on the angular distribution 
as reported in [54, 55]. 
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One notes the characteristic form of the angular distribution of the 
scattered clectrons. At low energics (~ 3 eV) forward and backward 
scattering predominate. This effect is duc to the cxchange interaction. 
As a result of exchange the phase shifts 5{ and 5} are much greater than 
the values 6, and 6, which are obtained without allowing for the exchange. 
Figure 20 shows the total clastic scattering cross-scction (or morc preciscly 
the sum of the three first partial cross-sections). The various experimental 
points are given for comparison [45]. 


Chapter 5 


EXCITATION OF THE HYDROGEN ATOM BY 
ELECTRON IMPACT 


18. APPROXIMATE METHOD FOR THE CALCULATION OF THE 
EFFECTIVE EXCITATION CROSS-SECTION 


THE various approximate methods can conveniently be discussed on the 
basis of the system of equations given by (10.55). If we confine our attention 
to direct transitions between the initial and final states and neglect virtual 
transitions to intermediate states and also the coupling between different 
real transitions in the case of multichannel processes, then the system given 
by (10.55) reduces to the following two equations: 


AF ,+(k2 —2V_a)Fa—2 | Wag Fa dr’ = 2V a9 Fo+2 | Woo Fo dr’ (18.1) 


AF y+(k2—2V9)Fo —2 i WyFodr' = 2(Voa Fat WF, dr’) (18.2) 
The quantities V,,, Wia. Voo, Woo and W,, fall off exponentially as r > oo. 
The quantity V,, falls off as r~'~’ when r — oo (except in the case of tran- 
sitions between s-states when V,, is also found to be exponential). If we 
retain only the slowly decreasing interactions represented by V,. and neglect 
all others we obtain 


AF, +k?F, = 2V,0F 5 (18.3) 
AF o+koFo = 2Vo.F (18.4) 


The simplest and crudest approximation is the Born approximation in 
which V,, is regarded as a small perturbation. The term Vp, F, on the right- 
hand side of (18.4) is then neglected and a plane wave is substituted for Fo 
on the right-hand side of (18.3). As a result we havet 


AF, +k2F, = 2Vi9 el” (18.5) 
The amplitude g, on the Born approximation 1s 
1 
Goo = -— oe tkan — kov)r V0 dr’ (18.6) 
20 


+ There exists also the so-called Born-Oppenheimer approximation in which (18.5) 
is augmented by the exchange term { Wao e'*o"-r’ dr’ where the function Fo is replaced by 
a plane wave. However, numerical calculations have shown that this leads to a con- 
siderable overestimate of the cross-section. In many cases the cross-section calculated on 
this approximation exceeds the theoretical Jimit by a large factor. A modification which 
removes much of the discrepancy is described in [105). 
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Let us introduce the vector 
K=k,n—kopv (18.7) 


which represents the momentum transfer. In this notation we have 


1 e Kr Vv dr (18.8) 


On 


We can now use the expression for V,,, given by (10.20) to transform (18.8) 
to the form (cf (5.25)) 


—iK.r 2 a 
= aa] eo Wo dr = Ge ac (18.9) 


The total excitation cross-section is given by 


doa 


Co. = an | |q.| sin 0 dO 
0 
O 


The integration with respect to @ may be replaced by integration with 
respect to K. From (18.2) we have 


K? = k?+k2—2k,k,cos0 (18.10) 
Since 
k? = kg—-2AE (18.11) 
where 
KES ESE. (18.12) 
it follows that 
K?* = 2(AE—k,k, cos 0) (18.13) 
and hence 
KdK = kyk,sin@ dé (18.14) 
Therefore, 
kotk, : 
oun ge | Keel (18.15) 
ko—ke 


When the momentum transfer is small, the matrix element (e'*-"),, can be 
expanded into a series in powers of A. The first term will vanish because 
the wave functions of the excited and ground states of the atom are 
orthogonal. For optically allowed transitions (e'*-"),, ~ K(r),9, which 
may be referred to as the dipole approximation. The quantity (r),) in its 
turn is related to the oscillator strength for the transition 0 > «& [56]. In 
general, the dipole approximation for (e'*-”),, cannot be used in the cntire 
range of integration in (18.15) since at the uppcr limit the quantity A may 
not, in fact, be small [57]. However, the cxprcession given by (18.15) may be 
modified by introducing instead of the uppcr limit 4 y+,, an cffective 
number k’ which is chosen so that the dipole approximation can be used in 
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the entire range of integration. The number k’ can, for example, be deter- 
mined from the condition that for sufficiently high energies the magnitude 
of the cross-section o calculated from (18.15) (without expansion in the 
powers of K) should be practically identical with the value of o calculated 
on the dipole approximation [58]. On substituting the dipole approximation 
for (e'*-")_, into (18.15) we have 


4 


Ko—k, 
Since, on the dipole approximation, the cross-section o is a logarithmic 
function of k’, the latter need not be determined to a high degree of accuracy. 
On the dipole approximation, therefore, the excitation cross-section do, 
can be expressed in terms of the oscillator strength. 

On the Born approximation the expression for the partial amplitude 
is obtained with the aid of the partial-wave expansion (18.8). On the total 
angular momentum representation we then have 


(18.16) 


, oT 
Tox = 53 (Ma0ln 


dw =—2 | ko wit; (kyr) dr (18.17) 
0 


We note that even at low energies (near the threshold), the amplitude 
(18.17) for a given transition is not a monotonically decreasing function of 
the angular momentum /, of the incident electron. In contrast to elastic 
scattering, when the s-wave has a dominating effect at low energies while 
the contribution due to partial waves with / # 0 1s small and decreases 
rapidly with /, it turns out that in the present case the dominant partial 
wave is that corresponding to a value of /, which is not the smallest value. 
This phenomenon was established, for example, in (56] as a result of calcu- 
lations of the excitation of p-levels in He. The reason why a partial wave 
corresponding to an angular momentum /, may give a bigger contribution 
than waves with /,—1 and /,+1 can be deduced by analysing the properties 
of the integrand in (18.17). For example, when the atom undergoes a transi- 
tion from the s-state to the p-state, the matrix element ¥%,,; vanishes at 
r = 0 and r— ©, while its modulus |“,,;| has a maximum for r of the 
order of a few atomic units of length. The functions w,, and w,; vanish at 
r = 0 and then monotonically increase with increasing r up to the first 
extremum (the position of the extremum depends on /, and k). As Fr Is 
increased still further, the two functions exhibit an oscillatory behaviour. 
At r = 0, and in the oscillatory region, u,, and w,; contribute very little to 
the integral. In order that the value of the integral as a whole should be 
sufficiently large, it is necessary that the first maxima of w,, and u,; should 
lic in the neighbourhood of the maximum of |%,,-,|.. This is achieved only 
at a particular optimum value of /,. 

A modification of the Born approximation has been given in [59]. In 
this work the idea was to calculate the reaction matrix K given by (12.11) 
on the Born approximation, 1.e. assume that 


K =—kit#pokitt (18.18) 
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where By is obtained from B by replacing the functions x,-,- in (10.72) by 
the terms u,d,,- outside the integral. The exchange interaction was 
neglected. We shall denote the approximate reaction matrix given by (18.18) 
by K©® and the approximate amplitude by q®. According to (12.9): 
gq = k7FKO "1 — 1K) ko3 (18.19) 
In point of fact this approximation is a generalization of the formula 
given by (5.18) which was derived for elastic scattering by a centre of force. 
If the matrix elements of K°® are small, then q? ~ k7?K(k7? and we 
have the Born approximation. However, if the interaction is sufficiently 
strong, the Born amplitude given by (18.17) may exceed the limit set by the 
unitarity relation for the accurate amplitude, while the expression given 
by (18.19) will never exceed this limit. Therefore, the latter expression may 
be looked upon as a reasonable extrapolation of the Born approximation 
to the region of strong interaction. It may be used to allow, at least in part, 
for the effect of coupling of the type a > a’, a > a’ ~ a’ and also the com- 
petition between the aa’ and a-— a’ transitions. In this method it is 
no longer possible to sum over all the partial waves by reducing the entire 
sum to a simple expression of the form of (18.8) as on the Born method. 
However, in some cases it is sufficient to take into account a few partial 
waves, and then the summation problem does not arise. When a large 
number of partial waves enter into the calculations it is occasionally possible 
to use the fact that, beginning with a certain value /,, the amplitude calcu- 
lated from (18.19) is practically identical with that found on the Born 
approximation. We then have the approximate relation 


lo 
g=6°+) (o,-¢/) 
=o 


When /, is very large, this method is not particularly productive and 
other methods of estimating this part of the sum must be found. In [60] 
the specific case of the excitation of the hydrogen atom is used to discuss 
in detail the various possible methods of summing over all the partial waves. 

The effect of the atomic field and the exchange interaction on the electron 
can be calculated on the distorted wave approximation. In this method, 
the quantities V,,, W,,, Voq and Wo, are retaincd in the system given by 
(18.1) and (18.2). Moreover, V,,. and W,, are retained in (18.1), but in 


(18.2) they are neglected. In this way one obtains the uncoupled equations 
AF, +(k2?—2V,,)F,—2 i WF, dr’ = 2V,,F,+2 { W.oFodr’ (18.20) 


AF. +(ki—2Voo)Fo—2 | WooFodr’ =0 (18.21) 
The asymptotic forins of F, and Fo are 
eikor 
Fo a5 efKovr sg 
: 
other (18.22) 


Fa = doz 
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Let G,(r, r') be the Green’s function satisfying the equation 
AG, +(k2—2V,,)G,—2 | W_(r,r")G,(r", r’) dr” = 5(r'—r) (18.23) 


and having the following asymptotic form for r — 00 
1 ikar 


G,~——F Lr’) (18.24) 
4n r 
where F, is a function to be determined. We can then write 
F.==2 | G(r, r’) (Vso Fo+ i W.oFo dr") dr’ (18.25) 
In view of (18.24) we have for large r 
1 iker 
F, ~-=" —|F, (YeoFo + + | MoFo ar’) dr (18.26) 
and hence 
1 
doe=— 5- | Fa(YoFe + | WaoFod") dr (18.27) 
Tt 


In order to elucidate the significance of the function ¥, let us express G 
in terms of Go, which ts given by 


1 eikalr—r’| 


Go=- (18.28) 


2n jr—r'| 
The connection between G and G, is contained in the integral equation 
G=G, -| Vid’ G(r, r")Go(r", r’) dr” 
-| Wi A(r’, "G(r, rv" )Go(r", r') dr’ dr” (18.29) 


For large r 


e 7 ikanr’ (18.30) 


Assuming that, for large r, G is given by (18.24) we have 
F (rv) =e ihemr’ _ 9 | V..GoF dr’ —2 i WG )¥ dr'dr” (18.31) 


We thus obtain an integral equation which is equivalent to the integro- 
differential equation 


AF, +(k?—2V,,)F,—2 | W,,F, dr =0 (18.32) 


with the asymptotic behaviour 
ika? 
e « 


F ~ ethan. Qa (18.33) 


As can be seen from these formulae the function F describes the elastic 

scattering of an electron by the excited atom where prior to the collision 

the electron has the momentum k,n. The distorted wave method involves 
E.A.C. 8 


102 THE THEORY OF ELECTRON-ATOM COLLISIONS 


the resolution of F and ¥ into partial waves. The partial amplitude q,,. 
in the distorted wave approximation is proportional to the integral 
- a 
[A. (7) Ky vif, 1”) +| Wis, vil’, rf ir’) dr’ | dr 
0 
where fy, describes the elastic scattcring with angular momentum /, by 
an atom in the ground state, where /,;, represents the elastic scattering 
with angular momentum /, by an cxcited atom in the state 7, /,. As in the 
simplest casc of the Born approximation, g,y 1s not a monotonic function 
of the angular momentum /, of the incident electron. There are spccial 
values of /, which play a dominant role. These values arc somewhat different 
from those which enter into the Born approximation, since the zeros and the 
extrema of the functions f,, are displaccd relative to the corresponding 
points for the functions w;,(r). 
One occasionally uses a simplified version of the distortcd wave approxi- 


mation in which the exchange interaction 1s neglected, so that (18.20) and 
(18.21) are replaced by 


AF, +(ki—- 2a) Fg a 2V,0F 5 


We note that unless factors other than the distortion arc also taken into 
account (for example, exchange and coupling through intermediate states), 
the final result may not in fact constitute an improvement on the Born 
approximation, but may on the contrary lead to an cven larger discrepancy 
between the approximate and true amplitudcs. 

A review of the various versions of these methods from the point of 
view of variational principles is given in [17]. 

When all the terms in (18.1)-(18.2) arc included, the corresponding 
approximation is known as the strong coupling approximation. However, 
this approximation takes into account only the strong coupling O—a@ and 
ignores the effcct of O—a’—«” ctc. coupling. Nevertheless, the cffect of the 
latter may be appreciable, so that it is better to consider a system of a number 
of equations which take into account the coupling within a group of states. 
Both the distorted wave and the strong coupling approximation involve 
extensive numerical work which is usually carricd out by computers. A 
detailed discussion of the various problems associated with the sclection 
of the best computational scheme is given in [61]. 

In addition to numcrical calculations, usc can also be made of variational 
methods. In their simplest form the variational methods are employed 
within the framework of the distorted wave approximation. As was shown 
above, in the distorted wave approximation the transition amplitude is 
cxpressed in terms of functions describing clastic scattcring by the atom 
in the ground and cxcited states. These functions may be reprcsented by 
approximate analytic functions containing paramcters to be varied. There- 
forc, the application of the variational method rcduces to the determination 


(18.34) 
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of the optimum parameters of trial functions describing elastic scattering— 
a problem with which we are already familiar [17]. If the excited state is 
an s-State, the trial function will in general have the same structure as in 
the case of elastic scattering in the ground state. But if the excited state has 
a non-zero angular momentum, then the variational calculation becomes 
so involved that it is much easier to carry out a numerical integration of 
the equation. 

The application of variational methods becomes more difficult in the 
case of the strong coupling approximation. Comparisons of different 
methods, and also the comparison of variational calculations with the exact 
solution (using a model which admits of an analytical solution on the strong 
coupling approximation) are given in [62]. 

Below we give the results of some calculations concerned with the excita- 
tion of the hydrogen atom to different states; these will be compared with 
experimental results. 


19. EXCITATION OF THE 25-STATE 


The excitation of the 2s-state of the hydrogen atom has been extensively 
investigated. It is the testing ground for the various methods which are 
available for such calculations. Thus, calculations have been carried out 
using the modified Born approximation [60], the distorted wave approxi- 
mation [63], the strong coupling approximation (ls—2s) [64], and finally 
calculations with allowance for the 1s—2s—2p coupling [61, 65]. 

It is evident from Table | that the latter calculation involves the solution 
of three equations for Z = 0 and four equations for each of the angular 
momenta L > 0. The calculations yield the excitation cross-section for the 
2s-state and for the 2p-state. Partial cross-sections for 0 < LZ < 4 are cal- 
culated in [65] and for 0 < L < 6 in [61]. The total cross-section can be 
obtained by adding the partial cross-sections for large L, calculated 
on the modified Born approximation. Figure 21, which is based on the 
data reported in [65], shows the excitation cross-section for the 2s-state 
averaged over the spins and computed on different approximations. The 
partial cross-sections for angular momenta L between O and 2 are shown 
separately for the triplet and singlet states in Fig. 22. Inspection of these 
figures shows that the results obtained when the ls—2s coupling is taken 
into account are not very different from those obtained when ls—2s—2p 
coupling is allowed for. The addition of the 2p coupling has little effect 
on the excitation cross-section in the triplet state, except for the region 
near the threshold; but it has a more pronounced effect in the singlet state. 
Moreover, the cross-section near threshold is increased when 2p is taken 
into account, whereas away from the threshold it is reduced. The partial 
cross-sections for large L exhibit a greater relative change when the coupling 
with 2p is added than do the cross-sections corresponding to small L. 

The excitation cross-section for the 2s-state has been determined experi- 
mentally [66]. Direct measurements of the ls—2s transition cross-section 
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[O 40) 30 
E (eV) 
Fic. 21. Cross-section for the excitation of the 2s-state (in units of mao?): solid curves— 
distorted wave approximation (sum of partial cross-sections for L = 0, 1, 2); crosses— 
with ls — 2s — 2p coupling (all partial waves); open circles—strong 1s — 2s coupling 
approximation (sum over L = 0,1, 2); triangles—with Is — 2s— 2p coupling (sum 
over L = 0, 1, 2) 
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Fic. 22. Partial cross-sections for the excitation of the 2s-state in units of zao?: (a)—sym- 
metric case, (b)—antisymmetric case (solid curves—distorted wave approximation; open 
circles—with strong ls — 2s coupling; crosses—with Is — 2s — 2p coupling) 


5- EXCITATION OF THE HYDROGEN ATOM BY ELECTRON IMPACT 105 


are not readily obtainable, and therefore in the experiment metastable 
hydrogen atoms in the 2s-state were subjected to the action of an electric 
field which induced their transitions to the 2p-state. The radiation due to 
the 2p— Is transition was measured, and from the intensity of this radiation 
an estimate was made of the magnitude of o(is—2p). In order to obtain 
the absolute magnitude of o it was necessary to introduce a calibration. 
For example, it was assumed in [66] that at 45eV the cross-section o is 
practically equal to the cross-section obtained on the Born approximation. 
Another method of calibration was used in [67]. 

The results reported in [66] and [67] differ by a factor of about 2. The 
reason for the discrepancy is obscure, and therefore the comparison between 
theory and experiment is somewhat preliminary in nature. Nevertheless, 
there is a good agreement between the calculations of [65] and the experiment 
described in [66]. Figure 23 shows the experimental curve and the theoretical 
results from [65]. 


E (ev) 


Fic. 23. Calculated cross-sections for the excitation of the 2s-state and experimental 
results (solid curve—experimental; crosses—calculation including ls — 2s — 2p coupling. 
The cross-section is given in units of 2ao? 


20. EXCITATION OF THE 2p-STATE 


The cross-section for the excitation of the 2p-state has been calculated 
on the Born approximation [60] by the distorted wave method [68] and with 
allowance for the 1s—2s—2p coupling [65]. As in the calculation of the 
excitation of the 2s-state, the partial cross-sections with allowance for the 
ls—2s—2p coupling were determined for 0 < L < 6 and were then aug- 
mented by the cross-sections calculated on the Born approximation for 
Le> 6. 

Figure 24 shows the partial cross-sections for L = 0,1,2 calculated 
by different methods. The formula given by (11.13) was used to determine 
the amplitudes for the excitation of the 2p) and 2p, , states. The excitation 
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Fic. 24. Partial cross-sections for the excitation of the 2p-state (in units of zap7): 

(a) symmetric case, (b) antisymmetric case (solid curves—distorted wave approximation; 

open circles—-with strong coupling Js — 2s; crosses—calculation with coupling. For 
L = 1, 2, the cross-sections arc reduced, taking 2, = L — 1) 
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Fic. 25. Polarization of the radiation FIG. 26. Comparison of measured and 
cmitted as a result of the 2p > Istransition calculated values for a1: 1—experiment; 
(curve — calculations; points — experi- 2—calculated with 1s — 2s — 2p coupling; 

mental) 3—Born approximation (all in units of a9?) 


cross-sections for these states, o(2p,) and o(2pp)), were then caleulated from 
these amplitudes. In addition the following quantities were calculated: 
1. The total cross-section for the excitation of the 2p-state: 


o(2p) = o(2po) + 20(2p;) (20.1) 

2. The cross-section ¢, determined experimentally from the number 

of photons emitted by the excited atoms in the direetion perpendicular to 
the direction of the electron beam. Thus cross-section is given [69] by : 


a, =0-9180(2p,)+0°2460(2po) (20.2) 


5S. EXCITATION OF THE HYDROGEN ATOM BY ELECTRON IMPACT 107 


3. The polarization of the radiation emitted at right angles to the 
direction of the electron beam [69]: 


= = o(2p;) 
7 E ae 5(2P,) ore 


The calculated cross-sections are given in Table 5 in units of zag. The 
calculated values of P and o are compared with experimental values in Figs. 
25 and 26. 


TABLE 5 

k o(2p1) o(2po) o(2p) omy P 

0:9 0:0683 0:2234 0:3599 0:3853 0:2009 
1 0:0997 0:6005 0:7909 0-8820 0:2835 
1:1 0:1395 0:8149 1:0939 1:2047 0:2799 
1:2 0:1631 0:8539 1:1801 1-2934 0:2667 
1°5 0:2244 0:6423 1:0911 1-1596 0-1800 
2 0:2395 0:3927 0:8716 0:8968 0:0855 


The excitation of states with n = 2, 3 and 4 has been calculated on the 


Born approximation in [70]. The effective cross-sections summed over 


21. EXCITATION OF STATES WITH n > 2 


/and yp for given 7 are given in Table 6. 


TABLE 6 


Cross-sections (in units of ao”) for the transitions 


E 
(eV) 1—2 1—3 1—4 1—5 1—6 

13-5 1-2818 0-1787 0-0509 0-0202 0-0092 
19-5 1-5354 0-2782 0-0999 0-0476 0-0265 
26-5 1-4993 0-2798 0-1021 0-0490 0-0274 
34-8 1-3886 0-2600 0-0951 0:0457 0-0256 


The interpolation method was used to determine the excitation cross- 
sections for states with nm = 5 and nm = 6. This method 1s based on the fact 
that for large quantum numbers nv an energy change of dE is related to the 
change dn in the quantum number by 


dn~nedE (21.1) 


108 THE THEORY OF ELECTRON~-ATOM COLLISIONS 


The excitation cross-section of an atom in the interval dn near n is 

do =o, dn (21.2 
or, on substituting (21.1) 

do =a,n°dE (21.3) 


In the case of ionization of the atom, the probability that the ejected 
electron will be found in an energy interval de may be written in the form 


do =a,dé (21.4) 
It is therefore expected that n°o, will be a continuous function of » which 


goes smoothly into ¢, when n > o0.° Figure 27 shows n°a, as a function 
of the energy of the atomic electron. Negative values correspond to discrete 


—-O2 -Oll O O'| O-e2 
E (eV) 


Fic. 27. Dependence of 3¢, on the energy of the atomic electron (¢ in units of 7@o?). 
The figures against the curves indicate the energy of the incident electron in eV) 


levels, positive values to ionization. In the region of ionization, o, was 
calculated on the Born approximation. The broken curves indicate the 
interpolation region. The excitation cross-sections calculated in this way 
for states with 1 = 5 and nm = 6 are also given in Table 6, but we only 
cover the low energy region. In [70] the cross-sections are calculated up to 
energies of about 1000 eV. 


22. TRANSITIONS OF AN EXCITED ATOM INTO HIGHER LEVELS AS A 
RESULT OF A COLLISION WITH AN ELECTRON 


In the preceding sections we considered transitions of an atom from 
the ground to the excited state produced by electron impact. In some cases 
it is necessary to know the cross-sections for the excitation of atoms from 
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excited states to higher excited states. A number of such transitions have 
been treated on the Born approximation. 

The cross-section for the transition from the 2s-state to all the levels 
with m between 3 and 10 summed over 7 and yp are calculated in [71]. The 
calculations were based on the general formula for values of 7 up to n = 7, 
whereas for n = 8, 9 and 10 use was made of the interpolation method 
described in § 21 in connection with the ls—n transition. The calculations 
in [71] are concerned with the energy range between 4°85 and 120 eV. 
Some of the results are given in Table 7. As is expected, the calculated 
cross-sections are much greater than the cross-sections for the corresponding 
transitions from the ground state. 


TABLE 7 

E Cross-sections (in units of zap) for the transitions 

OV) Dee Dyed Dea 5 Dp 6. eS]. Dees Dea De 16 
4°85 65-019 12-330 4:658 2°312 1°33] 0-843 0:569 0-404 
8-65 49-441 9-454 3:612 1-808 1-047 0:664 0-450 0-320 
13-5 37-667 7150 2°726 1-364 0-790 0-499 0-339 0:241 
26°4 23:703 4-446 1:687 0-842 0-487 0-309 0-209 0-149 
34°6 19-488 3-640 1:379 0-687 0-397 0-252 0-170 0-121 


Transitions from the 2p) and 2p,, states to levels with m between 3 and 
10 were calculated using an analogous method in [72] in the energy range 
between 4-85 and about 120 eV. Some of the results are given in Tables 8a 
and 8b. It is possible to use these data to determine the excitation cross- 
section from the 2p-state averaged over yp: 


o(2p) = to(2po) + $0(2P +1) (22.1) 
TABLE 8A 
E Cross-sections (in units of zao”) for the transitions 
(eV) 2Po — 3 2Po — 4 2Po — § 2Po — 6 2Po — 7 2Po — 8 2Po — 9 2Po — 10 
4:85 74:837 14:426 5:259 2-499 1-409 0-892 0-601 0-424 
8:65 64:529 12:323 4:749 2-347 .1:348 O-861 O0°586 0-417 
13-5 52-482 9°826 3:790 1-859 1:072 0-684 0-461 0:328 
26:4 35-403 6:352 2-408 1:176 0-677 420:-432 0:288 0:207 
34-6 29°727 5:246 1:974 0:963 0-553 0-348 0:234 0-165 


LE eR 
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TABLE 8B 
E Cross-sections (in units of za *) for the transitions 
(CV) 2p, —3 2p, —4 2p, —5 2p; —6 2p: —7 2p, —8 2pi — 9 2p, — 10 
4-85 76:515 13:346 4995 2:512 1:414 0:890 0:600 0-424 
8:65 57:44 10:692 4-052 2-040 1:166 0-740 0-505 0:358 
13:5 45-05 8-181 3:077 1:537 0:883 0-560 0-378 0-270 
26°4 29°278 5-102. 1:891 0-935 0:537 0-342 0-229 0-160 
34-6 24-313 4-173 1-538 0-758 0-435 0-268 0-180 0-128 

0-8 | 2 345 810 20 3040 
E (eV) 

Fic. 28. Cross-sections for the transitions 3s — 4p (curve 1), 3p — 4d (curve 2) and 


3d — 4f (curve 3). The cross-sections are given in units of zap? 


The cross-sections for the following transitions were calculated in a number 
of papers [73, 74, 75] (see Figs. 28-34)T: 

3s — (4s, 4p, 4d, 4f, Sp) 

3p — (4s, 4p, 4d, 4f, 5d) 

3d > (4s, 4p, 4d, 4f, Sf) 

4s — (Sp, 6p) 

4p 5d 

4d — d5f 

4f > (Sg, 6g) 

5(s, ps4, f,g) > O(p, df, 9, h) 


+ The broken curves in these figures represent the dipole approximation. 
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O66: { 2 345 8 10 
E (eV) 


Fic. 29. Cross-sections for the transitions 3s — 4s (curve 1), 4s — 4d (curve 2), 3p — 4p 
(curve 3), 3p — 4f (curve 4), and 3d — 4d (curve 5). The cross-sections are given in units 
of TMAg” 


0-6! 1 2 345 B10 2030 


E (eV) 


Fic. 30. The cross-sections (in units of mao*) for the transitions 3s — 4f (curve 1), 
3p — 4s (curve 2), 3d — 4s (curve 3), and 3d — 4d (curve 4). The scale for curve | is 
reduced bv a factor of 30 


1 2 3°45 8 10 [5 
E (eV) 
Fic. 31. The cross-sections (in units of zao*) for the transitions 3s — Sp (curve 1), 
3p — Sd (curve 2), and 3d — Sf (curve 3) 


05 :10 (582 3 5 810 1520 40 


E (eV) 


Fic. 32. Cross-sections for the transitions 4s — Sf (curve 1), 4p — 5d (curve 5), 4d — 5f 
(curve 3), and 4f — Sg (curve 4). The cross-sections are given in units of ao? 


500 
400 


300 


200 
ISO 


lOO 


E (eV) 


Fic. 33. Cross-sections for the transitions 4s — 6p (curve 1), and 4f— 6g (curve 2). 
The cross-Sections are given in units of ao? 
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10,000 


8000 


E (eV) 


Fic. 34. Cross-sections for the transitions 5s — 6p (curve 1), 5p — 6d (curve 2), Sd — 6f 
(curve 3), 5f— 6g (curve 4), and Sg — 6h (curve 5). The cross-sections are given in 
units of 7a? 


As can be seen from Figs. 28-34, the cross-sections for these transitions 
are much greater than the cross-sections for the corresponding transitions 
from the ground state. The cross-sections for optically allowed transitions 
involving an increase in / are greater than the cross-sections for optically 
allowed transitions involving a decrease in /. 


Chapter 6 


THE COLLISION OF AN ELECTRON WITH A HELIUM ION 


23. THE ELASTIC SCATTERING OF AN ELECTRON BY A HELIUM ION 
AND THE EXCITATION OF THE 25 STATE OF THE HELIUM JON 


THE elastic scattering of an electron by a helium ion differs from the elastic 
scattering of an electron by a hydrogen atom owing to the presence of the 
unscreened Coulomb field oc I/r acting on the scattered clectron. This 
Coulomb field modifies very significantly the scattering process and the 
form of the wave function, and its effect 1s not simply additive to that due 
to the hydrogen atom (the two effects are additive only on the Born approxi- 
mation). By analogy with the seattering of an electron by a hydrogen atom 
we shall begin with the approximation on which the helium ion is regarded 
as a centre of foree. We can then use the theory given at the end of § 2 and 
the problem is reduced to the determination of the phases ,1,. 

In [76] the phase jtg 1s calculated by a variational method, and also 
(at two energies) by a numerical integration of the equation 


d? 2 
Gat (+2 -2w)f=0 (23.1) 
dr r 


2 
wa-=4| Md — dr’ 
r lr—r'| 


where w = 8!x~* exp (—2r) is the wave function for the ground state of 
the He* ion, so that 


The formula for w is 


1 
w= (2 + = es (23.2) 
r 


In [76] the function fis normalized so that it has the asymptotic form 
f = F+tanyG, r+ © (23.3) 


where F is a regular Coulomb function having the asymptotic form 


1 
F~sin (ir ~5 +-In 2kr +1) (23.4) 
and G is a non-regular function with the asymptotie form 


1 
G ~ cos (41 — : + Fs In 2kr +1) (23.5) 
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According to (2.68) and (2.71) the scattering amplitude g may be written 
as the sum 
2in 


e | 
q=4ct+ ae (23.6) 


where q> is the Coulomb scattering amplitude. As a measure of the departure 
of the actual scattering from the purely Coulomb scattering, the authors 
of [76] used the quantity 


2 
Ra (23.7) 
dc 
From (23.6) 
a ee e 
= |}1 +— e*"—_(e**#—-1) (23.8) 
Qc 2k 
When exchange is taken into account, (23.1) is replaced by 
d*f 
a Pai + (ke +: + — - 2) f= 2] we r') fr) dr’ (23.9) 
In the present case 
t -—?2 t _—2r’ k? I , 
Wir,r') =+32re “'r’e or ee ak (23.10) 


Equation (23.9) was also solved numerically for two values of the energy. 
Moreover, the phase z+ was found with the aid of a variational method. 
The trial function was taken in the form 


f=F+(atbe*")(1-e-7)G (23.11) 


The trial function for the electron plus ion system containing the inter- 
electron distance r,;, was also employed. The quantities R, and R_, 
calculated with allowance for exchange, were averaged using the formula 


R=4R,+4R_ (23.12) 


The results of the variational calculation were found to be strongly dependent 
on the particular method employed, especially at low energies in the sym- 
metric case. This suggests that the trial function was not very satisfactory. 
The numerical values of y are given in Table 9, and a comparison of the 
results of the numerical integration and the variational calculation is given 
in Table 10. The function R(@) at an energy of 48-84 eV is shown in Fig. 35. 
A calculation of the excitation of 2s and 2p states of the helium ion with 
allowance for the 1s—2s—2p coupling is reported in [77] and [77a]. The 
calculation is quite analogous to that for the hydrogen atom. The partial 
wave equations are similar to the equations derived in Section 11, except 
that the functions u, and v, are the corresponding continuous-spectrum 
wave functions for the Coulomb field. The presence of the unscreened 
Coulomb field has an important effect on the behaviour of the amplitud 
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Fic. 35. <A plot of the function R(O) for E = 48-84eV: 1—without exchange and 
polarization, 2—with exchange, 3—with exchange and polarization 


TABLE 9 
ft ue 
E (eV) with oan Kohn’s Hulthén’s 
exchange method method 
3-28 0:571 0:959 —0:28] 0:419 
8:25 0:554 0:824 0-038 0-404 
15:74 0:531 0-807 0-073 0:365 
24:90 0:509 0:762 0-110 0:333 
48-4 0-460 0-645 0-240 0-315 
65°49 0-443 0-604 0:226 0:307 
TABLE 10 
ry . s 
without exchange : ie 
E _ Numeri- Numeri- 
(cy); OPENED: yee. cal Kohn's Hulthén's cal Varia- 
cal tional tional 
Gales sats integra- method method integra- male 
; tion tion 
48-4 0:468 0-460 0:288 0:240 0:315 0:648 0-645 


65:49 0-443 0-443 0-283 0-226 0-307 0-604 0-604 
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Fic. 35a. (Reproduced with permission from Burke et al. [77)) 


near the threshold. In particular, it is found that as the energy tends to the 
threshold value, the excitation amplitude increases without limit as k7}. 
This is a consequence of the behaviour of the normalization factor Co, 
given by (3.30), as k > o. It is the reason why the excitation cross- 
section ¢ = 4n(k,/ko)|g|* remains finite when k = 0. 

Figure 35a shows the excitation cross-section for the 2s and 2p levels 
of Het based on the data reported in [77]. 


Chapter 7 
THE COLLISION OF AN ELECTRON WITH A HELIUM ATOM 


24. THE WAVE FUNCTION FOR A THREE-ELECTRON SYSTEM 


THE form of the wave function with a given spin and the required symmetry 
with respect to the spin and coordinate variables for a system of three or 
more electrons is not so apparent as in the case of a two-electron system. 

There are a number of ways of constructing such wave functions in the 
form of a suitably symmetrized sum of products of coordinate functions 
and spin functions. The various methods of doing this are reviewed in [78] 
from the point of view of group theory. 

Consider the wave function ®(r,, 0,,..., ry, vy) Of an N-electron system. 
It should be antisymmetric with respect to the interchange of any pair 
(r;, 0;) and (r,, o,). We shall require that it should be the eigenfunction of 
the square of the spin angular momentum, so that 


$?@ = s(s+1)® (24.1) 
V. A. Fock [79] showed that an antisymmetric function satisfying (24.1) can 
be expressed in terms of the coordinate function ‘P(r,,..., rglfeag, -++> ty) 


and the spin function y(o,..., 0,|O,44,---;0y)- Henceforth we shall write 
WU, ...,A]K+1,..., NV) and v(l,..., A|/+1,..., N) instead of ‘Y(r,|ry) and 
v(o,|oy), or simply VW, y, whenever this does not introduce an ambiguity. 
The vertical line separating 1,...,k from k+1,..., N represents the division 
of the electrons into two swarms with opposite spin components. The 
number of electrons A which make up the swarm 1,..., & is related to the 
spin quantum number by 


N 
k=>-s (24.2) 


The functions and y have the following symmetry properties with respect 
to the exchange of the arguments: 

1. The function Y is antisymmetric with respect to the interchange of 
the arguments before or after the vertical line; this may be expressed in 
terms of the exchange operator P; ; in the form 


(1+P; )¥ =0, Lij<k or ij2ek (24.3) 
2. WI, ..., AJA +1,..., NM) exhibits cyclic symmetry, i.e. 


(1 - y Pry) ¥=0 (24.4) 


fuk 
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3. The function y is symmetric with respect to the interchange of the 
arguments before or after the vertical line, 1.e. 


(1—P, )x=0, ij<k or ijek (24.5) 


Consider the expression 


_ [ki(N—k)! 
® oF >» (—1)*¥(A,, sah sets veey Ay) 
N! Aiyerer> Ak 
K (Ag, ey Aglage rs +03 Aw) (24.6) 
where 4,,..., Ay is a set of N unequal numbers from among the sequence 
Li 2 cater AV 
1,..., N, and A is the parity of the permutation ( 1 ; ) In view 
“] 2 ooe ~N 


of (24.3)-(24.5), the function (24.6) is antisymmetric with respect to the 
interchange of r;, o; and r,, o,. Next, using the well-known formula 


S*?= N-4N?7+Y of (24.7) 
k<i 


which relates the operator $? to the spin exchange operator P{?, and bearing 
in mind (24.2) and (24.4), it can be shown [79] that '¥ is an eigenfunction 
of S? with the eigenvalue s(s+1). If x is chosen to be of the form 


2 us snas PO gd aonaea On (24.8) 
where 


x= Og, 1 and B = Og. es | (24.9) 


then ® will also be an eigenfunction of S, with the eigenvalue s. The 
coefficient [((N—A)!k!/N!]* in (24.6) ensures that when and y are nor- 
malized to unity, the function ® is also normalized to unity. 

In the case of the three-electron system there are two possible spin states, 
corresponding to s = 4 and s= 3%. According to (24.2) the coordinate 
function corresponding to the first case is ‘¥(1|23), while the wave function 
corresponding to the second case is ‘W(|123). The properties expressed by 
(24.3) and (24.4) can then be written in the following form: 


(1L+P,,9)'¥({123)=0, a8 =1,2,3 (24.10) 
(1+P3)¥(1[23) = 0 (24.11) 
(1—P,,— P,3)¥(1|23) = 0 (24.12) 


The problem now arises as to how one can construct the two functions 
Y¥(|123) and ‘Y(1|23) so that they will have the properties given by 
(24.10)-(24.12). It turns out that when an arbitrary function of three variables 
(123) is symmetrized in accordance with Young’s scheme, we obtain a 
function having the properties of ‘¥(1|23) or ¥(|123) depending on the 
form of the scheme.f This can be seen most easily in the case illustrated 


{+ This was shown for an arbitrary number of electrons in [80]. 
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in Fig. 36(a) [57]. This scheme involves the action of the operator 
(1-—P,,)(1—P,;)(1—P,3) on the function g(123). It is easy to see that the 
function 


Q, = (—P,2)U— P13) — P23) (123) (24.13) 
satisfies (24.10) since, for example, 
(1+P,)U-Pi2)=( — Pt) = 0 


It is also possible to carry out the symmetrization in accordance with the 
scheme illustrated in Fig. 36(b). This corresponds to the application of the 
operator (1—P,,)(1+ P,,). Consider the function 


@2 = (1—P23)(1+ P12) (123) (24.14) 
It satisfies the condition given by (24.11) since 
(1+ P,;)1—P23) =0 
It also satisfies (24.12) since (as can be shown by direct calculation) 
(1—Py,—P,3)(— P23) 4+ P12) 9123) = 0 (24.15) 


, (a) (b) 


Fic. 36. Young’s scheme: (a) for the function g, = yw(l1 23); (b) for the function 
G2 = w(I|2 3) 


It follows that the functions ‘(1 ]23) and ¥({123) can be obtained from an 
arbitrary function of three variables by means of the Young symmetrization 
scheme of Figs. 36(a) and 36(b) respectively, namely, ‘¥(|123) = @, and 
(1/23) = 3. 

The above analysis is concerned with the general case of a three-electron 
system. Since we shall be dealing with the collision of an electron with a 
helium atom, we shall consider three-electron functions (1 | 23) and ¥(| 123) 
of a special form, constructed by combining two-electron and one-electron 
functions. Since in the ground state the spin of the helium atom Is Zero, 
the only function which is of interest in practice is ¥(1]23). The function 
‘¥(|123) can only be useful in the case of transitions between triplet states 
of the helium atom as a result of a collision with an electron, and we shall 
therefore confine our attention to the function ‘Y¥(1 | 23). 

The function ¥(1|23) can include two forms of two-electron functions, 
namely, symmetric w(1|2) and antisymmetric y(|12). This is a consequence 
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of the law of addition of angular momenta. The state of a three-electron 
system corresponding to spin 4 may be formed by combining a one-electron 
state with both a singlet and a triplet state of a two-electron system. This 
can also be expressed in Young’s terminology: the scheme shown in 
Fig. 36(b) may be constructed in two ways, as shown in Fig. 37. Accordingly, 
there are two expressions for ‘¥(1|23), namely, 


 (1[23) = c,[W(1|2)F(3)—W(1|3)F(2)] (24.16) 
Y (1|23) = co[({12)F(3) — W(]13)F (2) + 2W(]32)F(1)] (24.17) 


3 


(a) (b) 


Fic. 37. Two ways of constructing Young’s scheme for the function g(1|2 3): (a) addition 
of a cell to the scheme for yw (|2 3); (b) addition of a cell to the scheme for y(1|2) 


It can be easily verified that these expressions exhibit all the properties which 
the function ‘¥(1 | 23) should show. A general formula for the construction 
of (2+1)-electron functions from n-electron and l-electron functions will 
be given in Chapter 8. The expressions given by (24.16) and (24.17) can be 
obtained from this formula as a special case. The coefficients c, and cy, 
are the normalization coefficients. Suppose that F refers to the continuous 
spectrum and 1s normalized as a plane wave, i.c. 


| F*(k, r)F(k’, r) dr = (22)°d(k—k’) (24.18) 


Suppose further that w(r,|r,) and W(|r,, r.) refer to the discrete spectrum, 
and are normalized to unity: 


[ worlra)? ar, drs,= 1, [ Wri ra)? dry dr, =1 (24.19) 


If it is required that ‘V(r, |r2, 73) should also be normalized as a plane wave, 
1.¢. 
i ww. dr, dr, dry = (2n)2d(k—k’) (24.20) 
then we must set 
l 


1 
ee Ty ee 


The connection between these coefficients and the normalization of the spin 
functions will be established in § 25. 


(24.21) 
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25. EQUATIONS FOR THE WAVE FUNCTION OF THE 
SCATTERED ELECTRON 
Starting from equations (24.16) and (24.17) we can write the coordinate 
wave funetion for the electron plus helium atom system in the form 


P(r, |ro, r3)= 3S (Wari fra) Farts) — Wari |rs)Fai(r2)] 


I 
+=) (WaClri, r2)F 42(r3)—Wa(|ry, r3)F42(r2) 
</6 A 


+2, (|rs, ro)F 42(r;)] (25.1) 


where w,(r, |r.) and w,(|r,,r2) are the coordinate functions respectively 
of the singlet and triplet statcs of the helium atom. In the simplest approxi- 
mation, they can be constructed from the one-electron functions: 


YW a(l|2) = YU e(2)+(2)e() 

Wa({12) = WD e(2)— 9(2)Y(1) 
The functions F,,(r) deseribe the state of the seattered elcetron. The sub- 
script A represents all the quantum numbers charaeterizing the state of the 
atom, and a = |, 2 represents the two possible combinations of the spins 
of the atom and the incident eleetron. 


The equations for F,, and F,, may be obtained from relationships 
analogous to (10.24): 


[ valra|r2)(H-E)(r |rasra)dr, dr, = 0 (25.2) 
[ vara r2)(H-S)¥ (ry |r2, 73) dry drz = 0 (25.3) 


where H is the energy operator for the electron plus helium atom system 
and is given by 


l I l 2 2 2 

H=— 4A, —4tA, —4A, +— +—4+-— avrg. Gy oe, pe . rx, = |ri-ry| 
Piz My3 0 o300dCOPy PnP 

(25.4) 


The energy & of the system is made up of the energy of the atom prior to 
the collision and the energy of the incident electron. 

Let us substitute cquations (25.1) and (25.4) into (25.2) and (25.3). 
We shall take into account the symmetry properties and normalization of 
War, fez) and Ww ({ry, 72), and will assume that the latter two functions 
satisfy the cquations 


2 Z 1 
[ vilra (= 4a. — 482 = - 2 + a) varifraddey dry = 0 
My OP, Ir,—rs| 
(25.5) 
Joan r,)dr, dr, = 0, 
A’ #A (25.6) 


Z Z l 
| vidoe (- JA, —4A,-— ---- + 


ry 5 Ir; —r| 
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These equations are satisfied if w, and w, are the exact atomic functions, 
but they may also be satisfied by approximate atomic functions. 
Let 


E =tkaat Eg (25.7) 


where E,,, is the energy of the atom and 3k4, is the energy of the incident 
electron.t After some simple rearrangement, we have the following system 
of equations for F,,, 


AF ga(¥3) + k4aF galt’) = 2» [Vaaralt's)a' F ggl¥3) 
+ | Waaraa(" 1) 3)F geq(¥') dr, | (25.8) 


where the potential energy matrix V, ,, 18 given by 


2 
Vases) =— Sy 42 | 24D ay, 25.9) 
3 


lrys—r,| 
with 


PAA‘'1 =| war [radar (ro) dry 
Pasa =| Wari ra Waris ta) dra 


The exchange operators W, 4/49’ are of the form 
Waarratasts) = — | WaCra|rs) (HS yba(ra|rs) dra (25.11) 
Waarratists) =— V3 { WiC |r2)(H- alts, ra) dr. (25.12) 
Waaaa('asts) =— V3 | Wari r2)(H-Oa(rs|rz)dr, (25.13) 


Waara(tsts) =— | WM|rist2)(H—-Sba(|rayraddr2 (25.14) 


The expressions for the exchange operators can be transformed further. 
To begin with, consider the operator —4A, in H. Since this operator 
is self-adjoint, we can replace w ,(r, |r.)A; F4,(7,) in equations (25.11)-(25.14) 
by Fy,(r)A,Wa(r;|r2). A similar replacement may be carried out in the 
terms containing the function w,(r,|r2). Next, if w4(r,[r2) and wy({ry, ro) 
are the exact wave functions for the helium atom, we can employ the con- 
ventional equations for them. As a result, the exchange operators assume 
the form 


I 2 
Wa araa(" 1; r3) = (34: + a =? =| Padtaa’(" 1; r3) 


31 1 
+ | Seeernete te) (25.15) 


ri2 


(25.10) 


+ When yw, are the approximate functions, the quantity E,, is the mathematical 
expectation of the energy of the atom calculated from this function. 
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where 
Paarsa(tasts) = | WAC |raWa(ra|rs) dre (25.16) 
paaaatists) = | Waris ra Wa |ras rs) drs (25.17) 
Paasaltists) = V3 | WRC |raWal|ras 03) drs (25.18) 


Paar, 13) = J3{ wary, ro )Wa(ra|r3) dry (25.19) 


Paans(hstast3) = WarifraWa(ra|rs) (25.20) 
Paan2hysF25%3) = V3 0X(ry [ra Wa |ra. rs) (25.21) 
Paari(has tests) = V3 We(|ry, r2)Warg|rs) (25.22) 
Paar? ¥as 03) = Waris te a|ras Fs) (25.23) 


In practice the equations must be simplified. In addition to the errors 
which appear as a result of the simplification and which are also encountered 
in problems concerned with the hydrogen atom we have here one further 
source of difficulty, namely, the fact that the atomic wave functions are only 
approximate. This leads to errors in V and W. 

The equations (25.8) must be augmented by the asymptotic conditions. 
We shall assume that 

eikaar 


F ,4(r) om ee Aer SA Aoay tT Aa» r—-> © (25.24) 

r 
Let us now consider the asymptotic behaviour of the complete wave 
function for the system (including the spin variables). We return to (25.1) 


and consider any one of the terms representing two discrete states of the 
atom in the sum over A: 


1 
Vary ro, r3)= V3 [Wari |ra Farts) —Wari [ts F412) ] 


1 
* 7g hall ro )F 42's) — Wally. ts)F 42(r2) 


+20 4(|r3, r2)F 42(r1)] (25.25) 


In addition to this term, the complete wave function also includes 
Y 4(r3[ry, r2) and V(r, [r3, r,). According to (24.6) we have 


1 
Dr, o1,%202,1303) = CV alr ara, t3)B 1224 
‘\ 


+P 4(r3|"1, r2)B 3% %2 +P 4(r|rs,1)B2%1 3] (25.26) 


Asr, — ©, only those terms remain in which r, cntcrs into F,,; the remaining 
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terms all vanish, 1.e. 
2 
Wry (Po, r3)~ V6 Wars, r2)F 42(r)) 


1 1 
W(rslri, ts) oo TeV Aalrsira)Fas(a) = gv allrs> r2)F 42(r1) 


1 1 
P 4('2|I35 ry)~ 9 W (ra |t3)F ai(11) 2 gv Alas r3)F 42(r1) 
If we collect these terms and substitute them into (25.26), we obtain 
(when r,; > o) 
a3 B,— a2 Bs 


J/3 Or, 01,202,303) ~ Fairy Wa(ra|rs)o la 


| = a 
+ Far Wallrs, r>) 23 2 %3 — B3 01 %2 ~ B2%1 &s (25.27) 


V6 


The functions 


1 
i= iri 3 — B32) (25.28) 
and 
= [2B exp ety —0r(0ts B+ 9 B,)] (25.29) 


are the spin functions for the electron + helium atom system which are 
normalized to unity (normalization is ensured by choosing c, and c, in 
accordance with (24.21)). The function y, describes the state of the system 
in which the atom has zero spin, and the z-component of the spin of the 
incident electron along a special direction is +4. The function yx, describes 
the superposition of the following two spin states: 

(a) The triplet state of the atom with the spin z-component equal to 1 
and the state of the incident electron with the spin z-component equal 
to —4. 

(b) The triplet state of the atom with the z-component equal to 0 and 
the state of the incident electron with z-component equal to +4. 

Therefore the wave function for the singlet state of the atom constructed 
in accordance with (24.6) is 


j 
Py2rl'2 02,1303) = Waltalts) 75s B.—%, B3) (25.30) 


The expression 
P42(12 02,1303) = Wa(\ro, 13) X3 (25.31) 


is the wave function for the triplet state of the atom, constructed in accordance 
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with (24.6) with the spin z-component of the atom equal to |. The expression 


3B y+ 
O'42(¥202, 1303) = Wal|r3,%2) ———= a : (25.32) 
s/f 2 


is also the function for the triplet state, but with zero spin z-component. 
In view of (25.30)-(25.32) and (25.24), we can rewritc (25.27) in the form 


eikairi 


J 350 (1101, %2 02,1393) ~ Dar dar % 


Fy 0 ry 
2. eikazrt | eikaqs 
+0 42942B Jp? + O424a2% J (25.33) 
3 Fy 3 Fy 


The current density along the r, axis in configuration space yields the 
cross-section for the excitation of the singlet state 


do, = ldas| dQ (25.34) 
0 


The excitation cross-section for the triplet state with spin z-component 
equal to | is 


dof) = 5942 Naya)? dO 


Whereas the excitation cross-section for the triplet state with zero spin 
z-component is 


1k 
do) = — 42 2 dQ 
2 3 ko Cre 
On summing over all the possible spin components in the final statc, we find 


that 


Kao 
do, =—~|qa2 
ko 


2 dQ (25.35) 


is the excitation cross-section for the triplet state with arbitrary spin z-com- 
ponent after collision. 


26. ELASTIC SCATTERING OF AN ELECTRON BY A HELIUM ATOM 


When polarization is neglected the elastic scattering of an electron by a 
helium atom is described by 


AF +(k?~2V)F = 2] WF dr’ (26.1) 


where V is defined by (25.9), (25.10) and W by (25.11) or (25.15) through the 
ground state function wW(r,|r,) (the subscripts of V, F and JW are omitted 
for simplicity). 

The function W(r,|r.) may be representcd approximatcly as a product 
of the onc-electron functions 


Wry |r2) = u(r, u(r) (26.2) 
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We then have 


t 2d t 
yen? 2 (HO (26.3) 
r |r—r'| 
The kernel of the exchange operator in (25.11) is 
2:2 1 
Wr’, r) = u(r) 40, + $A, += +—— 7 + 4ko + Eo 
pte [poe 
’ 4 Z 1 J a” 
+] u(r"){ 4A” + — — ——~ — ——— _ J u(r") dr" | u(r’) (26.4) 
r” |r —r"| [r—r"| 
Using (25.15), we obtain 
y) 1 2/.u d uw 
Wr) = u(r) [38,4 K3 +S — [OE ry 26.5 
[ror] J [eel 


The expressions given by (26.4) and (26.5) are, in general, different. However, 
they become identical in one important special case, namely, when «(r) 
satisfies the equation for the self-consistent field 


2 u* 
Au+2(-— | ——— dr’ Jut+22u =0 (26.6) 
r lr—r'| 


and if instead of the exact energy E of the helium atom, equation (26.4) 
contains the energy calculated by the self-consistent field method, i.e. 
2 2(..! 
P25) = | OE) ay (26.7) 
ror 
Equation (26.6) may be used to reduce (26.5) to the even simpler form 


W = u(r)u(r’) (34: roar ro) (26.8) 


which 1s analogous to the kernel of the exchange operator in the case of 
scattering of an electron by a hydrogen atom in the antisymmetric case [sce 
equation (13.2)]. 

If the function u(r) does not satisfy the equation for the self-consistent 
field, so that equation (26.4) cannot be reduced to (26.5), the problem 
arises as to which of them should be used in calculations. The variant 
represented by (26.4) was used in the calculations considered below [81, 82]. 
However, it may be argued that (26.5) is preferable. The point 1s that (26.4) 
and (26.5) differ considerably at small r: in (26.4) the function W(r'’, r) 
includes the term 2/r, whereas (26.5) does not include this term. When 
the replacement (r,,r3) > (r’,1r) is introduced into the exact expression 
(25.15), the function W(r,,r3) does not contain 2/r. Thus equation (26.5) 
is closer to the exact exchange operator than (26.4). The term 2/r in (26.4) 
appears only because the wave function 1s not exact. 

Equations for the partial waves can be obtained by separating the 
angular variables in (26.1). Numerical integration of the equations for the s 
and p waves using the atomic functions for the self-consistent field was 


128 THE THEORY OF ELECTRON-ATOM COLLISIONS 


carried out as far back as 1933 [83]. The dependence of the phases 6, and 6, 
on the wave number which was obtained in that paper is illustrated in 
Fig. 38, together with the values of these phascs without allowance for 
exchange. Morse and Allis used 69 and 6, to determine the angular dis- 
tribution and the total scattering cross-section; thcy compared these with the 
experimental data existing at that time (see Figs. 39, 40). It 1s evident from 
these figures that, in general, the theoretical calculations reproduce quite 
satisfactorily the observed dependence of the cross-section on the energy 


[80 


90 


6 (deg) 


0 1 2 3 


Fic. 38. Phases for an electron scattered by a helium atom: 1, 3—with exchange; 2, 4— 
without exchange 


and the angular distribution, including the characteristic backward scattering 
maximum at 180° at low energies, which is similar to that observed in the 
case of the hydrogen atom (§ 17). Figure 39 can be used to dctcrmine the 
scattering length. It is found to be a2 1-5 which corresponds to 
o = 4na* & 28ae. 

The scattering length was recently recalculated by numerical integration 
of the equation for s-waves of zcro energy [82]. The one-electron atomic 
function u(r) was chosen in the form 


u(r) = N(e~" + Ce 74") (26.9) 
where WN is the normalization factor and 
C= 0-6, Z = 1:4558 (26.10) 


This function is a good analytical approximation to the function for the 
sclf-consistcnt field: the calculation yields a = 1:442, corresponding to 
c = 26:lag. The cross-section at zero cnergy obtained by cxtrapolation of 
the experimental data was found to be o = 22. 

Thus the numerical calculations in [82] lead to a result which cxcccds 
the cxpcrimental valuc by only 18%. Part of this discrepancy is due to the 
fact that the polarization of the hclium atom by the scattcred clectron was 
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O {@) 20 30 40 
E (eV) 


Fic. 39. Elastic scattering cross-section for an electron scattered by a helium atom (in 
units of za 7): curve—theoretical, points—experimental 


not included. Some of the error is due to the approximate nature of the wave 
function for the helium atom, including the fact that the approximate 
functions W(r, |r.) given by (26.2) do not contain the distance r,, between 
the electrons and therefore the correlation between the atomic electrons is 
ignored. 

The variational method was used in connection with the scattering of an 
electron by a helium atom in [81]. The one-electron atomic function was 
taken in the simple form 


uy =Ne* (26.11) 


O 30 60 90 120 I50 I80 
@ (deg) 


Fic. 40. Angular distribution of electrons scattered by a helium atom at electron 

energies of 6eV (curve 1), 20 eV (curve 2) and 30 eV (curve 3). The broken curves 

represent theoretical calculations without exchange, the solid curves the theoretical 
values with exchange allowed for, and the points the experimental values 
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27 Z3\t 
Z=—, N= (=) (26.12) 


The exchange operator was determined with the aid of (26.4), using the 
approximate expression for the energy 


where 


E, = —2Z? (26.13) 
The trial function was chosen in the form 
f=sinkr+(at+tbe“)(1—e7-“) cos kr (26.14) 


The phase 6,(k) obtained from these calculations is in good agreement 
with the results of numerical integration in spite of the fact that different 
one-electron atomic functions u(r) were used in the two cases. 

The scattering length determined by the variational method was 
a = 1-46, quite close to the result obtained by numerical integration. 

Finally, the phase 6g was determined in [84] by an approximate solution 
of the integral equation for the s-wave. The resulting scattering length was 
found to be somewhat lower than that obtained by the variational method. 


27. EXCITATION OF THE 2!),S- AND 2°S-STATES OF THE HELIUM ATOM 


The excitation of the 2'S- and 2°?S-states of the helium atom was cal- 
culated in [85] by the distorted-wave approximation using the variational 
method; in [84] approximate expressions were obtained for the distorted 
waves with the aid of integral equations. 

Here we shall confine our attention to the calculations described in [85]. 
On this approximation, the equations for the wave function of the scattered 
electron are 


AF) +(k2—2Vo9)Fy —2 | WooFodr’ =0 OTD 


AF, +(k2—2V>,)F,—2 | W,,F,dr’ =2 (V20Fo “ | Wo Fo dr’) (27.2) 


AF, +(k}—2Vy1)F 1-2 | W.,F,dr' = 2 | We Fodr' (27.3) 
with the asymptotic behaviour 
; ikor eiki.2r 
Fy =e" +49 = rig dia (27.4) 


The function Fy describes the scattering of an electron when the target 
atom is in the ground state; F, and F, describe the excitation of the 2°S- 
and 2'S-states respectively. Following the distorted-wave method, we 
neglect W,,F, and W,,F,, which connect the excitation of the triplet and 
singlet states, and the effect of the functions F, and F, on Fy. As can be 
seen from (27.3), excitation of the triplet state 1s only possible as a result 
of the exchange interaction JV¥,,9. The quantities V,, and }¥’,, are determined 
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from (25.9) and (25.11). The wave function for the ground state is taken 
in the form 


Wo(r|r2) = Uo(T,)Uo(2) (27.5) 
where wy is given by (26.11). The wave function for the triplet 2°S-state 
is taken in the form 


Wilh, r2)= N [ey (ry )uo(r2) — wary) (r2)] (27.6) 
with uw, and w, given by [86] 
u,(r) =e 7 
u,(r)=re "’—ce7™" hy) 
where 
c= 2. pf =0-6l, eal eS 
and 
ety)” 
(B+a)* 


The coefficient c is chosen so that the function u,(r) is orthogonal to 
u(r), and the parameters f, y determined from the condition for minimum 
energy. The function for the singlet 2’S-state is taken as 


Wo(r,|r2) = Niu, (ry)ua(r2) +u2(r)uy(r2)] (27.8) 
where u,(r) and w,(r) are in the same form as for the 2°S-state [equation 
(27.7)] in order to simplify the calculations. The parameters of the one- 
electron functions u, and uw, are different for the triplet and singlet states 
when determined from the minimum energy condition. 

The function w,(r,|r.) given by (27.8) is not orthogonal to Wo(r, |r2). 
As a result VY, includes terms which fall off as r~* when r > oo. In reality, 
V9 should fall off exponentially when r > oo. For this reason, terms which 
behave as r~! when r > oo are rejected in V, 5, and this removes the most 
undesirable consequence of the fact that w,(r,|r.) and w,(r,|r2) are not 
orthogonal. There remains the error in the numerical value of the terms 
having the correct asymptotic behaviour, which 1s due to non-orthogonality. 
This error was not estimated in [85]. Only the s-waves were taken into 
account in the expansion of F, into the partial waves. Therefore the problem 
reduces to the determination of the three functions 


o(r) = rF (7) 
Y1,2(r) = rF¥ , 2(r) 


where F,(r) and #,(r) are the solutions of (27.2) and (27.3) without the 
right-hand side terms; /,; and /, have the asymptotic behaviour 


(27.9) 


1 
fy - e= sin (k,r+Nq) (27.10) 


The amplitudes qg, and qg, are expressed in terms of fo, 7; ,and ”, by formulae 
analogous to (18.27). 
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The latter three functions are determined by a variational technique. 
The trial funetion fp is chosen to be analogous to (26.14). When the 
asymptotic behaviour given by (27.10) is taken into account, the trial function 
becomes 

ino 


toe ma +a’) ?[sinkyrt+(at+tbe *)(l—e “)coskgr} (27.11) 
0 


where Z = 27/16. 
The trial functions for #, and /, are taken to be 
in 


yi = —— (1443+ Ce™) sin kyrtA,(t—e ")coskyr] (27.12) 
1 


and similarly for #,. Calculations show that (27.12) is a very poor approxi- 
mation for “, at low energies near k, = 0, and therefore 1, and q, are 
ealeulated for k, > 0-3. The corresponding effective cross-section is shown 
in Fig. 41. For 1, and g, the calculation is carried out with k, > 0-1, and 


4 


go (107 * 7708) 
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Fic. 41. Partial cross-sections for the excitation of the 2'S (curve 1} and 2°S (curve 2) 
states of helium at / = 0. The energy is measured from the threshold 


a determination is made of the scattering length for the 2° S-state of the atom 
(although the numerical result is not given in [85)). 

In order to illustrate the importance of exchange, the parameters n, 
and q, have been calculated without introducing this effect. Figure 42 shows 
the dependence of 4, 4, and 4,, on & without exchange. The broken 
eurve shows the part determined on the basis of the estimated scattering 
length, i.e. the slope of the 7(A) curve for & = 0. As can be scen from the 
figure, even with relatively small changes corresponding to k 2 0-5, the 
order of magnitude of the phase and its dependence on & is very nearly the 
same with and without exchange. However, ncar the threshold (k ~ 0-1-0-2), 
the inclusion of exchange has a radical effect on the variation of the phases: 
the scattering length in thc 2?S-state increases very rapidly as a result of 
exchange. The phase 4, approaches 32/2 when k = 0-05, i.e. very near 
the threshold. In othcr words, seattering by a hclium atom in the 2?S-state 
is very close to resonancc scattering. This is associated with the existence 
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of a negative metastable He state which lies about 0-1 eV above the 
2°S-state of He. 

The calculated excitation cross-section in the 2°.S-state near the threshold 
is shown in Fig. 43, which also gives the experimental values [87]. As can 
be seen, the calculations show a much sharper peak near the threshold than 1s 
observed in experiment. This is not surprising when it is recalled that 


Fic. 42. The effect of exchange on the phases of trial functions: 1—phase 4,2 with 
exchange, 2—72 and 3—”, 
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Fic. 43. The variation of the cross-section for the excitation of the 2°S-state of helium 
near the threshold: 1—calculated [75], 2—experimental [77] 


scattering by a helium atom in the 2°S-state is very near to resonance 
scattering. Under such conditions, the phase 7, and therefore the calculated 
cross-section for the 2°S-state, are very sensitive to small changes in the form 
of the wave function inside the atom, so that the variational method using 
simple trial functions is unsatisfactory. No exact numerical calculations 
have as yet been carried out. 


E.A.C. 10 
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A possible interpretation of the behaviour of the 2?S-state excitation 
cross-section near the threshold is given in [88], where it 1s suggested that 
inelastic electron scattering giving rise to excitation of the 2°S-state is 
analogous to nuclear reactions which proceed through the formation of a 
compound system (in the present case He in the 25, state). If this is so, 
the effective cross-section for s-scattering near the threshold should be 
described by the Breit-Wigner formula 


7 ror, 
CSS ee Se Ee 
ko (E-E,)? +409 +1) 
where E, is the energy of the compound system, I’, is the width for the decay 
of the ion into an atom in the /S-state and an electron, and I, is the width 


for the decay of the atom into an atom in the 2°S-state and an electron. 
The widths [) and T, may be written in the form 


ie = Done rj — 2k, 7 (27.14) 
IT, and I, are the so-called reduced widths and are independent of energy. 


(27.13) 
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Fic. 44. Cross-section for the excitation of the 2°S-state of helium on the compound 
system model: curve—Breit-Wigner formula with the parameters given by (27.15), 
points—experimental values [77]. Excitation threshold Ej, = 19-79 eV 


Figure 44 shows a plot of a(£) as given by (27.13), with the following 
empirical parameters 


E, = 20:19 eV 
yg = 22.107? eV.cm (27.15) 
y= 25.107? €V-cin 
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These values of the parameters are not the only ones to give adequate 
agreement with experiment. Another selection of the parameters is also 
discussed in [88]. 

For comparison, Fig. 44 also shows the experimental points, and as can 
be seen the agreement is quite good. The magnitude of yg , corresponds 
to a width T ~ IL eV or a lifetime ~ 107 '° sec, which is of the same order 
as the half-life for the auto-ionization of He in the 2s7'S) state (about 
2:5x 107 '° sec). The experimental o(£) curve exhibits one further feature. 
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Fic. 45. Cross-sections for the excitation of the 2°S-state of helium: 1—experimental, 
2—calculated for/=Oand 1. E is the energy measured from the threshold 


It is found that near the threshold the values of o lie very nearly on a straight 
line and therefore in this region o may be regarded as a linear function of 
the energy measured from the threshold: 


o ~(E—E)) 
However, it follows from the general expression o = (k,/ky) J |q| 7dQ that 


if g remains constant for k, = 0, and this is certainly the case for neutral 
atoms, then o should be proportional to k,, i.e. 


oO™~ Ce 
It appears that the region near threshold where q is approximately indepen- 


dent of energy and o ~ ,/E—E, is very narrow and is not detected owing 
to the low resolving power of the instruments. 
Direct comparison of the calculations of the partial cross-section for 
= 0 with the experimental data at energies of a few eV above the threshold 
has no meaning because the partial wave with / = 1 introduces a considerable 
contribution at such energies. The contribution due to the p-wave was 
estimated in [85] with V,, and IV,, neglected. The results are given in Fig. 45 
together with experimental data. 


28. THE 2!S—2°,S TRANSITION DURING A COLLISION WITH AN ELECTRON 


The transition from the metastable 2'S-state to the lower 2°S-state in 
which the excess energy is given to the incident electron was discussed in 
10—2 
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[89], in which the system of equations for the s-waves was numerically 
integrated on the strong coupling approximation. The equations were 


q2 ie 6) ® 
4 + (k?— =2% 22)f2-2 | Wo2f,dr’ =2 | Wash dr’ (28.1) 
0 0 
d*f, 2 ’ r r 
ip? > + (Ki - 2% yy )fi- af wi ,dr =2] W,.f,dr (28.2) 
ar . 
0 0 
with the asymptotic behaviour 
l ; 
tf, a 7, mn kor+q, ear ti ad qy eikir (28.3) 


2 


where f, describes the incident electron with energy 343 and f, describes 


the electron after collision, with energy 4k] = 4k3+AE. The partial cross- 
section for the process 1s 


Kk, 
CG = 4nlqi| (28.4) 

k, 
Y 4, and #,, were calculated from the general formulae (25.9) and (25.11) 
using the analytical expression (27.7) for the atomic functions. Table 11 
gives the results of the calculations, with the cross-sections o in units of ao. 


TABLE 11 
ky je (eV) Co 
0:044 0:026 180 
0-050 0:034 141 
0-060 0-048 113 
0-100 0-134 64-7 
0-300 1-210 14-6 


The calculated cross-section at 0-026eV is smaller than that observed 
by a factor of about six. This discrepancy may be due in part to the 
fact that the functions for the 2'S- and 1'S-states are not orthogonal. 
However, this would hardly lead to an error of this magnitude. Another 
possible explanation of the discrepancy is the formation of the metastable 
He™ ion. However, no quantitative estimates of these effects are available. 


29. EXCITATION OF SINGLET 7'P-STATES (11 = 2, 3, 4) 


Excitation of the singlet P-states was analysed in [90] on the strong- 
coupling approximation, but without allowance for exchange effects. 

In order to elucidate the relative importance of the various partial waves, 
a calculation was made of the cross-sections corresponding to total angular 
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momenta L = 0, 1,2 and 3. In the group of indices v’ = (n’, /{, /3) referring 
to the initial state, /; = O for all the transitions under consideration while 
> = L. In the group of indices v = (n,/,,/,) referring to the final state, 
/, = 1, while /, can assume the values |L—1|</, < ZL+1. When in 
accordance with the laws of composition of angular momenta there were 
more than two equations, they were split into independent pairs. Tables 
12-14 give the calculated values of the effective excitation cross-sections for 
the levels 1'S—3'P (Ey, = 23-308 eV), 1'S+2'P (Ey, = 21-211 eV), and 
11S = 4'P (E,, = 23-734 eV) in units of 2a5/100. As can be seen from these 


TABLE 12 
ky E (eV) F710 Oo. G21 O12 O32 023 o 
0:05 23:342 0-446 0-0°20+ 0-172 0:0°15 0-:0446 — 0-649 


0-10 23-444 0-245 0-0743 1-52 0-0°35  0:0756 0:0°37 1-808 
0:20 23-852 0-166 0-210 3°52 0:0712 0-606 0-0412 4-503 
0-40 25-484 0-132 0-430 22 0-0726 0-745 0-0712 4-031 
0-60 27:204 0-102 0:366 2°12 0-0732 = =1-135 0-0747 3-731 
0-80 32:012 0-074 0-298 I52 0:0728 1-286 0-:07978 3-191 
1-0 36:908  0:0503 0:233 1-02 0:0717 = 1-220 0:0139 2°539 


TABLE 13 
ky E (eV) C10 To. O21 Oj2 
0:05 21-245 1-13 0:0°43 0-036 — 
0:10 21°347 0:976 0:0713 0:29] 0:0524 
0:20 21-755 0-601 0:024 3:72 0:0°68 
0:60 25-107 0-274 1-28 9-16 0-0°43 
1-00 34°811 0-128 0-810 4-44 0:07420 

TABLE 14 
ky E (eV) T10 Co1 O21 O12 
0:05 23:°768 0-139 0:0°25 Q:250 0:0°61 
0-10 23°870 0-092 0-011 1-075 0:0710 
0:20 24-278 0-075 0:541 1-210 0:0723 
0:60 27:°630 0-050 0-172 0-820 0:0727 
1-00 37°334 0-026 — 0:406 0:0717 


1 Superscripts indicate number of zeros after the decimal point, for example 
0-0°43 = 0-00043. 
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tables, at the energy of only 1 eV above the threshold, the main contri- 
bution to the excitation cross-section is due to the partial cross-section ¢ 2. 
In some regions the latter is much greater than all the remaining cross- 
sections. The origin of this effect was discussed earlier in § 18. 


TABLE 15 
To1 T10 
ky EX) S.-C. d.-w. S.-C. d.-w. 
0-1 23-444 00°43 0-044 0-245 0-233 
0:6 27°204 0°366 0:346 0-102 0-096 
1:0 36:908 0-233 . 0-186 0-050 0-040 


The relative importance of strong coupling may be estimated from 
Table 15 which gives the cross-sections for the excitation of the 3'P-level 
(in units of 2a$/100) calculated by the distorted-wave (d.-w.) and the strong- 
coupling (s.-c.) approximations. 


Chapter 8 


GENERAL EQUATIONS DESCRIBING THE COLLISION OF AN 
ELECTRON WITH AN ARBITRARY ATOM OR ION 


30. THE WAVE FUNCTION FOR THE ATOM + ELECTRON SYSTEM 


CONSIDER an (1+1)-electron system consisting of an n-electron atom (ion) 
and an incident electron. According to § 24, the coordinate part of the 


wave function P(r), ..., | Peaps+++s "n41) IS antisymmetric, i.e. 
(+P, )¥ =0 (30. 1) 
ifi,j <k or i,j > k; it is also cyclically symmetric, i.e. 
n+l 
(1 — >» Px. } T= 0 (30.2) 
i=k+1 


The number of terms on the left of the vertical line is related to the spin s 
(which is assumed to be fixed): 


k = a = (30.3) 


By analogy with the treatment of the helium atom, we shall construct 

the function 

WC isceosl Fey qgteainay) 
from the products of n-electron atomic functions Ww and one-electron functions 
F, so that (30.1) and (30.2) are satisfied. 

For a given square of the spin angular momentum of the (7+ 1)-electron 
system, there are two possible spin states of the sub-system of n-electrons, 
namely 

Ss, =s—4, Sy=st4 
The corresponding functions are 
Wis Sifiiti.n) and W2(1,....jalj2+1,...0) 
According to equation (30.3), 
Jp =tn-G-HHk,  jg=jn—-Gti)=k-i (30.4) 

We shall therefore consider the products of the n-clectron function w 
and one-electron functions F of the form 

Wil, ..., kK[k+1,...,n)Fy(nt+1) 
and 
Wal, ....k—1k,...,.mFi(n4+ 1) 
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Consider the sums 


n+l 
Pi= Yo a: Pina Wi, k|R+1, mF y(n 41) (30.5) 
i= 1 
and 
k nti 
Y =) Db, iP, bP ine V2(1,. K— Ak, . nF, (n+ 1) (30. 6) 


t=1 1=k 
with the coefficients a; and 5, ; chosen so that (30. 1) and (30.2) are satisfied. 
It 1s understood that application of the exchange operator results in the 
following transformations 


Pp fs 2 eres ee eee) SS Ig) poets Dp res Oy sno, 1D) 


and 


PJ 1s sees cindy Di weg NOs seag onus Dyce) 
=o | 0 errno epee (| are 7 (7) (ener a meee / emma: 1g 
i.e. the a-th and b-th arguments can be exchanged independently of their 
position within the brackets. 
We note that all the steps necessary to obtain a; and 6, ; may be performed 
in compact form by using the commutation relations 
PoP ee = = Pie P as = pal be (30.7) 
These relations can easily be verified. For example, 
Pap Py-f(abe) = P,,f(acb) = f(bea) 
Pac Pas f(abe) = Pa. f(bac) = f(bca) 
The coefficients a; and 5; , are defined unambiguously to within a common 
factor which is unimportant and can be included in the normalization. 
We shall not go into details but will merely write down the final result and 
then verify (30.1) and (30.2) [91]. The expressions for ‘’, and , are of the 
form 


eG; (1 = Panes) Vall vy k[K + l,.... mF (n4t]) (30.8) 


k-1 n 
W2=¢2(1- 5 Pus }[t- y Punsrt(1=2k+2)P rats 
1=1 i=k+1} 
x wo, Lig k— 1k, 5 M)F3(n+1) (30.9) 


The coefficients c, and c, can be defined from the normalization conditions. 
Suppose that w, and w, belong to the discrcte spectrum while F, and F, 
belong to the continuous spectrum, and that w, , arc normalized to unity, 
but that F, and F, are normalized as a planc wave to (2n)°?d(kK+4’). If 
we require that VY, and Y, should also be normalized to (2x)*d(A +4’), 
then 


c, =(n—k+1)73 


30.10 
cy =k 4 (n—2k +2)? +(n-—k+1)] 7? ( ) 
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We shall now show that , satisfies (30.1) and (30.2) (the normalization 
factors will be omitted from the proof). 

1. The anti-symmetry condition 1s clearly satisfied up to the vertical 
line, since the arguments of ‘¥, and , up to the line are identical. 

2. In order to establish the anti-symmetry condition after the vertical 
linc let us apply the operator (1+P;, ;) and use equation (30.7). The result 
is 


(1 - > Puner)(L+ Pi Wa Fa Lj AntI 
¢ +P, i), _ l=k+1 
— 2 Pinel + Pid Fi, j=ntl 
which is zero because 
(+P; )¥, =0 
3. The cyclic symmetry may be established as follows. From (30.7) it 
can be shown that 


n+l n n n 
(1= S) Pas)(1- 3 Panes) =(1- Ye Parse )(1- 3 Pa) 
j=kti i=k+1 jak j=kt+i1 


it follows that 


n+l 
(1 - \ Pas) Wi =0 


j=kt+1 
The proof for ‘V, is quite similar. 
1. The anti-symmetry condition up to the vertical line can be proved by 
operating with (1+ P,, ;) on Y,: 


k-1 n 
(1-5 Pa] 1 - » Pine +(0—2k+2)Prnei 
t=1 joHkt+1 
XL +P; Wh 2 Fo, ljAnt+l 


(1+P; )'P, = k-1 
(1 — yy Pa] |(1 = 2k -+2)Pp we (hPa, 
| t=1 | 
| —Y Paneih + Pid] ¥aF jJ=nt+l 


which is zero because 
(1+P, bh. =(A+P, 2 = 0 
2. When /, 7 # A, the anti-symmetry condition after the vertical line is 
obviously satisfied. When j = & we have 


(1+ P; 4)?» = = > cool) 7 > Pan Ok) Enea (+P, )W2F2 
ttl] i 


Which is equal to zero because 
(+P, )W2=0 
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3. The cyclic symmetry condition may be verified by using the properties 
of W,. After some transformations we have 


(= Si nuee(Hof ron) [l- 8) 


=), Pe ‘(1 — d ie : — (n=2k-+2) WF, 

jJzk+ 
Next, it may be shown that 

n k-1 n 

(1 - > Pus) Pua(t y Pas) Ww. = (n—-2k4+2)y, 
jukt+i t=1 j=kti 
and therefore , 
n+l 
(1 = Pas) be — 0) 
J=kt+1 
From the general formulae (30.8) and (30.9) we have, as a special case 

of the expression used in § 24, 


W (123) = i CHC |2)F 3) — W(1|3) F2)] 


(123) = == Lv(|12)F(3)—¥(|13)FQ) +20((32)FUY] 
N 


We can use the above results to write down an expression for the coordi- 
nate function of the atom + electron system in terms of the atomic and 
one-electron functions. Let A be the set of numbers indicating the energy 
of the atom, its orbital angular momentum and the z-component of the 
angular momentum. We shall write Y as the sum 


~ =F (Pat ha) (30.11) 


Each of the terms V4, 2 willconsist of expressions of the form of (30.8) or (30.9). 
If we omit the common factor of ¥,, 3, which will be unimportant in 
the ensuing analysis, we have 


Yat a (1 = Y Paaet)barFan 


i=k 


(n—2k+2) “( K- 
P53 = Jie ee k {—)>P 
o ( n—k+l1 : d cs 
«fi- > Passe 2420s] VaFa (30.12) 
i=k+1 


31. EQUATIONS FOR THE FUNCTIONS F,, 


In order to establish the equations for Fy,, we shall use a relation analo- 
gous to (10.24), Le. 


[ WiCH-S)¥dr,,..,dr,=0,  a@=1,2 (31.1) 


in which // is the energy operator for the electron + atom system. 
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We shall suppose that w%, is the exact atomic function so that 


(Ho — Eg Aa = 0 (OE2) 

The exact functions will be replaced by the approximate atomic functions 
only in the final expressions which will be of the density-matrix type. 

There is another way of deriving the equations for F,, [92, 93], in which 
it 1S supposed from the outset that w,, are the functions for the self-consistent 
field. 

Let us now write the energy operator for the atom + electron system 
in the form 


H = Ho—4taAyi1t+V(r1, -s traa) (31.3) 
where V is the energy of interaction of the (7+1)th electron with the 
nucleus and the remaining n electrons; it is given by 

i | 


— we rrei—ti| 


(31.4) 

Mn+. i=l 
Substituting the expression for ‘Y and H given by (30.11), (30.12) and (31.3), 
(31.1) and using the orthogonality and symmetry of w, we have, in view 
of (31.2) 


(A+k4)F 4a = 2 y VaA'a Fata Sa Wied F ara’ (31.5) 
A’ A'‘a’ 
where 
Vaara(tn+ 1) =| WaaVWaadry,...,dr, (31.6) 
is the interaction energy matrix, 
ka = 2(6 — E4,) (31.7) 


and W is the exchange operator of the form 


W jaean’F ya Paes) = Poa | Wall 1 age |g, AG +k4g—2V) 


KW ara lng isco [ers hE a(t) dry, dt, 
+ Goo | Waal! 1, gad is r,) (Ay + k4rq— 2V) 
XWarg(1 1, a nee Pog dE geal 6) ai sna UP, (31.8) 


Pao’ 2Nd gaq are given by 
Pir = Pi2 = Par = 9 


(31.9) 
Poo = (K—1)(n—2k +3)(n—2k42)7? 
Qy=n—k 
di. = AA(n—k)(n—2k+3)(n—k 4-1)? [(n—2k 4-2)? tn-—k4+1]73 
n—k+1 7? 
dar = (nk) Lee rarreesd 


do, = —(n—kKk+ 1)(n—2k+2)7! 


Using the explicit expression for ”, the symmetry properties of w,, and the 
fact that A is self-adjoint, the potential energy matrix and the exchange 
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Operator may be reduced to the forms 


Z  [ Paarar dd = 
Vig = |--+| “48 | 
AAta +{ 1 (31.10) 
Wasree Fara = [Waren DF yele) dr’ (31.11) 


, / Z J S 
Vos (34 ee ) Paaee( 1” 
ee 
-| er rra Cake r’, r) 


ir 


dr” (31.12) 


where the p represent expressions of the density-matrix type. 

In order to write these expressions in the most compact form, we shall 
introduce a special symbol for the integral of the product W4,W,’,° over all 
the variables with the exception of a specified few. For example, the integral 
with respect to all the variables r,,...,7, except for r; will be denoted by 
the symbol 


(ri adldata =| Waal! Ara’ dry, ..., 0-1 ANi4 4,-0.5, (31.13) 


If, after integration, r; in W,, 1s replaced by, say, r, while in w,-, it is replaced 
by r’, this will be denoted by 
+ r’ 
sn ae 
In this notation we have 


r’ r’ 7) r’ 
alr jyosk +(n—k 
Paaiit) (") () mn 1 an ). 


(31.14) 

HN . (7) to-ken(T) (7) (31.15) 

Paani yr an(" ). (; )., (31.16) 

Paar2"'s sonal ), (; .. (31.17) 

Paarl. N=aai(} ) . )., (31.18) 

neriednon()(")..tra("),()n ens 
pacer =aule() (05), 


" y’ r” r 
Ok DA... ") ae 
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;! yr” r 
Paanr yr’ = 412 |- D( anos am 


ne eo ee ad 
Pe Pat ar \leR Un/ a2 r PaJ A1\ln-1 n/ A’'2 


pasts st)= 001 [&-D( es m 


” ; 7 4 ‘f 
r’ or r’ or r’ r r r 
Pe Va A2\le Pn Py-1 TaJa2\"n-1 "n/a 
yy 


y! a r 
ror n= k—-1 
Paar", ) = 422 ( (c »( 4m 
+0" 4 (" "| | 
Va-1 i A2 Pa-1 rn A‘2 
rr" ror 
4: k—2 ( 
P22 [ 14 ay My P27 a2 


yr’ re" r r” 
+n ke1)(" an eS ao! (31.23) 


Using the expression for the exchange operator we can reduce (31.5) to 
(A+ kia Faq = 2), VanaFarat2 > | Waaraat WF yar’) dr’ (31.24) 
A’ A’a’ 


The asymptotic expression for F,, in the case of a collision with a 
neutral atom 1s 


eikaar 
ik gars 
VA, ~e - Ong Anas UG Acaoda r (31.25) 
while in the case of a collision with an ion it 1s 
: Zo Zo 
if kaa + —In(kK4al — Kaar) i( kaart = In 2k ar) 
KaAa 
E ig gee < G AgagAa a ne (31.26) 


where k,, = K,,V. 

The system of integro-differential equations (31.24) subject to the 
asymptotic conditions (31.25) or (31.26) may be transformed to a system of 
integral equations. Let G,, be the Green’s function which in the case of 
a collision with a neutral atom is of the form 

1 eikaalr-r'l 


Gua —_— (31.27) 


2m lr—r'| 
while in the case of a collision with an 1on of charge Z it satisfies the equation 


(34 4+ 4k? 4 2) G,,=—dr—r’) (31.28) 
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By analogy with (10.53), consider the kernel 


Z 
K 4 a'aa’ ee Gadr, r’) (Ze aoe 54) Oa’ -| G(r, r’) Wada’ r’) dr” 
r 
(31.29) 
(in the case of collisions with neutral atoms Z = 0). We then have 
Fae = P4aSstoSay + | Kavaoltsr Fae(e)de’ (31.30) 
A’a’ 


where @,, represents either exp (/k,,r) in the case of collision with a neutral 
atom, or else the Coulomb function with the asymptotic behaviour 


Z 
[ck 2 ee i(kaar+ “In 2kaar 
Dag Auge [kA0 ee (K Aa Kaa | Be g®, : (31.31) 
It is known [57] that this function can also be written in the form 
le ALN. se Een. 
p= eT I-iF e'" F ~ 2 1 ikr—ke) (31.32) 


where F is the degenerate hypergeometric function. 

From equation (31.30) it is possible to obtain the expression for the 
amplitude g4..,44. One must consider the asymptotic behaviour at large r. 
The Green’s function G,,(r, r’) for large r and finite r’ has the following 
asymptotic properties: for collisions with neutral atoms 

ik Ag? 
G~- oe ee (31.33) 
4n r 


whereas for collisions with ions 


k 


G~——e Aa Daa (31.34) 
4n 


l i(k sar ap 2k Aer ) 


where 
nZ ; 
nZ TX... Z 
7 = ek r(1 +if) a | - —< I, (kr + kr) (31.35) 


Substituting (31.33) or (31.34) into (31.30) and bearing in mind (31.25) 
and (31.26) we have for collisions with neutral atoms: 


4 AgagAa — ~ ae cae 5. Vaara F giglt ‘) 
+ de Wanaalt', rFawl) dr’ (31.36) 


and for ee with ions: 


GW AgagAa — a — PaAt ) Maes Oaa’ FB yg(t) 


p Wasa st’ )F yodr’) ir’ dv +qQ),, 40a,544o 5009 (31.37) 


Sd 
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Consider the asymptotic form of the complete wavc function, including 
the spin variables. By analogy with § 25, we shall scparate out the terms 
which are finite when, say, r,4, 2 00 (we shall thercfore be concerned with 
that part of (30.11) which contains the atomic functions for the diserctc 
spectrum). After some transformations 


‘ oF, ( Vn—k+1 w’,, 

~) Matin SSS SSS ) Tati )Xn 

= AL: Ail +1 $40 Jin ak 2) " iad a F 42 +1 at | 
n—2k+2 


ey 
Qi ey =e > Dao Far ns DBat1 (31.38) 
where a numerical factor has been omitted. In these expressions, ®,, is 
the atomic wave function with spin s — 4 made up in accordance with 
(24.6) from W,, and x = fy, ..., Bus Oa ps ++s Gs Py. 1S the atomic wave 
function with spin s+4 made up in accordance with (24.6) from w,, and 

= Bi, ..-5 Burt, Us ++) M3 and @,, is made up in accordance with (24.6), 
also from 2, but with y given by 


X = (n= 1y4(1 + d Pay) Bo nas ig Onan sony O,, 
fakt+ 


(Di. corresponds to spin s+34 but the spin z-component for this function 
is s—4 and not s+}). 

As an cxample, consider the collision of an electron with an atom whosc 
spin prior to the collision is zero. In view of (30.3) we have from (31.36) 


ikaw 
els aor € 
D~ O,,1e me +) Mar Daorai 0+ [ALY ea tioa ret 


4 etka 
ta > V4 Gacta2z—— Fe pi (31.39) 


k+3 4 


where we have omitted the subscript »+1. Henee, when there is no change 
in the spin of the atom, the aris cross-section Is 


UO Ast A1 = 


: dai at}? dQ (31.40) 
Aol 
while if the spin change is allowcd, we have on summing over the final 
spin states 


k 

U6 Ag. a2 = 7 |daoraa|? dO (31.41) 

Agl 
If the atom has a non-zero spin, say, S9, prior to collision, there are two 
possible values for the spin of the system, namely, s, = s9+4 and 
S» = Sg—}. This leads to two different problems. In one of them the 
initial state of the atom is described by a function ®,, while in the other 
it is described by ®, 5. In order to obtain the observed cross-section, the 
result must be averaged over the spins. Considcrable simplification is achieved 
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by considering the valency electrons of the atom and the incident electron 
as moving in the field of the nucleus, and fixing all the internal electrons. 
In this approximation we shall usc the symmetrized function for the system 
consisting of the valency electrons of the atom plus the incident electron. 
The stmplest example of this kind is the monovalent atom. According to 
V. A. Fock, the wave function for the valency electron satisfies the equation 


[-4A+ UK) Wet Wu, = EW, (31.42) 


where U is the potential energy of interaction of the valency electrons with the 
nucleus and the inner electrons and ie given by 


Z 2p(r,r ) 4 
dr’ 


U=-——+ (31.43) 
r \r—r'| 
while Wis the exchange operator defined by 
Ww, = [Se uN? 5, (31.44) 
Equation (31.32) may be written in the shorter form 
Hew, = EW. (31.45) 


The coordinate part of the wave function for the system consisting of 
the valency electron plus the incident electron will satisfy the equation 


H(r) +H,(r') + —| Y= EV (31.46) 


yr" 


By analogy with (10.46) we shall write Y in the form 
i Sl d [Fa War’) + Fa(r War) (31.47) 
From 
[veo [Hat + Hate) + a | Yoda’ =O 
we have the following equation for F, 
(1H, — =) Fr) + 2, Ve ver) + > | Hr Wh Fr) dr’ =0 (31.48) 


where 


V.. =] var Wale’) (31.49) 


aan 


We = (E, 4+ TY oxa (31.50) 


The opcrator }V,,- describes the exchange interaction between the 
incident elcctron and the valency electron. The cxchange interaction with 
the inner electrons is described by the operator Win /,. 
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If we let F have the asymptotic form 
gikar 


Pg butte (31.51) 


then by analogy with the case of the hydrogen atom, the scattering cross- 
section assumes the form 

do* = |qz,|° dQ (31.52) 

Separation of the angular variables in (31.24) or (31.48) is carried out 

exactly as in the case of a collision of an electron with a hydrogen atom 


($ 11). 


The equations for the partial waves are of the form 


d? +1) 
pa aa ae Ki 7 < Saat =2 > V para Sarat 
r r rue 
+2 a W gata sT)faav(r dr’ (31.53) 
A’a'l’ 
0 
where ¥, yay’ aNd W 4 gaa’? aLe given in terms of Vgyrao aNd Wy ying BY 
V gaa = | Y(O'O' Wawel’ ¥ry(0'G') AQ (31.54) 


W anroatt 0) = 00" | YEO" OWar aol # Yr y(0' 9) AQ dQ” (31.55) 


Owing to the presence of the term [r’ +r”"|~* in Wy 4-aq(r', r”), the expression 
given by (31.55) will be different for r <r’ andr > vr’. The two cases will 
be represented by W“) and W, respectively. The exchange term on the 
right-hand side of (31.53) may be written in the form 


{ Wf dr’ = | (WO) — WR) fdr! + i WF, dr’ (31.56) 
Q Q 0) 


where y represents the group of subscripts Aa/. Equation (31.53) cannot 
readily be solved numerically because of the presence of the integral 


\° 6) 


i 2 , : é - 
| ye J, dr’, and successive iterations are necessary. However, if we use 


Q 
approximate atomic functions which are expressed in terms of the one- 


electron functions, we can write 
WY i » Ri) Qiv(7’) (31.57) 


This may be used as in § 14 in order to transform equation (31.53) into a 
form which can be solved without iteration. 


Chapter 9 
COLLISION OF AN ELECTRON WITH A SODIUM ATOM 


32. ELASTIC SCATTERING OF ELECTRONS BY SODIUM ATOMS 


As an example of the application of the general theory developed in the 
preceding chapter, let us consider one of the simplest problems, the collision 
of an electron with a sodium atom. In accordance with the discussion in 
§ 31 thts problem may be looked upon as involving the motion of two 
electrons, i.e. the incident and the valency electrons in the field of the nucleus 
and of the inner electrons. The latter comprise ten electrons in the ls, 2s 
and 2p shells. 

Among existing calculations of the inelastic scattering of electrons by 
sodium atoms, the most accurate is that reported in [94]. This work is 
concerned with the scattering length, making allowance for exchange and 
3s—3p coupling. It turns out that this coupling plays a very important 
role. Using experimental data on oscillator strengths, the authors of this 
paper conclude that 99-4°% of the polarizability of the atom is due to the 
3s—3p coupling. Therefore, when the single virtual transition 3s—3p is 
included, this accounts almost completely for the polarization of the atom 
by the incident electron. 

The equations discussed in [94] are, in fact, identical with (31.48), 
although they are derived in a different way. The corresponding equations 
for the partial waves on the total angular momentum representation are 
of the form 


ee) 


(HE) Ake fr ae » Wry Iver aT > | Woyet Sve dr’ _ 0 (32.1) 
v y’ “ 
1d? 1(l+1) : 
(hoes eS ‘ / 
HOS et aa OCW (32.2) 
and WS y ty ds) 


The function U(r) and the exchange operator Jf’f can be expressed in 
terms of the densities p(r, r) and p(r, r’) constructed from wave functions 
for the closed shells using (31.33) and (31.34). In the present case 


] | 
pr, r’) = — Ry (r)R,,07") + — Ror) R2(7) 
An An 


1 
1: Dy Rap t)Rapl") Y*(0,~)Y,,(0',0') (32.3) 


p= 
where R,,, R, and R2, are the radial wave functions for the atomic electrons 
in the closed shells. 
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Using (32.3), we have the following expressions for U and W 


11 

U=—-~—+4?2 [ (RE, + RI +3RE Dior r')r’? dr’ (32.4) 
r a 

0 
Ve | (rr) Ry(r')F(r) dr’ Ry 
= r')r r 
PI+1 VIa r, Is r r 1s r) 
0 


2B | Mr)’ RCP) F(r') dr’ R2,(r) 


0 


a 


31 
Fy | Manas Rag VFO) dr’ Rap”) 
2 
0 


(315 +3) 
(2l,+3) J 


Yi +iC Tr Ra(r')F(r') dr’ R2,(r) > (32.5) 


where 


= 


The quantities ¥,,, and ,,-, which enter into (32.1) can be expressed 
in terms of the wave functions for the valency electron in the ground 3s 
state or the excited states. 

As was pointed out above, calculations of the scattering length can, 
in practice, be confined to the single 3p excited state. 

The wave functions for the electrons in the closed shells R,,, R,, and R2, 
and the function for the valency electron in the 3s and 3p states were taken 
from [95, 96]. 

The system of equations given by (32.1) can, on the present approxi- 
mation, be reduced to two equations for f399,9 and /311,9 (the subscript A 
which characterizes the state of the atom is now replaced by the triad of 
quantum numbers n, /, 4). In order to find the scattering length we must 


set k3y = 0. We then have k3, = —x? = —0-137. The asymptotic behaviour 
Of fyo0,0 and f3144.9 aSr > is 
feo, o(r) ~ r—a*; fait,orvwe O (32.6) 


In [91] the equations for f3o99,9 and f314.9 were solved by the method 
of successive approximations. At first the exchange operator W was neglected 
11—2 
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in the operator H, so that Hy; became the energy operator Hf, on the central 


io @) 


field approximation. The cxchange term | wfar’ was also neglected. 
0 


The zero-order approximation, f°), to the functions f was then substituted 
into the exchange term. 
The solution of the equations with the exchange term | Wf dr gave 


0 
the first approximation f“? which was again substituted for fin the exchange 


term, and so on. The values of a* and a™ given in Table 16 were obtained 
after five successive approximations (# represents the order of the 
approximation). 


TABLE 16 
n 0 l 2 3 4 
at —2°4 10-0 14-0 13-9 _— 
a7 —2°4 —9-6 — §-3 9-5 --6'6 


Next, the effect of the exchange interaction with the ls and 2s shcll on 
the function f3499,9 was determined in a similar way by substituting the 
function /$$o,.9 into the opcrator IW’. 

The final values of a= were found to be a* = 12,a~ =—9. The cross- 
section averaged over the spins could be determined from known a™ and a7: 


{o> = n[3(a-)* +(a*)*] 


It was found to be 380za3. Moreover, it is possible to use these values of 
a* and a™ to determine the amplitude when the scattering process is accom- 
panicd by the spin flip of the valency clectron. 


Na( t )+e({) = NaC} )+e(t) 


Let ,, and ji,, be the spin z-components of the valency and incident clcctrons 
respectively. Spin flip involves the transition 


(fs, Its) i (fly, — [ls2) (32.7) 


Since [l.,;+f,. must be conserved, this transition can only occur when 
Ms + Hls2 = 0. 

Lct us cxpress the scattcring Icngth a(figy. flea, Shey. —Ps2) for the 
process defincd by (32.7) in terms of a* and a~. Instcad of the symbols 
+ and — we shall write (0) and (1) in accordance with the magnitude of 
the total spin. We then have 


Ads 4, — fst, = fiaviGa) — Cake C1 nian (32.8) 
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Substituting the values of the Clebsch-Gordan coefficients we obtain 
a(fig4s — Meas — Hats Her) = F(a — a) (32.9) 


The corresponding effective cross-section is ¢ = 440nag, which is con- 
sistent with the experimental result ¢ > 260za6 [97]. 


33, EXCITATION OF THE SODIUM ATOM 


The excitation of the 3p level of sodium was considered in [94] on the 
modified Born approximation (§ 18). In the equations for the partial waves, 
only the matrix element V, 5, tyen.4)’.4. Was taken into account, Le. the 
3s—3p interaction. In this matrix element only the dipole part was retained. 
The latter 1s proportional to 


l 
(Yi)35. 3p = oe } R3,R3,r° dr + ff R3,R3,rdr (33.1) 
0 r 
This quantity was then replaced by its asymptotic form 
| 3 
(y)35, 3p a ey R3,R3,r dr (33.2) 
0 


It would appear that only the long range part of the intcraction was 
taken into account 1n this particular calculation. However, this long range 
part plays the dominating role in the case under considcration. 


80 
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E (eV) 


Fic. 46. Cross-section for the excitation of the 3p-state of the sodium atom (in units 
of mao?): 1—calculated, 2—experimental 


Figure 46 shows the numerical results of the calculation and the expcri- 
mental data as given in [94]. 


Chapter 10 


EXCITATION OF MULTIPLY-CHARGED IONS 


34, EXCITATION OF THE 25 AND 3s STATES OF HYDROGEN-LIKE IONS 


WHEN an electron collides with a multiply-charged ion, the unscreened 
Coulomb field has a considerable effect on the electron. It is then convenient 
to use an approximation in which the Coulomb field is taken into account 
exactly, whereas the effect of the electrons in the electron shell is regarded 
as a small perturbation. This approximation involves the replacement of 
Fy,’ in (31.37) by oy. If, in addition, the exchange interaction is also 
neglected, then the excitation amplitude becomes 


1 es 
LAca > | Pa VjdoPAo dr (34.1) 


Calculations based on this formula were carried out in [98] for transitions 
to the 2s and 3s states in hydrogen-like ions and for optically allowed 
transitions on the dipole approximation. 

Let us consider the excitation of hydrogen-like ions. The matrix element 
V sq, 18 then given by 


WiaoW 


rr 


dr’ (34.2) 


awe = 


where Wy, and yy, are the wave functions for the initial and final states of 
the ion. We separate out in (34.1) and (34.2) the dependence on the nuclear 
charge Z. The functions w,, and w, contain the quantity 


p=Zr (34.3) 


as an argument, while the normalizing factor is proportional to Z?. There- 
fore, V,,, may be written in the form 


Vaas = ZU 44,(P) (34.4) 


The function @m, contains as its arguments the quantities Z9/k, and kyr 
where Z, is the charge on the ion, Le. 


Zs 

PaA= Pa ne ky. (34.5) 
“A 

If we now express r in terms of p in accordance with (34.3) and substitute 

_k 

2 
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p (34.6) 
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we have 
Zo 
Pa=Pa 7? Pa-P (34.7) 
PA 
For hydrogen-like ions Z) = Z—1, so that 
J 
eens 
ZL 
Pa = Pa\—_—>Pa-P (34.8) 
Pa 
Substituting (34.4) and (34.8) into (34.1) we have 
1 a 
Lod Paos Pa) =~ 5-33 | PaU aa Pao dP (34.9) 
The cross-section is then given by 
1 p, l 
I 4oa = 73 7, 2 |M|7dQ (34.10) 
where 
J | 
1—— 1—— 
x Z Z 
M = | 3. » Pa P U 449 PAg\ > Pay: P| ap = (34.11) 
: A AAo Pao 
| OF 
sat 
N 
31S 0-5 
Coe ee 
O 90 180 
6 (deg) 


Fic. 47. Angular distribution of electrons in the case of the excitation of the 2s-state 
of the helium ion 


When Z > | we can neglect Z~' in the arguments of Y, and y,,, and M 
is then independent of A. According to (34.10) the cross-section is then 
inversely proportional to Z*. The integration in (34.11) is carried out numeri- 
cally. The angular distributions of the scattcred electrons are anisotropic, 
and when the encrgy of the incident clectron is close to threshold, the 
maximum of the angular distribution lies near 180°. As an example, Fig. 47 
shows the angular distribution for Z = 2 and the Is—2s transition at 
different encrgies. Figure 48 shows the total excitation cross-scction (obtained 
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by numerical integration of the differential cross-section) for the ls—2s 
transition for different Z. The results of calculations using the Born approxi- 
mation and the distorted wave approximation [99] are given for the five-fold 
ionized carbon atom. In the units chosen in Fig. 47 (p = k/Z along the 
abscissa axis) all the thresholds for the excitation of the 2s state lie at the 
single point p? = 0-75. We note that the excitation cross-section at the 
threshold is finite. The reason for this was discussed in § 23. 


k 


rs 


Fic. 48. Excitation cross-sections for the 2s-state in hydrogenic ions of different nuclear 
charge: points—distorted wave approximation (in units of zao*); broken curve—Born 
approximation 


35. OPTICALLY ALLOWED TRANSITIONS ON THE DIPOLE APPROXIMATION 


Let us consider now the excitation of optically allowed transitions. 
If we replace the matrix element Vy,, by its asymptotic expression for 
large r, we have 


; 
VaAAg ~ r3 De (35.1) 


where Daag, 1s the matrix clement for the dipole moment. The quantity 
Daa, iS in its turn related to the oscillator strength for the particular transition. 
The excitation cross-section may be expressed in terms of the oscillator 
strength saq, as follows [98]: 
or AE E 
Ce Og = a 35.2 
ae (47 5) alee & AE an 


where AE is the transition energy, F is the energy of the incident electron; 
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Jo IS given by 


AE E AE erm d 
Se i Vie ee ee 5 
v0( a3 a) i> (NTA a | (in, in’ A, xo)[" (35.3) 
where 7 
1] SS SS 
/2AE 
. Z 
yf = — 
\/2(E— AE) 
din’ 
eee iLL 35.4) 


Xo=-> 
(1y' — 11)” 
and F is the hypergeometric function. 


4 


O 
‘ 2 3 
£ 
AE 
Fic. 49. A plot of go) as a function of E/AE (points represent the distorted wave 
approximation) 


Figure 49 shows go in units of za) as a function of E/AE for different 
AE/Z*. This figure also shows (AE)*¢4,4/84,4 Calculated from the data in 
(99], which gives the calculated cross-section for the transition ls—2p 
in five-fold ionized carbon ions using the distorted wave approximation. 
A detailed table of o9 for different values of the arguments is given in [98]. 


CONCLUSION 


THE aceount presented in this book is in some respeets incomplete. This is 
unavoidable beeause the theory of eleetron-atom eollisions is in the process 
of rapid development, and in many eases only preliminary results are avail- 
able. Nevertheless, it is possible to discern general trends in the development 
of the subject, and this enables us to indicate certain problems which will 
remain of current interest. 

Firstly, we must mention the gencral analysis of cross-seetions near 
threshold with allowance for eoupling between channels and for polari- 
zation of the atom. Studies of resonanees in the excitation eross-seetion for 
a particular atomic state near the excitation threshold for the next state (with 
higher energy) are of particular interest. The M-matrix formalism is econ- 
venient in such investigations (Section 12). Seeondly, we must reeall studies 
of the asymptotic behaviour of the wave function for the p + e + e system 
when both eleetrons are in the continuous spectrum. This is necessary for 
the elucidation of the ionization eross-section near threshold. Caleula- 
tions of eross-seetions for the elastie and inelastic scattering of electrons by 
hydrogen atoms show that partial waves with large L play an appreciable 
part even relatively near threshold. As indicated in {61}, the develop- 
ment of approximate methods of ealeulation whieh do not require the 
partial wave expansions will be of major importanee in this connection. 

In calculations of eross-sections for the scattering of electrons by com- 
pheated atoms, one eneounters the unavoidable problem of the rational 
ehoiee of approximate atomic funetions, and also the problem of the general 
strueture of expressions describing the electron-atom interaction when 
approximate atomie functions are employed (this was considered in Seetion 
26 in connection with the form of the exchange operator). 

Among the current problems in collision theory, there is the dcvelop- 
ment of approximate formulae for estimating the cross-seetions (possibly 
with the aid of some cxperimentally obtained parameters), and also approxi- 
mate formulae which are convenient in the integration of the velocitics of 
the incident electrons over the spectrum. Simple formulae of this nature 
may be obtained on the basis of elassieal mechanics. In this connection we 
must mention the work reported in [100] and [101]. 

Finally, the sueeessful applieation of the methods of quantum ficld 
theory in a number of branehes of thcoretieal physies suggests that the 
introduetion of these methods into the theory of atomie collisions will turn 
out to be fruitful. 
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